NAVAL  POSTGRADUATE  SCHOOL 
Monterey,  California 


THESIS 


ANALYSIS  OF  LOW  PROBABILITY  OF  INTERCEPT  (LPI) 
RADAR  SIGNALS  USING  CYCLOSTATIONARY  PROCESSING 

by 

Antonio  F.  Lima,  Jr. 

September  2002 

Thesis  Advisor: 

Phillip  E.  Pace 

Thesis  Co-Advisor: 

Flerschel  FI.  Loomis 

Approved  for  public  release;  distribution  is  unlimited 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


REPORT  DOCUMENTATION  PAGE 


Form  Approved  OMB  No.  0704-0188 
Public  reporting  burden  for  this  collection  of  information  is  estimated  to  average  1  hour  per  response,  including 
the  time  for  reviewing  instruction,  searching  existing  data  sources,  gathering  and  maintaining  the  data  needed,  and 
completing  and  reviewing  the  collection  of  information.  Send  comments  regarding  this  burden  estimate  or  any 
other  aspect  of  this  collection  of  information,  including  suggestions  for  reducing  this  burden,  to  Washington 
headquarters  Services,  Directorate  for  Information  Operations  and  Reports,  1215  Jefferson  Davis  Highway,  Suite 
1204,  Arlington,  VA  22202-4302,  and  to  the  Office  of  Management  and  Budget,  Paperwork  Reduction  Project 
(0704-0188)  Washington  DC  20503.  _ _ 

I.  AGENCY  USE  ONLY  (Leave  blank)  2.  REPORT  DATE  3.  REPORT  TYPE  AND  DATES  COVERED 

September  2002  Master’s  Thesis 

4.  TITLE  AND  SUBTITLE:  5.  LENDING  NUMBERS 

Analysis  of  Low  Probability  of  Intercept  (LPI)  Radar  Signals  Using  Cyclostationary 
Processing 

6.  AUTHOR(S) 

Antonio  F.  Lima,  Jr. _ 

7.  PERFORMING  ORGANIZATION  NAME(S)  AND  ADDRESS(ES)  8.  PERFORMING 

Naval  Postgraduate  School  ORGANIZATION  REPORT 

Monterey,  CA  93943-5000 _  NUMBER _ 

9.  SPONSORING  /MONITORING  AGENCY  NAME(S)  AND  ADDRESS(ES)  10.  SPONSORING/MONITORING 
Office  of  Naval  Research  AGENCY  REPORT  NUMBER 

II.  SUPPLEMENTARY  NOTES  The  views  expressed  in  this  thesis  are  those  of  the  author  and  do  not  reflect  the  official 
policy  or  position  of  the  Department  of  Defense  or  the  U.S.  Government. 

12a.  DISTRIBUTION  /  AVAll.ABll.lTY  STATEMENT  12b.  DISTRIBUTION  CODE 

Distribution  unlimited. _ 

13.  ABSTRACT  (maximum  200  words) 

LPI  radar  is  a  class  of  radar  systems  possessing  certain  performance  characteristics  that  make  them  nearly 
undetectable  by  today’s  digital  intercept  receivers.  This  presents  a  significant  tactical  problem  in  the  battle  space. 
To  detect  these  types  of  radar,  new  digital  receivers  that  use  sophisticated  signal  processing  techniques  are 
required. 

This  thesis  investigates  the  use  of  cyclo stationary  processing  to  extract  the  modulation  parameters  from  a 
variety  of  continuous-wave  (CW)  low-probability-of-intercept  (LPI)  radar  waveforms.  The  cyclostationary 
detection  techniques  described  exploit  the  fact  that  digital  signals  vary  in  time  with  single  or  multiple  periodicity, 
owing  to  their  spectral  correlation,  namely  non-zero  correlation  between  certain  frequency  components,  at  certain 
frequency  shifts.  The  use  of  cyclo  stationary  signal  processing  in  a  non-cooperative  intercept  receiver  can  help 
identify  the  particular  emitter  and  aid  in  the  development  of  electronic  attack  signals.  LPI  CW  waveforms 
examined  include  Frank  codes,  PI  through  P4,  Frequency  Modulated  CW  (FMCW),  Costas  frequencies  as  well  as 
several  frequency-shift-keying/phase-shift-keying  (FSK/PSK)  waveforms.  This  thesis  show  that  for  signal-to-noise 
ratios  of  0  dB  and  -6  dB,  the  cyclostationary  signal  processing  can  extract  the  modulation  parameters  necessary  in 
order  to  distinguish  between  the  various  types  of  LPI  modulations. 


14.  SUBJECT  TERMS  Low  Probability  of  Intercept  (LPI)  Radars,  Electronic  Support  Measures 
(ESM),  EFT  Accumulation  Method  (FAM),  Direct  Frequency  Smoothing  Method  (DFSM),  Binary 

NUMBER  OF 
PAGES  186 

Phase  Shiit  Keying  (BPSK),  Frequency  Modulated  Continuous  Wave  (FMCW),  Polyphase  Codes  (P4, 
P3,  P2,  PI  and  Frank  Codes),  Combined  FSK/PSK  (Frequency  Shift  Keying  and  Phase  Shift  Keying) 

16.  PRICE  CODE 

17.  SECURITY 

18.  SECURITY 

19.  SECURITY 

20.  LIMITATION 

CLASSIFICATION  OF 

CLASSIFICATION  OF  THIS 

CLASSIFICATION  OF 

OF  ABSTRACT 

REPORT 

PAGE 

ABSTRACT 

Unclassified 

Unclassified 

Unclassified 

UL 

NSN  7540-01-280-5500  Standard  Form  298  (Rev.  2-89) 

Prescribed  by  ANSI  Std.  239-18 


1 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


11 


Approved  for  public  release;  distribution  is  unlimited. 


ANALYSIS  OF  LOW  PROBABILITY  OF  INTERCEPT  (LPI)  RADAR  SIGNALS 
USING  CYCLOSTATIONARY  PROCESSING 

Antonio  F.  Lima,  Jr. 

Captain,  Brazilian  Air  Force 
B.S.,  Brazilian  Air  Force  Academy,  Brazil 


Submitted  in  partial  fulfillment  of  the 
requirements  for  the  degree  of 


MASTER  OF  SCIENCE  IN  SYSTEMS  ENGINEERING 


from  the 


NAVAL  POSTGRADUATE  SCHOOL 
September  2002 


Author:  Antonio  F.  Lima,  Jr. 


Approved  by:  Phillip  E.  Pace 

Thesis  Advisor 


Herschel  H.  Loomis 
Thesis  Co-Advisor 


Dan  C.  Boger 

Chairman,  Information  Sciences  Department 
iii 


THIS  PAGE  INTENTIONALLY  LEET  BLANK 


IV 


ABSTRACT 


LPI  radar  is  a  class  of  radar  systems  that  possess  certain  performanee 
characteristics  that  make  them  nearly  undeteetable  by  today’s  digital  intereept  receivers. 
This  presents  a  significant  tactical  problem  in  the  battle  space.  To  deteet  these  types  of 
radar,  new  digital  receivers  that  use  sophisticated  signal  processing  techniques  are 
required. 

This  thesis  investigates  the  use  of  cyclostationary  processing  to  extract  the 
modulation  parameters  from  a  variety  of  continuous-wave  (CW)  low-probability-of- 
intercept  (LPI)  radar  waveforms.  The  eyclostationary  detection  techniques  deseribed 
exploit  the  fact  that  digital  signals  vary  in  time  with  single  or  multiple  periodicities, 
beeause  they  have  speetral  eorrelation,  namely,  non-zero  correlation  between  certain 
frequeney  components,  at  certain  frequency  shifts.  The  use  of  cyclostationary  signal 
processing  in  a  non-cooperative  intercept  receiver  can  help  identify  the  partieular  emitter 
and  can  help  develop  eleetronic  attaeks.  LPI  CW  waveforms  examined  inelude  Frank 
eodes,  polyphase  codes  (PI  through  P4),  Frequeney  Modulated  CW  (FMCW),  Costas 
frequeneies  as  well  as  several  frequency-shift-keying/phase-shift-keying  (FSK/PSK) 
waveforms.  It  is  shown  that  for  signal-to-noise  ratios  of  OdB  and  -6  dB,  the 
cyclostationary  signal  processing  can  extract  the  modulation  parameters  necessary  in 
order  to  distinguish  among  the  various  types  of  LPI  modulations. 
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EXECUTIVE  SUMMARY 


LPI  radar  is  a  class  of  radar  systems  possessing  eertain  performance 
eharaeteristics  that  make  today’s  digital  intercept  reeeivers  virtually  unable  to  deteet 
them.  This  presents  a  significant  tactical  problem  in  the  battle  space.  To  detect  these 
types  of  radar,  the  military  requires  new  digital  receivers  that  use  sophisticated  signal 
processing  techniques. 

This  thesis  investigates  the  use  of  eyclostationary  processing  to  extraet  the 
modulation  parameters  from  a  variety  of  eontinuous-wave  (CW)  low-probability-of- 
intereept  (LPI)  radar  waveforms.  The  eyclostationary  detection  techniques  described 
exploit  the  faet  that  digital  signals  vary  in  time  with  single  or  multiple  periodieity,  owing 
to  their  speetral  correlation,  namely,  non-zero  correlation  between  certain  frequency 
components,  at  eertain  frequeney  shifts.  The  use  of  eyelostationary  signal  processing  in  a 
non-eooperative  intereept  reeeiver  ean  help  identify  the  partieular  emitter  and  can  help 
develop  eleetronic  attaeks.  LPI  CW  waveforms  examined  include  Frank  codes,  polyphase 
eodes  (PI,  P2,  P3  and  P4),  Frequeney  Modulated  CW  (FMCW),  Costas  frequeneies  as 
well  as  several  frequeney-shift-keying/phase-shift-keying  (FSK/PSK)  waveforms.  This 
thesis  shows  that  for  signal-to-noise  ratios  of  OdB  and  -6  dB,  the  eyelostationary  signal 
processing  can  extract  modulation  parameters  such  as  carrier  frequency,  chip  rate,  code 
rate  and  bandwidth,  neeessary  in  order  to  distinguish  between  the  various  types  of  LPI 
radar  modulations. 

Two  computationally  efficient  methods  of  eyelostationary  processing  were 
implemented:  Time  Smoothing  FFT  Aeeumulation  and  Direet  Frequeney  Smoothing.  It  is 
possible  to  verify  that  the  time  smoothing  method  is  more  eomputationally  efficient  than 
the  frequeney  smoothing  for  signals  with  higher  eomplexity  (polyphase  eodes,  Frank 
codes,  Costas  and  FSK/PSK).  The  results  from  both  methods  were  compared  and 
discussed  for  various  LPI  modulation  types. 

The  results  due  to  variations  of  modulation  characteristics  are  compared  and  the 
effieiency  of  both  eyelostationary  methods  for  each  modulation  is  measured  with  relation 
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to  the  original  parameters.  This  thesis  also  ineludes  eomments  on  whieh  LPI  radar  signals 
were  more  suitable  for  eyelostationary  analysis  and  suggestions  for  future  elassifieation 
systems  for  these  signals,  using  eombined  teehniques. 
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I.  INTRODUCTION 


A,  LPI  RADARS 


Standard  surveillance  radars  always  faced  the  problem  of  being  detected  by  third- 
party  intercept  receivers.  To  mitigate  this  problem,  a  low  probability  of  interception  has 
become  a  very  common  tactical  requirement  for  new  radars.  The  new  task  of  these  radars 
is  essentially  two  fold;  the  radar  must  meet  its  specified  detection  performance  with 
minimum  radiated  power  and  must  measure  target  characteristics  with  a  modulation  that 
is  difficult  for  an  intercept  receiver  to  identify.  [1] 

Because  of  their  low  power,  wide  BW  (BW)  and  frequency  variability,  LPI  radars 
make  detection  by  modern  intercept  receivers  difficult.  The  resulting  basic  definition  is 
that  LPI  radars  can  detect  targets  at  longer  ranges  than  the  modem  intercept  receivers 
detect  the  radars. 


This  main  characteristic  of  LPI  radars  may  be  described  by  the  sensitivity 
advantage  (S),  defined  as  the  ratio  between  the  signal  power  needed  at  the  intercept 
receiver  to  that  needed  at  the  LPI  radar,  which  can  be  expressed  mathematically  as  [2]; 


P^T  V  y  V  y 


(1.1.1) 


where  Pm  is  the  power  required  at  the  intercept  receiver  to  detect  a  signal;  Prt  is  the 
power  required  at  the  LPI  radar  receiver  to  detect  a  target;  <jj  is  the  target’s  RCS  (Radar 

Cross  Section);  Gt  is  the  bore  sight  gain  of  the  LPI  radar’s  directive  transmit  antenna;  G, 

is  the  gain  of  the  LPI  radar’s  transmit  antenna  side  lobe  in  the  direction  of  the  intercept 
receiver;  Gi  is  the  gain  of  the  intercept  receiver’s  antenna;  Gr  is  the  radar’s  receive 
antenna;  Ljr  is  the  loss  in  the  intercept  receiver;  Lrr  is  the  loss  between  the  radar’s 

antenna  and  receiver;  Rj.  is  the  radar  to  target  range  and  Rj  is  the  radar  to  intercept 

receiver  range.  The  sensitivity  advantage  ( (5' )  depends  on  the  intercept  receiver 
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characteristics  and  should  be  a  high  value,  on  the  order  of  50  dB,  for  a  ease  where  we 
have  a  simple  reeeiver  against  an  LPI  radar. 

The  sueeess  of  LPI  radars  depends  on  how  hard  it  is  for  a  reeeiver  to  deteet  the 
radars’  emission  parameters.  The  proeessing  eapabilities  of  modern  ES  (Eleetronie 
Support)  equipment  are  inereasing,  leading  to  more  speeifie  EPI  requirements.  On  the 
reeeiver  side,  better  results  in  speetral  analysis,  for  non-eooperative  deteetion  and 
elassifieation,  may  be  obtained  if  these  radar  signals  are  modeled  as  eyelostationary. 

All  digitally  modulated  signals  are  eyelostationary,  meaning  that  their 
probabilistie  parameters  vary  in  time  with  single  or  multiple  periodieities.  One  property 
that  extends  from  this  is  that  digitally  modulated  signals  have  a  eertain  amount  of  speetral 
eorrelation.  In  other  words,  the  signal  is  eorrelated  with  frequeney-shifted  versions  of 
itself  (auto-eorrelation)  at  eertain  frequeney  shifts.  Analyzing  LPI  radar  waveforms  using 
eyelostationary  modeling  is  advantageous  beeause  non-zero  eorrelation  is  exhibited 
between  eertain  frequeney  eomponents  when  their  separation  is  related  to  the  periodieity 
of  interest  (e.g.,  the  symbol  rate  or  earner  frequeney).  The  value  of  that  speetral 
separation  is  referred  to  as  the  cycle  frequency. 

Two  main  algorithms  stand  out  as  eomputationally  effieient  tools  for 
eyelostationary  signal  proeessing.  The  first  is  the  time  smoothing  Past  Pourier  Transform 
(PET)  Aeeumulation  Method  and  the  other  is  the  Direet  Prequeney  Smoothing  Method 
[3].  Both  tools  are  implemented  in  MATLAB®  6.1  for  this  thesis.  [4] 

B,  PRINCIPAL  CONTRIBUTIONS 

The  objeetive  of  the  researeh  deseribed  in  this  thesis  was  to  implement  in 
MATEAB®  [4],  two  eomputationally  effieient  eyelostationary  algorithms  known  as  the 
Time  Smoothing  PET  Aeeumulation  Method  and  the  Direet  Prequeney  Smoothing 
Method,  defined  in  [3]  and  investigate  them  as  an  ES  reeeiver  for  proeessing  EPI  radar 
signals. 

The  first  step  was  to  generate  the  EPI  signals  in  a  standardized  way.  The  eode 
used  was  ealled  the  “EPI  Signal  Generator,”  developed  by  Pernando  Taboada  [5],  and 
also  ineludes  eontributions  by  the  author  of  this  thesis.  The  generated  modulations  were 
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•  BPSK  (Binary  Phase  Shift  Keying); 

•  FMCW  (Frequeney  Modulated  Continuous  Wave); 

•  Polyphase  Codes  (P4,  P3,  P2  PI  and  Frank  Codes); 

•  Costas  Codes  (Frequency  Hopping  -  FH); 

•  FSK/PSK  (Combined  Frequency  Shift  Keying  and  Phase  Shift  Keying)  with  a 
Costas  frequency  distribution;  and 

•  FSK/PSK  with  a  target  matched  frequency  distribution. 

Once  the  signal  test  matrix  was  completed,  simulations  to  verify  the  implementation  of 
each  algorithm  were  performed  in  MATLAB®  [4]  and  the  results  were  compared  with 
other  receiver  signal  processing  techniques,  such  as  Higher  Order  Statistics  [5], 
Quadrature  Mirror  Filter  Banks  [6]  and  the  Wigner  Distribution  [7]. 

A  Graphic  User  Interface  (GUI)  was  developed  in  MATLAB®  [4]  to  simplify  the 
analysis  of  the  simulation  results.  The  output  obtained  from  the  cyclostationary  signal 
models  were  then  used  to  determine  the  various  characteristics  of  each  modulation  in 
question. 

Previous  work  has  been  done  to  analyze  phase  modulation  techniques  such  as 
BPSK  (Binary  Phase-Shift  Keying)  and  QPSK  (Quaternary  Phase-Shift  Keying)  using 
time-smoothing  techniques  [8,  9].  In  this  thesis,  both  frequency  and  time-smoothing 
techniques  are  used  to  analyze  various  LPI  radar  modulations  and  to  evaluate  the 
measurement  of  the  modulation  parameters. 

C.  THESIS  OUTLINES 

The  purpose  of  this  thesis  is  to  document  the  software  implementation  of  a  non- 
cooperative  cyclostationary  receiver  for  LPI  radar  waveforms.  The  remainder  of  this 
thesis  supports  this  purpose  and  is  organized  as  follows. 

Chapter  II  presents  a  brief  description  of  Low  Probability  of  Intercept  (LPI) 
waveforms  and  their  spectral  properties.  Two  FSK  and  PSK-combined  modulations  are 
discussed  and  analyzed.  BPSK,  P4,  FMCW  and  Costas  Codes  are  described  in  depth  in 
[5]. 
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Chapter  III  presents  the  Cyclostationary  signal  processing  algorithms,  a  brief 
description  of  the  cyclostationary  processor,  the  MATLAB®  [4]  tools  and  the  extracted 
parameters  description. 

Chapter  IV  shows  the  analysis  of  the  different  modulation  types,  their  parameters, 
as  well  as  the  simulation  results. 

Chapter  V  summarizes  the  results  of  this  thesis  and  also  makes  recommendations 
for  future  research. 

Appendix  A  contains  the  MATLAB®  [4]  M-files  used  for  implementation  of  both 
algorithms  from  [3]. 

Appendix  B  contains  the  MATLAB®  [4]  M-files  for  the  LPI  Generator  blocks  for 
FSK/PSK  Costas  and  Target  signals. 

Appendix  C  contains  a  table  of  all  LPI  radar  signals  analyzed. 
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II.  LPI  WAVEFORMS  DESCRIPTION 


A,  BACKGROUND 

LPI  radars  are  especially  designed  to  oppose  external  third-party  receivers  that 
attempt  to  identity  the  system  characteristics  and  emitter  location.  Modern  electronic 
support  (ES)  receivers  can  easily  detect  the  high  peak  power  transmitted  by  pulsed  radars. 
The  use  of  CW  (Continuous  Wave)  modulations  and  the  ability  to  manage  the  transmitted 
power  limiting  emission  to  the  minimum  power  required  to  detect  typical  targets,  at  the 
required  range,  make  LPI  radar  signals  much  less  detectable.  [5]  Besides  power 
management,  LPI  radars  modulate  their  transmissions  spreading  the  energy  in  frequency 
so  that  the  frequency  spectrum  of  the  transmitted  signal  is  wider  than  required  to  carry 
the  signal’s  information.  Spreading  the  signal  energy  reduces  the  signal-strength-per- 
information  BW. 

LPI  waveforms  investigated  in  this  thesis  include  BPSK,  LMCW,  P4,  P3,  P2,  PI, 
Prank  Codes,  Costas  Codes,  PSK/PSK  with  a  Costas  frequency  distribution,  and 
PSK/PSK  with  a  target  matched  frequency  distribution.  Refer  to  Pernando  Taboada’s  [5] 
thesis  work  for  a  detailed  description  of  the  other  LPI  modulations  analyzed  in  this  thesis. 
The  complete  matrix  of  analyzed  signals  is  shown  in  Appendix  C.  This  thesis  presents  the 
analysis  of  one  signal  example  per  modulation  type.  The  analysis  of  the  rest  of  the  signals 
is  included  in  a  Technical  Report  to  be  published.  [10] 

This  chapter  specifically  discusses  the  two  modulation  types  that  combine 
frequency  and  phase-shift  keying  (PSK/PSK  combined).  One  modulation  type  is  a 
combination  of  a  Costas  frequency-hopping  technique  and  binary  phase  modulation  using 
Barker  sequences  of  different  lengths.  The  second  is  a  frequency-hopping  technique  that 
uses  the  characteristic  frequency  distribution  of  a  desired  target,  creating  a  matched  PSK, 
which  is  then  modulated  with  a  random-phase  keying. 
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B,  FSK/PSK  COMBINED  USING  A  COSTAS-BASED  FREQUENCY¬ 
HOPPING  (FH)  TECHNIQUE 

This  modulation  technique  is  the  result  of  a  combination  of  frequency-shift 
keying  based  on  a  Costas  frequency-hopping  matrix  and  phase-shift  keying  using  Barker 
sequences  of  different  lengths.  A  thorough  description  of  the  implementation  of  a  Costas 
frequency-hopping  technique  is  in  [5].  The  purpose  of  this  section  is  to  describe  briefly 
the  phase  encoding  applied  to  a  Costas  signal,  generating  the  FSK/PSK  combined 
waveform. 


In  a  Costas  frequency-hopped  signal,  the  firing  order  defines  what  frequencies 
will  appear  and  with  what  duration.  Since  we  are  discussing  CW  radars,  the  usual 
terminology  does  not  apply  to  this  case.  Instead  of  a  “bursf’  of  pulses,  we  have 
frequencies  being  continuously  emitted  during  a  defined  period  of  time.  This  period  may 
be  divided  into  sub-periods,  labeled  TV  for  each  frequency.  The  length  of  each  sub-period 
depends  on  the  sampling  interval.  During  each  sub-period,  the  signal  frequency  (one  of 
the  Np  frequencies)  is  modulated  by  a  binary  phase  sequence  according  to  a  Barker 


sequence  of  length  five  (-h  -h  -h  -  +),  seven  {+  +  +  --  +  -),  eleven  (+  +  + - +  --  -!--)  or 

thirteen  For  example,  the  FSK/PSK  signal  defined  by  S  =  1"^, 

1^  1^,  r,  1^  2+,  2^  2+,  2',  2^,  3^  3^  3^  3',  3^,  4^  4+,  4^  4',  4+  5^  5^  5^  5',  5^ 
represents  a  waveform  comprised  of  Np  =5  different  frequencies,  that  are  each 


subdivided  into  five  phase  slots,  labeled  TV,  according  to  the  Barker  sequence  of  length 
five  (+  -h  -h  -  -h).  The  final  waveform  may  be  seen  as  a  binary  phase-shifting  modulation 
within  each  frequency  hop,  resulting  in  5  phase  slots  equally  distributed  in  each 
frequency  slot,  giving  a  total  of  25  phase  slots. 


As  illustrated  in  Figure  1,  if  we  consider  Np  as  the  number  of  frequency  hops  and 
Np  as  the  number  of  phase  slots  of  duration  TV  (Chip  Period)  in  each  frequency  sub¬ 
period  TV,  the  total  number  of  phase  slots  in  the  FSK/PSK  waveform  is  given  by: 


N  =  Np  *  Np 


(2.1.1) 
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The  block  diagram  in  Figure  1  describes  the  MATLAB®  [4]  implementation.  The 
user  defines  which  sequence  of  Costas  frequency  hops  to  be  used  and  also  how  long  the 
Barker  sequence  is  (5,  7,  11  or  13).  The  number  of  frequency  hops  is  pre-defmed  to  be 
“seven”  and  the  user  may  select  from  two  different  frequency  sequences,  varying  from 
IkHz  to  8kHz.  The  Costas  matrix  used  in  the  implementation  is  the  following; 

•  Costas  Sequence  1— >  4716523 

•  Costas  Sequence  2— >  2  63875  1. 

Frequencies  (kHz) 


Voltage 


+1. 

1 

* —  Frequency  1  — ► 

■* —  Frequency  2  — ► 

■* —  Frequency  5  — ► 

Time 

0- 

-1 

+ 

+ 

+ 

- 

+ 

2 

Tp... 

NpTp=  Tp  2NpTp=2Tp  (NprOTp  NpTp 

b) 

Figure  1  a)  FSK/PSK  Costas  LPI  Generator  MATLAB®  [4]  code  block 

diagram  and,  b)  general  FSK/PSK  signal  containing  Nf  frequency 
hops  with  Np  phase  slots  per  frequency. 

The  Barker  sequence  is  generated  and  the  frequency-hopping  signal  is  then  phase- 
modulated  accordingly.  For  example,  if  the  first  Costas  sequence  is  selected,  after  a  phase 
modulation  using  a  Barker  sequence  of  length  5,  the  final  waveform  becomes  S  =  4^,  4^, 
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2^,  3"^,  3^,  3"^,  3',  3^.  Figure  2  shows  the  Power  Spectral  Density  (PSD)  plots  that  reveal 
the  spread  spectrum  characteristic  of  these  signals.  The  Costas  sequence  used  in  the 
following  example  is  always  the  same  and  the  seven  frequency  hops  are  4,  7,  1,  6,  5,  2, 
and  3kHz  .  The  sampling  frequency  was  15kHz  ,  satisfying  the  minimum  Nyquist  rate 
(fs»2*f)  for  the  largest  frequency  value.  All  plots  were  generated  in  MATLAB®  [4] 
using  the  routines  “fsk_psk_costas.m”  and  “PAF_FSK_PSK.m”;  both  listed  in  Appendix 
B. 


Figure  2  PSD  plot  for  a  Costas  FH  waveform  with  no  phase  modulation. 


Figure  3  shows  the  Costas  frequency-hopping  waveform  time  domain,  before  the 
phase  modulation. 


Time  -  Seconds 

Figure  3  Time  domain  plot  for  a  Costas  FH  waveform  with  no  phase 
modulation. 

Figures  4  to  11  show  the  PSD,  Periodic  Ambiguity  Function  (PAF)  and  Phase 
Plots  for  the  same  frequency  sequence  as  in  Figure  2  but  now  with  a  phase  modulation 
using  a  Barker- 11  sequence.  Other  signals  that  may  be  generated  are  Costas  Frequency- 
Hopping  with  Barker-5,  Barker-7  and  Barker- 13  phase  modulation.  These  signals  are 
going  to  be  analyzed  later  in  Chapter  IV. 

The  firing  orders  generated  using  Costas  arrays  are  designated  as  optimum  in 
reference  to  the  side-lobe  behavior  of  the  PAF.  As  we  increase  the  length  of  the  Barker 
sequence  (e.g.  from  five  to  eleven),  we  notice  a  decrease  in  the  side-lobe  level  as  well. 

The  waveforms  generated  are  just  examples  and  various  combinations  of 
frequency-shift  keying  and  phase-shift  keying  may  be  applied  to  obtain  similar  results. 
Other  examples  of  these  kinds  of  waveforms  are  presented  in  [12]. 
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The  PAF  plots  were  performed  both  for  a  eomplete  period  for  all  Costas 
frequencies  in  the  sequence  and  for  only  one  period  in  one  frequency  hop.  Therefore,  the 
Costas  PAF  characteristics  as  well  as  the  BPSK  PAF  characteristics  may  be  compared. 

Figure  4  shows  the  PSD  of  the  Costas  sequence  1  but  now  with  the  phase 
modulation  using  a  Barker- 11  sequence  and  with  one  carrier  cycle  per  phase  {cpp).  We 
may  verify  the  “flat  top”  characteristic  of  the  spectrum  of  LPI  signals  and  also  the  seven 
“frequency  hops.” 


Figure  4  PSD  for  FSK/PSK  Costas  FH  phase  modulated  with  a  Barker- 1 1 
sequence  and  1  cpp. 

The  phase  plot  in  Figure  5  reveals  the  Barker- 11  sequence  phase  change.  Figure 
6,  Figure  7  and  Figure  8  show  the  PAF  “thumbtack”  [11]  characteristic  of  these  types  of 
Costas  signals.  These  PAF  plots  were  generated  for  one  period  of  the  whole  Costas 
sequence. 
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Figure  5  Phase  plot  for  FSK/PSK  Costas  FH  phase  modulated  with  a  Barker- 
1 1  sequence. 
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Figure  6  PAF  contour  plot  for  FSK/PSK  Costas  FFl  phase  modulated  with  a 
Barker- 1 1  sequence  (plot  of  one  period  for  all  frequencies  in  one 
Costas  sequence). 
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Figure  7  PAF  delay  axis  cut  for  FSK/PSK  Costas  FFI  phase  modulated  with  a 
Barker- 1 1  sequence  (plot  of  one  period  for  all  frequencies  in  one 
Costas  sequence). 


Normalized  v 

Figure  8  PAF  Doppler  axis  cut  for  FSK/PSK  Costas  FH  phase  modulated  with 
a  Barker- 1 1  sequence  (plot  of  one  period  for  all  frequencies  in  one 
Costas  sequence). 
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Figure  9,  Figure  10  and  Figure  11  show  the  modulation  for  one  period  of  one 
frequency  hop  of  the  Costas  sequence  used  in  Figure  2.  Each  frequency  hop  has  similar 
PAF  characteristics  of  a  single  carrier  frequency  BPSK  signal.  The  plots  of  these  Figures 
were  done  for  the  first  carrier  frequency  of  the  first  Costas  sequence.  Other  modulation 
examples  are  analyzed  later  in  Chapter  IV. 


Figure  9  PAF  contour  plot  for  FSK/PSK  Costas  FH  phase  modulated  with  a 
Barker- 1 1  sequence  (plot  of  one  period  for  one  frequency  in  the 
Costas  sequence). 
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Figure  10  PAF  delay  axis  cut  for  FSK/PSK  Costas  FH  phase  modulated  with  a 
Barker- 1 1  sequence  (plot  of  one  period  for  one  frequency  in  the 
Costas  sequence). 


Normalized  v 


Figure  1 1  PAF  Doppler  axis  cut  for  FSK/PSK  Costas  FH  phase  modulated  with 
a  Barker- 1 1  sequence  (plot  of  one  period  for  one  frequency  in  the 
Costas  sequence). 
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C.  FSK/  PSK  COMBINED  USING  A  TARGET-MATCHED  FREQUENCY 

HOPPING 

Instead  of  spreading  the  energy  of  the  signal  equally  over  a  broad  BW,  this  type 
of  technique  concentrates  the  signal  energy  in  specific  spectral  locations  of  most 
importance  for  the  radar  and  its  typical  targets.  The  signals  have  a  pulse  compression 
characteristic,  and  therefore  they  can  achieve  a  low  probability  of  intercept. 

The  implementation  starts  with  a  simulated-target  time-radar  response.  The  block 
diagram  in  Figure  12  describes  the  signal  generation  in  detail.  The  target  signature  data  is 
Fourier  transformed  and  the  frequency  components,  their  correspondent  magnitudes,  and 
their  initial  phases  are  collected.  A  random  selection  process  chooses  each  frequency 
with  a  probability  distribution  function  defined  by  the  spectral  characteristics  of  the  target 
of  interest  obtained  from  the  Fast  Fourier  Transform  (FFT).  That  is,  the  frequencies  at  the 
spectral  peaks  of  the  target  (highest  magnitudes)  are  transmitted  more  often.  Each 
“frequency  hop,”  transmitted  during  a  specific  period  of  time,  is  also  modulated  in  phase, 
having  its  initial  phase  value  (from  FFT)  modified  by  a  pseudo-random  spreading-phase 
sequence  code  of  values  equally  likely  to  be  zero  or  n  radians.  [11]  The  matched 
FSK/PSK  radar  will  then  use  a  correlation  receiver  with  a  phase  mismatched  reference 
signal  instead  of  a  perfectly  phase  matched  reference.  This  allows  the  radar  to  generate 
signals  that  can  match  a  target’s  spectral  response  in  both  magnitude  and  phase.  [11] 
Only  a  single  target  test  signal  is  generated  and  serves  our  purpose  of  testing  the 
performance  of  a  digital  cyclostationary  receiver  against  these  kinds  of  signals. 
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Generate  Random 

Frequency  Firing  Order 

Random  Phase 
Generation 

I  =  cos(2rfinT  +  O  + 

Q  =  sin(2rfinT  -t  0  +  where  i  =  1  to  256 

Figure  12  Block  diagram  of  the  MATLAB®  [4]  implementation  of  FSK/PSK 
target  matched  waveform. 


Figure  13  shows  the  64  complex  points  target  range  radar  response  plot.  Figure  14 
reveals  the  64  frequency  components  that  will  be  selected  randomly  256  times.  Figure  15 
illustrates  the  frequency  firing  order,  Figure  16  illustrates  the  PSD  and  Figures  17,  18  and 
19  illustrate  the  PAF  properties  of  these  signals.  Figure  15  shows  the  histogram  of  the  64 
frequency  components  and  shows  the  number  of  occurrences  of  each  frequency.  Note 
that  this  is  similar  to  the  FFT  output  or  probability  distribution  shown  in  Figure  14.  The 
following  figures  show  one  signal  example  with  5  carrier  periods  per  phase  and  256 
frequency  hops.  Figure  16  shows  the  PSD  plots  and  reveals  the  highly  spread-spectrum 
characteristics  of  this  type  of  modulation.  Note  the  noise-like  behavior  due  to  the  random 
phase  modulation. 
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Figure  13  FSK/PSK  target  64  complex  points  radar  range  simulated  response. 


Frequency 

Figure  14  FSK/PSK  target  frequency  probability  distribution  of  64  frequency 
components. 
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Frequency  (Hz) 

Figure  15  FSK/PSK  target  64  frequency  components  histogram  with  number  of 
occurrences  per  frequency  for  256  frequency  hops. 


Figure  16  PSD  for  FSK/PSK  target  with  64  frequency  components,  256 
frequency  hops,  random  phase  modulation  and  5  cpp. 
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Figure  17  PAF  contour  plot  for  FSK/PSK  target  with  64  frequency  components, 
256  frequency  hops,  random  phase  modulation  and  5  cpp. 


- 1 - 1 - 1 

- 1 - 

- 1 - 1 - 1 - 1 - 

-0.8  -0.6  -0.4  -0.2  0  0.2  04  0.6  0.8 

Normalized  t 


Figure  18  PAF  delay  cut  for  FSK/PSK  target  with  64  frequency  components, 
256  frequency  hops,  random  phase  modulation  and  5  cpp. 
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Figure  19  PAF  Doppler  cut  for  FSK/PSK  target  with  64  frequency  components, 
256  frequency  hops,  random  phase  modulation  and  5  cpp. 


An  extensive  discussion  regarding  the  PAF  “thumbtack”  characteristic  of  these  types  of 
waveforms,  as  shown  in  the  last  three  Figures,  are  presented  in  Donohoe  et  al  in  [1 1], 

In  this  chapter  we  discussed  the  implementation  of  two  complex  LPl  Radar 
signals  using  FSK  and  PSK  techniques  combined.  A  brief  theoretical  and  practical 
tutorial  on  cyclostationarity  processing  and  its  implementation  using  the  FFT 
Accumulation  Method  and  the  Direct  Frequency  Smoothing  Method  is  given  in  the  next 
chapter. 
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III.  CYCLOSTATIONARY  SIGNAL  PROCESSING 
ALGORITHMS  AND  TUTORIAL 


This  chapter  briefly  explains  the  eyelostationary  proeesses,  the  time-smoothing 
(FAM)  and  frequeney-smoothing  (DFSM)  algorithms  and  how  they  were  implemented. 
A  thorough  deseription  on  eyelostationarity  and  its  properties  may  be  found  in  [12],  [13] 
and  [14], 


A,  CYCLOSTATIONARY  THEORY 

The  eyelostationary  theory  for  signal  proeessing,  as  deseribed  by  William  A. 
Gardner  [12],  involves  three  main  properties: 

•  Generation  of  speetral  lines  by  quadratieally  transforming  a  signal; 

•  The  statistieal  property  ealled  “second-order  eyelostationarity,”  namely  the 
periodie  fluetuation  of  the  auto-eorrelation  funetion  with  time; 

•  The  eorrelation  property  for  signal  eomponents  in  distinet  speetral  bands. 

The  eyelostationary  attribute,  as  it  is  refleeted  in  the  periodieities  of  the  seeond 
order  moments  of  the  signal,  ean  be  interpreted  in  terms  of  the  generation  of  speetral  lines 
from  the  signal  by  putting  the  signal  through  a  quadratie  non-linear  transformation.  This 
property  explains  the  link  between  the  speetral-line  generation  property  and  the  statistieal 
property  ealled  “speetral  eorrelation”,  eorresponding  to  the  eorrelation  that  exists 
between  the  random  fluetuations  of  eomponents  of  the  signal  residing  in  distinet  speetral 
bands.  The  eorrelation  integral  is  very  important  in  theoretieal  and  praetieal  applieations 
and  may  be  defined  as 

00 

h{x)=  ^  f{u)g{x  +  u)du  (3.1.1) 

-00 

Applying  an  FFT,  it  forms  a  Fourier  transform  pair  given  by: 

3{/2(x)}=F(5)G*(5)  (3.1.2) 
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If  f(x)  and  g(x)  are  the  same  funetion,  the  integral  above  is  normally  ealled  the 
autocorrelation  funetion  and  ealled  cross-correlation  if  they  differ.  The  autoeorrelation 
funetion  is  a  quadratie  transformation  of  a  signal  and  may  be  interpreted  as  a  measure  of 
the  predietability  of  the  signal  at  time  t  +t  based  on  knowledge  of  the  signal  at  time  t. 
[13] 

When  eonsidering  a  time  series  of  length  T,  the  autoeorrelation  funetion  beeomes 
the  time-average  autoeorrelation  funetion  given  by 


(3.1.3) 


The  non-zero  eorrelation  (seeond-order  periodieity)  eharaeteristie  of  a  time  series  x(t) 
exists,  in  the  time  domain,  if  the  equation. 


^0 


(3.1.4) 


where  a  is  the  eyele  frequeney.  R“  (r)  is  the  eyelie  auto-eorrelation  funetion,  also  known 
as  the  “time-frequeney  limit  autoeorrelation  funetion”.  The  derivation  of  (3.1.4)  from 
(3.1.3)  using  a  non-probabilistie  approaeh  is  developed  in  [10].  Sinee  (3.1.4)  is  a 
generalization  of  (3.1.3),  when  «  =  0,  the  DC  eomponent  of  (3.1.4)  yields  the  time- 
average  autoeorrelation  funetion  of  (3.1.3).  Therefore,  the  proeess  defined  by  (3.1.4)  is 
able  to  extraet  more  information  from  the  signal  than  the  proeess  defined  by  (3.1.3).  [13] 

It  is  well  known  that  the  PSD  may  be  obtained  from  the  Fourier  Transform  of  the 
autoeorrelation  funetion  (3.1.3).  [14] 

00 

SAf)=\KO)e-‘^’'’dT  (3.1.5) 


In  the  same  manner,  it  is  shown  in  [3]  that  the  Speetral-Correlation  Density  (SCD),  or 
Cyelie-Speetral  Density,  may  also  be  obtained  from  the  Fourier  Transform  of  the  eyelie 
autoeorrelation  funetion  (3.1.4) 
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(3.1.6) 


=  lim  -  X  J/  +  -  V;  f  /  - 

2)  2) 

where  a  is  the  cycle  frequency  and: 

T 

X,{f)=\x{u)e-^^’^^^du  (3.1.7) 

-T 

which  is  the  Fourier  Transform  of  the  time  domain  signal  x(u).  The  additional  variable 
a  leads  to  a  two-dimensional  representation  S"  (/) ,  in  the  bi-frequency  plane  ox  if,  a  ) 
plane.  [12] 

The  speetral  correlation  exhibited  by  eyelostationary  or  almost  cyclostationary 
processes  is  completely  charaeterized  by  the  eyelie  speetra  (5*“)  or  charaeterized 
equivalently  by  the  cyclic  autocorrelations  (7?“).  [12]  In  praetiee,  the  cyclie-spectral 

density  must  be  estimated  beeause  the  signals  being  eonsidered  are  defined  over  a  finite 
time  interval  ( A  f),  and  therefore  the  cyelic-spectral  density  eannot  be  measured  exactly. 
Estimates  of  the  cyelie-spectral  density  ean  be  obtained  via  time  or  frequency-smoothing 
techniques.  In  this  work  we  will  be  able  to  compare  both  methods  when  analyzing  LPl 
radar  signals. 

An  estimate  of  the  SCD  can  be  obtained  by  the  time-smoothed  eyelie 
periodogram  is  given  by  [10]; 

At 
t-\ - 

j  s^juj)du  (3.1.8) 

At 

t - 

2 

where 


I  R-;(T)e 

—X) 


A  (u,f)  =  —X^  \u,f  +  -]x*  \uj-- 


(3.1.9) 


and  At  =  total  observation  time  of  the  signal,  Tw=  short-time  EFT  window  length,  and; 
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(3.1.10) 


Xj,^(m,/)=  I  x{u)e 


2 


is  the  sliding  short-time  Fourier  Transform,  and  is  a  viable  solution  for  computing  the 
SCD  estimations.  Using  a  graphical  explanation,  in  Figure  20,  for  some  signal  x(t)  the 
frequency  components  are  evaluated  over  a  small  time  window  (sliding  FFT  time 
length),  along  the  entire  observation  time  interval  At.  [12]  The  spectral  components 
generated  by  each  short-time  Fourier  Transform  have  a  resolution.  A/  =  .In  Figure 


20,  L  is  the  overlapping  factor  between  each  short-time  FFT.  In  order  to  avoid  aliasing 


and  cycle  leakage  on  the  estimates,  the  value  of  L  is  defined  as  F  < 


k 


At 


i  L  \ 

I _ I 


Short  FFT’s 


J 


J 


Figure  20  Pictorial  illustration  of  the  estimation  of  the  time-variant  spectral 
periodogram  (adapted  from  [12,  17]). 
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Figure  21  shows  that  the  spectral  components  of  each  short-time  FFT  are 
multiplied,  still  providing  the  same  resolution  capability  A/  =  jZ  ,  for  the  cyclic- 

/  -'(T 

spectrum  estimates.  Note  that  the  dummy  variable  u  has  been  replaced  by  the  time 
instances  At  each  window  {Tw),  two  components  centered  about  some  fo  and 

separated  by  some  are  multiplied  together  and  the  resulting  sequence  of  products  is 
then  integrated  over  the  total  time  (At),  as  shown  in  (3.1.8). 


Figure  21  Sequence  of  frequency  products  for  each  short-time  Fourier 
transforms  (adapted  from  [12,  17]). 


The  estimation  5*“  (/) «  5*“  (t,  /)^,  can  be  made  as  reliable  and  accurate  as 

desired  for  any  given  t  and  Af,  and  for  all  f  by  making  At  sufficiently  large,  provided  that 
equation  (3.1.4)  exists  within  the  interval  At  and  that  a  substantial  amount  of  smoothing 
is  carried  out  over  At,  which  leads  to  the  Grenander’s  Uncertainty  Condition  At*  A/"  »  1 
[12].  This  Uncertainty  Condition  means  that  the  observation  time  (At)  must  greatly 
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exceed  the  time  window  {Tw),  which  is  used  to  compute  the  spectral  components.  A  data 
taper  window  is  also  used  to  minimize  the  effects  of  cycle  and  spectral  leakage 
(estimation  noise),  introduced  by  frequency  component  side-lobes.  [12] 


If  we  consider  the  fact  that  the  cycle  frequency  estimate  is  Aa  w  ,  it  results 

that  the  estimation  of  some  (^,  )  represents  a  very  small  area  on  the  bi-frequency  plane 

as  shown  in  Figure  22  and  since  one  needs  a  significant  number  of  estimates  to  represent 
the  cyclic  spectrum  adequately,  it  follows  that  obtaining  estimates  may  become  very 
computationally  demanding.  [12] 


Figure  22  Bi-frequency  plane,  frequency  and  cycle  frequency  resolutions  on 
detailed  area  (adapted  from  [12,  17]). 
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B, 


DISCRETE  TIME  CYCLOSTATIONARY  ALGORITHMS 


The  computationally  efficient  algorithms  for  implementation  of  time  and 
frequency-smoothing  techniques  are  discussed  in  [3].  These  are  the  FFT  Accumulation 
Method  (FAM)  and  the  Direct  Frequency-Smoothing  Method  (DFSM)  as  described 
below.  The  temporal  and  spectral  smoothing  equivalence  is  also  addressed  in  [12],  The 
computationally  efficient  algorithms  for  implementation  of  time  and  frequency¬ 
smoothing  techniques  are  extensively  discussed  in  [3]. 


1,  The  Time-Smoothing  FFT  Accumulation  Method  (FAM): 

The  time-smoothing  FFT  Accumulation  Method  was  developed  to  reduce  the 
number  of  computations  required  to  estimate  the  cyclic  spectrum.  [3]  This  technique 
divides  the  bi-frequency  plane  into  smaller  areas  called  the  channel-pair  regions  and 
computes  the  estimates  a  block  at  a  time  using  the  Fast  Fourier  Transform.  Describing  the 
estimated  time-smoothed  periodogram  from  Equations  (3.1.8)  and  (3.1.9),  in  discrete 
terms,  yields 


Sl{n,k) 


N 


N 


N-\ 

z 


n~0 


—  X 
N' 


N' 


f  A 

n,k  +  — 

I  2j 


r 

n,k- 

V 


(3.2.1) 


where 

-jlxkn 

X^,(«,A:)=  ,  (3.2.2) 

«-0 

is  the  Discrete  Fourier  Transform  of  x[n],  w[n]  is  the  data  taper  window  (e.g.  Hamming 
window)  and  the  discrete  equivalents  of  /  and  a  are  k  and  y  respectively.  Figure  23 
represents  a  block  diagram  [13]  used  in  the  implementation  of  this  method  in 
MATLAB®.  [4] 
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Figure  23  FAM  block  diagram  (adapted  from  [3,  13]). 

The  algorithm  consists  of  three  basic  stages:  computation  of  the  complex 
demodulates  (divided  into  data  tapering,  sliding  N'  point  Fourier  transforming  and  base 
band  frequency-downshift  translation  sections),  then  computation  of  the  product 
sequences  and  smoothing  of  the  product  sequences.  Making  a  parallel  between  the 
variables  in  equations  (3.1.8),  (3.1.9)  and  (3.2.1),  we  have: 


NAME 

CONTINOUS  TIME 

DISCRETE  TIME 

SCD 

5;, «./)., 

Short  EFT  Size 

Tw 

N' 

Observation  Time 

At 

N 

Time 

t 

n 

Frequency 

f 

k 

Cycle  Frequency 

a 

Y 

Grenander’s  Uncertainty 

Condition 

M  =  %»1 

II 

Table  1 .  Comparison  of  the  estimated  time-smoothed  Periodogram  expressed  in 
continuous  and  discrete  time. 
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The  parameter  N  represents  the  total  number  of  discrete  instances  within  the  oservatio 
time  and  represents  the  number  of  pints  within  the  discrete  short-time  sliding  FFT.  In 
the  FAM  algorithm,  spectral  components  of  a  sequence,  x[n],  are  computed  using  (3.2.2). 
Two  components  are  multiplied  (3.2.1)  to  provide  a  sample  of  a  cyclic  spectrum  estimate 
representing  a  finite  area  on  the  bi-frequency  plane  called  a  “channel  pair  region,”  as 
shown  in  Figure  24.  There  are  channel  pair  regions  in  the  bi-frequency  plane.  Note 
the  sixteen  small  channel  pair  regions  corresponding  to  a  value  of  A  =  4  in  Figure  24. 

N=4 


Figure  24  Division  of  bi-frequency  plane  in  channel  pair  regions  (adapted  from 
[3,  15]). 

A  sequence  of  samples  for  any  particular  area  may  be  obtained  by  multiplying  the  same 
two  components  of  a  series  of  consecutive  short-time  sliding  FFT's  along  the  entire 
length  of  the  input  sequence.  After  the  channelization  performed  by  an  A' -point  FFT 
sliding  over  the  data  with  an  overlap  of  L  samples,  the  outputs  of  the  FFT’s  were  shifted 
in  frequency  in  order  to  obtain  the  complex  demodulate  sequences  (see  Figure  23).  [3] 
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Instead  of  computing  an  average  of  the  product  of  sequences  between  the  complex 
demodulates,  as  in  (3.1.8),  they  are  Fourier  transformed  with  a  P-point  (second)  FFT.  The 

computational  efficiency  of  the  algorithm  is  improved  by  a  factor  of  L  since  only 
samples  are  processed  for  each  point  estimate.  With  the  sampling  frequency,  the  cycle 
frequency  resolution  of  the  decimated  algorithm  is  defined  as  (compare  to 


),  the  frequency  resolution  is 


(compare  to  Af  = 


)  and  the 


Grenander’s  Uncertainty  Condition  is  (compare  to  At-A/’»1).  Figure  25, 

clearly  reveals  that  the  estimates  toward  the  top  and  the  bottom  (yellow  areas)  of  the 
channel  pair  region  do  not  satisfy  the  Uncertainty  Condition. 


Figure  25  Cycle  frequency  and  frequency  resolutions  and  the  Grenander’s 
Uncertainty  Condition  (adapted  from  [3,  13]). 


In  order  to  minimize  the  variability  of  these  point  estimates,  we  can  retain  only 


K  / 

those  cyclic  spectrum  components  that  are  within  of  y  =  ±  from  the  center  of  the 
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channel  pair  region.  [15]  A  solution  to  resolve  the  entire  area  of  the  ehannel  pair  region 
without  leaving  gaps  is  to  apply  a  data  taper  window  on  the  frequeney  axis  to  obtain 
better  eoverage.  A  Hamming  window  is  applied  in  this  implementation. 


2,  Direct  Frequency-Smoothing  Method: 

Direct  frequency-smoothing  algorithms  first  compute  the  data  frequency 
components  and  then  exeeute  spectral-correlation  operations  direetly  on  the  frequency 
components.  Generally,  the  Direct  Frequency-Smoothing  Method  (DFSM)  is 
computationally  superior  to  indireet  algorithms  that  use  related  quantities,  sueh  as  the 
Wigner-Ville  Distribution.  [12]  But  DFSM  is  normally  less  effieient  than  a  time¬ 
smoothing  approach. 


Time  and  frequeney-smoothing  equivalenee  is  discussed  in  details  in  [10]  and 
derives  from  the  fact  that  randomly  fluctuating  statistical  spectra  (in  both  time  and 

frequeney),  sueh  as  S"  (compare  to  3.1.8,  where  7]^  =  > 


converge  in  the  limit  (At  ^  GO,  A/  O)  to  the  non-random  limit  spectrum  S'®  (/)  if  the 

limit  autocorrelation  (z")  exists.  Employing  discrete  frequency  averaging  and 

following  the  same  analogy  as  in  the  time-smoothed  cyclic  periodogram,  the  basis  for  the 
DFSM  is  the  discrete  time  frequency-smoothed  eyclie  periodogram  represented  by 
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f 

n,k  +  — 

A* 

f  A 

2) 

N 

1  2) 

(3.3.1) 


where 


v-i 


—jlTrkn 

N 


n=0 


(3.3.2) 


is  the  Discrete  Fourier  Transform  of  x\n\,  w\n\  is  a  reetangular  window  of  length  N  that  is 
the  total  number  of  points  of  the  FFT  related  to  the  total  observation  time  At ,  y  is  the 
cycle  frequency  discrete  equivalent,  the  frequency-smoothing  ranges  over  the  interval 

\m\<^y2  Ak~^  frequency  resolution  diserete  equivalent. [12]  All 
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considerations  made  previously  in  Section  A  (Cyclostationary  Theory)  are  still  valid  and, 
for  M  » 1 ,  the  frequency  and  time  smoothing  methods  will  yield  very  similar  results. 
The  block  diagram  on  Figure  26  [13]  illustrates  the  implementation  of  the  “DFSM.m” 
MATLAB®  [4]  routine. 


Figure  26  DFSM  algorithm  block  diagram  (adapted  from  [3,  13]). 

In  order  to  provide  full  coverage  of  the  bi-frequency  plane  at  a  minimal 
computational  expense,  Equation  (3.3.1)  is  computed  along  a  line  of  constant  cycle 

frequency  spacing  the  point  estimates  by  Ak  =  ^  method  is  easier  to 

implement  and  generally  is  used  to  validate  the  time-smoothing  “FAM.m”  code,  but  this 
method  may  become  more  computationally  demanding,  mainly  in  the  last  block  in  which 
the  complex-demodulates  product  sequences  are  summed.  Considerations  on  the  parallel 
processing  of  both  time  and  frequency  algorithms  are  discussed  in  [13].  Combinations  of 
both  methods  are  also  advantageous  for  certain  applications. 

Tthe  FAM  routine  took  considerably  less  time  than  the  DFSM  routine.  The 
computing  time  for  the  DFSM  was  two  or  three  times  longer  than  the  FAM  for  signals 
with  a  large  number  of  samples.  Therefore,  for  very  long  signals  and  large  M  =  AtA/" , 
FAM  should  be  used  to  calculate  the  complete  cyclic  spectrum. 
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3. 


GUI  Implementation: 


A  Graphical  User  Interface  (GUI),  as  shown  in  Figure  27,  was  developed  using 
the  MATLAB®  GUIDE  Version  6  [4]  in  order  to  facilitate  the  execution  of  the 
mathematical  model  processing. 


Hit  <Retum>  after 
all  the  entries  bellow: 

1)  Enter  the  data  file 
directory  (same  format); 

2)  Enter  file  name 
(any  format); 

3)  Enter  sampling 
fi'equeney; 

4)  Choose  the 
frequeney  resolution; 

5)  Choose  the  value  of 
M;  and 

6)  Hit  the  desired 
method  button. 


Figure  27  Cyclostationary  processing  GUI  schematic  tutorial. 


In  this  GUI,  generated  by  the  MATLAB®  [4]  routines  “cyclo.m”,  “FAM.m”  and 
“DFSM.m”,  the  user  defines 

•  The  directory  where  the  data  file  is  to  be  located; 

•  The  name  of  the  file  (without  the  extension); 

•  The  sampling  frequency; 

•  The  desired  frequency  resolution  ( );  and 

•  The  Grenander’s  Uncertainty  Condition  M  =  be  chosen  to  be 

equal  to  2  or  4. 

These  variables  are  then  used  to  calculate 
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Desired  eyele  frequeney  resolution: 


N'  = 


fs. 


-  Number  of  input  ehannels  in  the  ehannelization  matrix: 


(rounded  to  the  nearest  multiple  of  2).  For  example  if 
=  700^  =  875,  then  N'  =  1024; 


Overlap  Parameter:  the  short  time  FFT’s  will  be  performed  and  shifted  by 
L  number  of  samples  for  every  eomputation  (instead  of  at  every  sample): 


Number  of  Columns  in  the  Channelization  Matrix:  P  = 


Total  number  of  points  in  the  input 


proeessed:  N  =  PxL=M  ■ 


fs 


data  to  be 


The  results  of  the  simulations  and  a  tutorial  on  how  they  should  be  interpreted 
will  be  diseussed  later  in  this  ehapter.  We  will  also  be  able  to  eompare  the  results  of  this 
algorithm  and  the  Frequeney-Smoothing  approaeh. 


Sinee  the  number  of  points  to  be  proeessed  N  is  related  to  the  sampling  frequeney, 
the  frequeney  resolution  and  the  value  of  M  are  limited  by  MATLAB®  [4]  funetions 
memory  eonstraints  related  to  the  size  of  the  matrix  being  handled.  The  suggested 
optimum  value  of  N,  in  order  to  avoid  memory  errors,  may  vary  between  2,048  and  4,096 
points  depending  on  the  signal  analyzed  and  may  be  given  by  using  the  suggested  values 
on  the  table  below  when  entering  data  in  the  GUI. 


On  the  Cyelostationary  Proeessing  GUI,  for  an  optimum  N  value  of  2048  points 

II 

Sampling  Frequeney 

7000  Hz 

15000  Hz 

M(%) 

2 

4 

2 

4 

Frequeney  resolution  ( AA:) 

8 

16 

16 

32 

Table  2.  Reeommended  variables  values  for  GUI  users. 
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The  platform  for  all  the  simulations  was  a  Pentium  IV,  2  GHz  CPU  clock  and 
1GB  RAM.  The  most  critical  memory  constraints  appeared  while  using  the  “contour.m” 
MATLAB®  [4]  built-in  function,  but  since  this  is  the  plot  with  the  best  visualization,  we 
decided  to  keep  an  average  maximum  number  of  points  of  2,048  whenever  possible. 

Moreover,  by  increasing  the  value  of  N  from  2048  to  4096  or  more  only  improved 
the  results  a  little.  And  since  the  data  files  were  standardized  for  comparison  with  the 
results  from  the  methods  described  in  [5],  [6]  and  [7],  the  input  signals  were  generated  to 
have  characteristics  that  would  facilitate  measuring  for  all  four  processing  types 
simultaneously. 

All  the  files  in  the  Cyclostationary  MATLAB®  [4]  package  (“cyclo.m,  FAM.m 
and  dfsm.m”)  can  be  edited  very  easily  to  include  new  capabilities.  The  function 
“cyclo.m”  triggers  the  GUI  and  after  entering  the  data,  the  user  may  choose  between  time 
or  frequency-smoothing  method  by  clicking  the  appropriate  button.  Additional  features 
and  options  to  the  GUI  (“Cyclo.fig”)  may  be  added  by  using  the  MATLAB®  GUIDE  [4] 
and  by  editing  “cyclo.m”  accordingly. 


C.  PROCESSING  TUTORIAL 


In  order  to  give  a  general  idea  of  how  the  results  of  both  processing  methods 
should  look,  this  section  provides  some  examples  of  the  signals  analyzed. 

1.  Test  Signals: 


For  this  example,  the  input  complex  signals  characteristics  are  as  following; 


Tab 


Carrier  Frequency  (fc) 

1000  Hz 

Sampling  Frequency  (fs) 

7000  Hz 

e  3.  Test  signal  characteristics. 


This  signal  is  composed  of  a  single  carrier  sampled  at /j=1000Hz;  using  the  real 
part  of  the  signal  we  expect  that  the  plot  on  the  bi-frequency  plane  will  have  the  carrier 
frequency  at  twice  its  value,  as  shown  in  Figure  28  (DFSM  SCD  for  a  Test  Signal  with  N 
=  1024,  df=  16  and  M=  2). 
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•3000  *2000  -1000  0  1000  2000  3000 

Cycle  Frequency  (Hz) 


Figure  28  DFSM  generated  SCD  bi-frequeney  for  a  test  signal  (lOOOFIz  single 
earrier). 


We  see  that  the  frequeney  eontent  is  very  simple,  exeept  for  some  estimation  noise 
around  the  highest  values  of  the  eorrelation  (red  dots).  This  Speetral-Correlation  Density 
was  generated  using  the  Frequeney-Smoothing  Method  (DFSM).  The  test  signal  in  Figure 
29  (FAM  SCD  for  a  Test  Signal  with  A  =  1024,  df=  16  andM=  2)  shows  that  the  Time- 
Smoothing  Method  generates  similar  results  that  are  a  little  bit  noisier  beeause  more 
estimation  is  involved. 


-3000  -2000  -1000  0  1000  2000  3000 

Cycle  Frequency  (Hz) 


Figure  29  FAM  generated  SCD  bi-frequeney  for  a  test  signal  (lOOOHz  single 
earrier). 
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2. 


BPSK: 


For  the  BPSK  (Barker  Code)  example,  the  signals  eharacteristies  are  as  following: 


Carrier  Frequency  (fc) 

lOOOHz. 

Sampling  Frequency  (fs) 

7000Hz. 

Barker  Code  Length  (Np) 

1 1  bits 

SNR 

Only  signal 

Number  of  Carrier  Cycles  per  Barker  Bit  (NPBB) 

1 

Table  4.  BPSK  signal  eharacteristies. 

In  Figure  30,  “k”  is  the  frequency,  “y”  is  the  cycle  frequency,  ‘yd”  is  the  carrier. 


‘X”  is  the  sampling  frequency, is  the  code  rate 


,  “Ay>”  is  the  number 


of  phase  changes,  related  to  the  length  of  the  Barker  sequence  used  in  the  phase 
modulation,  and  “7p”  is  the  chip  period. 
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Figure  30  Pictorial  generic  illustration  of  a  BPSK  signal  SCD  result. 
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By  zooming  in  on  the  plot,  as  shown  in  Figure  31  and  in  the  examples  on  Figure 
32  and  Figure  33,  we  can  measure  the  code  rate  ^fh  \  in  the  cycle  frequency  axis.  It  is 
possible  to  see  a  pattern  of  equally  spaced  points  that  will  eventually  fade  away  as  the 
magnitude  decreases  (from  red  to  blue).  By  changing  the  background  color  of  the  contour 
plots  in  order  to  facilitate  establishing  a  visual  threshold  (around  0.4  on  the  normalized 
values,  or  “light  blue”  on  the  color  bar)  to  eliminate  the  estimation  errors,  we  can 
approximately  determine  the  BW  of  the  signal  being  analyzed  and  this  value  will  then 
divided  by  the  code  rate  giving  the  value  of  ( bandwidth  =  Npxfj^). 


Figure  3 1  Zoomed  in  pictorial  generic  illustration  of  a  BPSK  signal  SCD  result 
with  estimation  of  BW,  number  of  phases  and  code  rate. 
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Figure  32  DFSM  generated  SCD  plot  for  a  BPSK  signal  (lOOOHz  single  earner 
and  1 1  bits  Barker-code  phase  modulation)  with  estimated  BW. 


Cycle  Frequency  (Hz) 


Figure  33  Zoomed-in  DFSM  generated  SCD  plot  for  a  BPSK  signal  (estimated 
code  rate  measurement). 
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The  estimation  of  the  BW  using  this  method  is  somewhat  subjeetive  and  may  lead 
to  still  aeeeptable  errors  on  the  order  of  .  From  the  plots  above  we  may  eonsider  the 
values: 

•  Carrier  frequeney  ifc)  =  1000  Hz  (direetly  from  plot); 

•  Code  rate 

•  BW  (direetly  from  plot)  « 1000//z  ±  90Hz  ; 

•  Number  of  Barker  Bits  (A^)  ~  =  1 1-1 1 «  1 1  bits . 

This  proeedure  of  zooming-in  the  plots  and  measuring  the  distanees  between  the 
highest  magnitudes  (red)  points  will  be  used  for  all  the  signals  analyzed  in  the  next 
ehapter. 


from  zoomed-in  plots, 


-  1/ 


\TpNp) 


=  90  Hz 


J 


3.  FMCW: 


For  the  FMCW  ease,  the  signal  eharaeteristies  are  as  following: 


Carrier  Frequency  (/c) 

1000  Hz. 

Sampling  Frequency  (fs) 

7000  Hz. 

Modulation  Band  Width  ( AF) 

250  Hz 

SNR 

Only  signal 

Modulation  Period  (tm) 

20  ms 

Table  5.  FMCW  signal  eharaeteristies. 

The  same  behavior  showed  by  BPSK  appears  here.  Figure  34  and  Figure  35  show 
that  the  earrier  frequeney  appears  at  twiee  its  value  on  the  eyele-frequeney  axis  and  that 
the  Modulation  BW  and  the  modulation  period  may  also  be  estimated  by  measuring  the 
distanee  from  the  eenter  of  the  SCD  result  to  the  edge,  on  the  A:  =  0  axis  and  by  measuring 
the  distanee  (“t/e/to”)  between  two  eonseeutive  high  eorrelation  points  in  the  eyele 
frequeney  axis. 
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Figure  34  Pictorial  generic  illustration  of  a  FMCW  signal  SCD  result. 


Figure  35  Pictorial  generic  illustration  of  a  zoomed-in  plot  for  a  FCMW  signal 
SCD. 
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Frequency  (Hz) 


This  secondary  variable  called  “delta”  is  related  to  the  modulation  period  as 
shown  in  Figure  35  and  in  the  examples  of  Figure  36  and  Figure  37.  From  Figure  36 
(FAM  SCD  with  N  =  2048,  M  =  4  and  AA:  =  16)  and  Figure  37  (Zoomed  in  plot  of 


Figure  36),  we  get  the  values; 

•  Carrier  frequency  ifc)  =  1000  Flz  (directly  from  plot); 

•  BW  ( AF  ±  delta  )  «  250//z  ±  25Hz ; 

•  ^in~y2  delta  ~  Yl  25~  ~  2t)ms  ,  where  delta=25Hz; 


1700  1800  1900  2000  2100  2200 

Cycle  Frequency  (Hz) 

Figure  36  FAM  generated  SCD  plot  for  a  FMCW  signal  (lOOOHz  carrier  and 
estimated  modulation  BW  of  230  Hz). 
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Figure  37  Zoomed-in  FAM  generated  SCD  plot  for  a  FMCW  signal  ( “delta  ” 
value  of  25Hz). 

4.  P4: 


The  example  signal  eharaeteristies  are  as  following: 


Carrier  Frequency  (/c) 

1000  Hz 

Sampling  Frequency  {fs) 

7000  Hz 

Number  of  Phases  {Np) 

16 

SNR 

Only  Signal 

Number  of  Carrier  Cycles  per  Phase  (cpp) 

1 

Table  6.  P4  signals  eharaeteristies. 

Figure  38  and  Figure  39  show  FAM-generated  SCD  plots.  We  may  see  that  the 
earner  still  appears  at  twiee  its  value  on  the  eyele  frequency  axis,  where  the  pattern  of 
equally  spaced  points  begins  to  become  visible.  Note  the  white  box  drawn  around  the 
data  showing  the  general  shape  of  the  SCD  plot  divided  in  channel  pair  regions.  In  Figure 
39  we  can  see  in  detail  the  estimation  of  the  code  rate  value  {ft).  Figure  40  shows  the 
estimation  of  the  number  of  phases,  as  shown  in  a  pictorial  representation  on  Figure  31. 
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Frequency  (Hz) 
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Figure  38  FAM-generated  SCD  plot  for  a  P4  signal  (1  \25Hz  carrier  and 
estimated  BW  of  1000//z). 


2240  2260  2280  2300  2320 


Cycle  Frequency  (Hz) 

Figure  39  Zoomed-in  FAM-generated  SCD  plot  for  a  P4  signal  (with  estimated 
code  rate  ifh)  of  66Hz). 
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Figure  40  Zoomed-in  FAM-generated  SCD  plot  for  a  P4  signal  (with  estimation 
of  BW  and  number  of  phases). 


Therefore 


•  Carrier  frequeney  (/c)  =  1 125Flz  (direetly  from  plot); 

•  Measured  BW  1^5 IF  =  « 1000//z±66//z  (therefore  c/7p=i); 

•  Number  of  Phases  Used  in  the  P4  Modulation  {Np): 

'^66Hz  ~  ~  phases  ; 


Np  =  BW/  =1000//z/ 

/  Jb 


Code  Period  (7»;  Tp  =  ^PP '  =  1 '  ^^qqq  =  0.01 65 ;  and 


Code  Rate  /,  -  -  (%  gio)  “ 


=  62.5//Z . 
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C.  CHAPTER  SUMMARY 

This  chapter  briefly  discussed  the  theory  related  to  the  cyclostationarity  concept. 
Two  computationally  efficient  cyclostationary  processing  algorithms  (FAM  and  DFSM) 
were  discussed  and  their  implementation  in  MATLAB®  [4]  was  described. 

The  visual  feature  extraction  was  explained  for  some  of  the  modulation  types  that 
will  be  analyzed  in  the  next  chapter.  The  contour  bi-frequency  plots  were  chosen  due  to 
their  greatest  contrast  with  the  figure  background  and  the  convenient  zooming  capability 
of  MATLAB®  [4]  was  explored  to  extraet  parameters  such  as:  carrier  frequency,  code 
rate  and  code  period,  BW,  number  of  phases  or  number  of  Barker  bits  present  on  the 
modulation  and  modulation  period  for  the  FMCW  case.  In  the  next  section,  the  signal 
matrix  shown  in  Table  1  is  analyzed  and  the  theoretical  values  are  compared  with  the 
measurements  from  the  SCD  plots,  according  to  the  tutorial  discussed  in  this  chapter. 
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IV.  DESCRIPTION  OF  LPI  SPECTRAL  PROPERTIES  AND 
CYCLOSTATIONARY  PROCESSING  RESULTS 


In  this  chapter  we  will  analyze  a  set  of  LPI  signals,  with  and  without  the  addition 
of  White  Gaussian  Noise,  and  extract  their  main  characteristics.  All  signals  were 
generated  using  the  “LPI  Generator”  code  from  [5].  A  detailed  description  of  all 
waveforms  may  be  found  in  [5].  Table  8  shows  the  signal  examples  analyzed  in  this 


chapter. 


SIGNAL  NAME 

CARRIER  OR 
HOPPING 
SEQUENCE 
(kHz) 

BW 

(Hz) 

NPBB/ 

CPP 

NUMBER 
OF  BITS  / 
PHASES 
AND 
MOD. 
PERIOD 

MOD.  TYPE 

B  1  7  7  1  s.mat 

1 

1000 

1 

7 

BPSK 

F  1  7  250  20  s.mat 

1 

250 

X 

20 

FMCW 

PI  1  7  4  1  s.mat 

1 

1000 

1 

16 

PI 

P2  1  7  16  1  s.mat 

1 

1000 

1 

16 

P2 

P3  1  7  16  1  s.mat 

1 

1000 

1 

16 

P3 

P4  1  7  16  1  s.mat 

1 

1000 

1 

16 

P4 

FR  1  7  16  1  s.mat 

1 

1000 

1 

16 

FRANK 

C  1  15  10  s.mat 

4,  7,  1,  6,  5,  2,  3 

- 

10 

1 

COSTAS 

FSK  PSK  C  1  15  5  1  s.mat 

4,  7,  1,  6,  5,  2,  3 

1000 

1 

5 

FSK/PSK  COSTAS 

FSK  PSK  T  15  128  5  s.mat 

- 

4200 

5 

128  Hops 

FSK/PSK  TARGET 

Table  7.  List  of  signal  examples  analyzec 

1  in  this  thesis. 

The  signal  analysis  will  consist  of  a  brief  theoretical  description  of  the  signals,  an 
illustration  of  the  speetral  properties  using  plots  of  the  PSD,  and  the  resulting  plots  from 
the  Cyclostationary  processing  algorithms  discussed  previously.  Both  time-smoothing 
and  frequency-smoothing  methods  will  be  used  and  the  method  that  presents  the  best 
visualization  will  be  chosen  for  parameter  extraction.  The  features  extracted  from  the 
signals  will  be  presented  in  tables  comparing  the  measured  results  with  the  given 
parameters. 
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A,  TEST  SIGNALS 

1,  Description 

The  Test  Signals  were  chosen  to  be  single  carrier  signals  with  no  modulation 

•  lOOOHz  Single  Carrier  with  7000Hz  Sampling  Frequency; 

•  lOOOHz  and  2000Hz  Double  Carrier  with  7000Hz  Sampling  Frequency. 

These  signals  give  an  insight  on  how  the  carrier  frequency  information  may  be  extracted 
from  each  plot. 

2,  Spectral  Properties  and  Results  (T_l_7_l_s  and  T_12_7_2_s) 

Figure  41  shows  the  PSD  of  both  Test  Signals. 


Figure  41  PSD  plots  for  both  Test  signals;  a)  lOOOHz  single  sarrier  and  b)  1  and 
2000Hz  double  carrier,  no  modulation. 
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The  single  earner  ease  is  diseussed  on  Chapter  III.  From  Figure  42  and  Figure  43 
we  ean  see  that  the  SCD  plots  for  the  real  part  of  the  double  eanier  signal  show  the  inter¬ 
modulation  produet  at  the  average  of  both  earner  frequeneies  (1500Hz)  and  at  500Hz. 


-6000  -4000  -2000  0  2000  4000  6000 
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Figure  42  DFSM  generated  estimated  SCD  for  a  Test  signal  (lOOOHz  and 
2000Hz  double  earner). 
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Figure  43  Zoomed  in  DFSM  generated  estimated  SCD  for  a  Test  signal  (1000 
and  2000Hz  double  eanier). 
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Table  9  compares  the  measured  and  original  characteristics. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier  Frequency  (/c) 

lOOOHz 

2000Hz 

lOOOHz 

2000Hz 

1500Hz  and  500Hz 

(intermodulation 

products) 

Sampling  Frequency  {f^ 

7000Hz. 

given 

Table  8.  Comparison  between  measured  and  original  characteristics  for  the  Test 
signals. 

The  Frequency-Smoothing  technique  generated  easier  to  visualize  results,  and  for  such 
simple  signals  the  performance  is  equivalent  to  the  Time-Smoothing  algorithm.  A  Time- 


Smoothing  generated  plot  is  shown  in  Chapter  III  for  the  first  Test  Signal  (lOOOHz).  The 
time-smoothing  results  for  the  second  test  signal,  with  double  carriers,  are  not  shown 
since  they  are  very  similar  to  the  DFSM  plots  of  Figure  42  and  Figure  43. 


B,  BPSK 

1,  Description: 

Binary  Phase  Shift  Keying  codes  are  not  considered  LPI  radar  waveforms  in  the 
strict  sense  of  the  word,  but  they  are  here  included  as  a  guide  or  academic  tool  to 
facilitate  the  visualization  of  the  resulting  outputs  from  each  algorithm. 

This  type  of  signal  consists  basically  of  using  Barker  Sequences  of  various 
lengths  to  perform  a  binary  phase  modulation  of  the  carrier  (0  and  n  ).  The  signal  is  then 
sampled  and  the  output  of  the  generator  is  saved  as  two  vectors  (I  and  Q)  representing  the 
in  phase  and  quadrature  components  of  the  signal.  The  modulation  is  performed 
following  the  block  diagram  on  Figure  44  [5]. 
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Figure  44  Block  diagram  for  BPSK  modulation  (from  [5]). 


In  this  work,  we  are  going  to  analyze  a  waveform  modulated  with  a  Barker 
Sequence  of  length  7  bits,  at  various  signal-to-noise  ratios.  All  MATLAB®  [4]  plots  were 
generated  with  the  “cyclo.m”  generated  GUI  and  using  a  frequency  resolution  of  16  Hz 
that  influence  directly  on  how  precise  the  measurements  of  the  code  rate  are  made. 

The  FAM  SCD  plots  were  estimated  using  A=2048  and  M=4,  the  DFSM  SCD  plots  were 
estimated  using  A=1024  andM=2. 

2.  Spectral  Properties  and  Results  (B_l_7_7_l_s) 

The  first  signal  to  be  analyzed  is  a  7-Bit  Barker  Code  phase  modulated  signal, 
with  carrier  frequency  at  lOOOHz  and  with  NPBB  =  1 .  No  noise  is  added  at  this  moment 
and  the  carrier  presents  a  lower  magnitude  in  the  PSD  plot  meaning  that  this  may  be 
called  a  “Suppressed  Carrier  Dual  Side  Band”  (SC-DSB)  Modulation.  The  first  signal 
analyzed  has  its  characteristics  defined  on  Table  10. 


Name 

Carrier 

Sampling 

Number  of 

Number  of  Carrier 

Frequency 

Frequency 

Barker  Bits 

Periods  per  Barker  Bit 

(fc) 

{fs) 

{Np) 

(NPBB) 

B  1  7  7  1  s 

1000  Hz 

7000  Hz 

1 

1 

Table  9.  B_l_7_7_l_s  signal  characteristics. 

The  following  sequence  of  Figures  give  an  overview  of  the  frequency  content  (Figure  45) 
and  the  Estimated  SCD  (Figure  46,  Figure  47,  Figure  48,  Figure  49  and  Figure  50)  of  the 
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Frequency  (Hz) 


signal  analyzed.  The  blue  background  of  the  SCD  estimation  contour  plots  helps  in 
visually  filtering  out  the  noise. 


Figure  45  PSD  for  a  Barker  signal  (lOOOHz  carrier,  7-bits  Barker  sequence  and 
1  NPBB). 
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Figure  46  Estimated  FAM  SCD  contour  plot  for  a  BPSK  real  signal  with 
lOOOHz  carrier,  Barker-7  code  and  1  NPBB. 
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Figure  47  Estimated  FAM  SCD  contour  plot  for  a  BPSK  real  signal  with 

lOOOHz  carrier,  Barker-7  code  and  1  NPBB,  with  estimated  carrier  of 
lOOOHz  and  estimated  BW  of  lOOOHz. 
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Cycle  Frequency  (Hz) 

Figure  48  Zoomed-in  FAM  SCD  contour  plot  for  a  BPSK  signal  with  lOOOHz 
carrier,  Barker-7  code  and  1  NPBB,  with  estimated  of  141  Hz. 
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Figure  49  Estimated  DFSM  SCD  eontour  plot  for  a  BPSK  signal  with  lOOOHz 
earner,  Barker-7  eode  and  1  NPBB,  with  estimated  BW  of  lOOOHz. 
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Figure  50  Zoomed-in  estimated  DFSM  SCD  eontour  plot  for  a  BPSK  signal 
with  lOOOHz  earrier,  Barker-7  eode  and  NPBB  =  1,  with  estimated 
of  142Hz. 
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Frequency  (Hz) 


Figures  51  to  55  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  51  PSD  for  a  Barker  signal  (lOOOHz  carrier,  7-bits  Barker  sequence,  1 
NPBB  and  0  dB  SNR). 
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Figure  52  Estimated  FAM  SCD  contour  plot  for  a  BPSK  signal  with  lOOOHz 
carrier,  7-bits  Barker  code,  1  NPBB  and  0  dB  SNR,  with  estimated 
BW  of  lOOOHz. 
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Figure  53  Zoomed-in  estimated  FAM  SCD  for  BPSK  with  lOOOHz  carrier, 
Barker-7  code,  1  NPBB,  0  dB  SNR,  with  estimated  of  143Hz. 


Figure  54  Estimated  DFSM  SCD  for  BPSK  with  lOOOFIz  carrier,  Barker-7  code, 
1  NPBB,  0  dB  SNR,  with  estimated  BW  of  lOOOHz. 
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Figure  55  Zoomed-in  estimated  FAM  SCD  for  BPSK  with  lOOOHz  earrier, 
Barker-7  code,  1  NPBB,  0  dB  SNR,  with  estimated  of  144Hz. 

The  following  table  shows  a  comparison  between  the  given  and  the  measured 


parameters. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier  Frequency  (fc) 

1000  Hz 

1000  Hz 

Band  Width  (BW) 

1000  Hz 

997.5  Hz 

Carrier  Periods  Per  Barker  Bit 
(NPBB) 

NPBB  =  y^^  =  l 

NPBB  =  5  =  1 .0025 

Code  Rate  (/, 

^100(^  =  142,8 

142.75  Hz  (average) 

Number  of  Barker  Bits  (Np) 

^"100%42.75"  7-005 

Table  10.  Comparison  between  measured  and  original  characteristics  for  This  BPSK 
signal. 


The  BPSK  signals  analyzed  helped  as  an  example  in  determining  how  to  make 
measurements  in  more  complex  signals  such  as  FMCW,  Polyphase  and  frequency 
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hopping  signals.  Table  11  shows  a  summary  of  all  measurements  for  the  BPSK 
modulation.  Table  12  and  Figure  56  show  the  detection  effectiveness  of  the 


cyclostationary  processing  for  all  BPSK  cases,  comparing  with  the  original  values. 


Carrier  (Hz) 

Bandwidth  (Hz) 

Bits 

Code  period  (ms) 

SNR 

1000 

990.5 

7 

7.06 

Only  signal 

1000 

1004.5 

7 

6.968 

0 

1000 

0 

0 

0 

-6 

1000 

995.5 

11 

11.049 

Only  signal 

1000 

1045 

11 

10.526 

0 

1000 

968 

11 

11.363 

-6 

1000 

196 

7 

35.714 

Only  signal 

1000 

192.5 

7 

36.363 

0 

1000 

196 

7 

35.714 

-6 

1000 

187 

11 

58.82 

Only  signal 

1000 

187 

11 

58.82 

0 

1000 

0 

0 

0 

-6 

Table  1 1 .  Summary  of  all  measurements  for  the  BPSK  modulation. 


E 

iPSK  Detection  Ef 

Fectiveness 

Carrier 

Bandwidth 

Code  period 

Bits/code 

Only  signal 

100% 

98% 

103% 

100% 

OdB 

100% 

99% 

102% 

100% 

(-)  6  dB 

75% 

49% 

51% 

50% 

Table  12.  Detection  effectiveness  for  the  BPSK  modulation. 


BPSK  effectiveness  on  detection 


Parameters 


□  Only  signal 

□  OdB 

□  (-)6dB 


Figure  56  Graphic  demonstration  of  detection  effectiveness  for  the  BPSK 
modulation. 
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C.  FMCW 

1.  Description 

The  triangular  modulation  of  a  Frequency  Modulated  Continuous  Wave  (FMCW) 
is  very  popular  as  a  LPI  application  for  a  Continuous  Wave  signal.  The  emitter  uses  a 
continuous  100%  duty-cycle  waveform  and  therefore,  target  range  and  Doppler 
information  can  be  measured  unambiguously  while  retaining  the  LPI  characteristics.  The 
FMCW  signal  is  usually  easier  to  implement  than  phase  code  modulation  signals  as  long 
as  there  is  not  strict  demand  on  linearity  over  the  modulation  BW.  [5]  The  triangular 
modulation  consists  of  two  linear  frequency  modulation  sections  with  positive  and 
negative  slopes.  With  this  configuration,  the  range  and  Doppler  frequency  of  the  detected 
target  can  be  extracted  unambiguously  by  taking  the  sum  and  the  difference  of  the  two 
beat  frequencies.  These  characteristics  are  showed  in  Figure  57. 


Figure  57  Linear  Frequency  Modulated  Triangular  Waveform  and  Doppler 
Shifted  Signal  [5] 
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The  frequency  of  the  transmitted  signal  may  calculated  by 


^  AF  AF 


(4.3.1) 


for  0  <  t  <  tm  and  zero  elsewhere.  Here  fo  stands  for  RF  carrier,  AF  stands  for  transmitted 
modulation  BW,  and  tm  stands  for  modulation  period.  The  rate  of  the  frequency  change  or 
chirp  rate  F  is  [5] 


F  = 


AF 


(4.3.2) 


The  following  section  shows  that  we  can  measure  these  characteristics  of  an  FMCW 
signal  in  the  bi-frequency  plane  generated  by  the  cyclostationary  algorithms 
implemented. 


2,  Spectral  Properties  and  Results  (F_l_7_250_20_s) 


This  FMCW  signal  to  be  analyzed  has  a  modulation  period  of  20  ms,  with  carrier 
frequency  at  lOOOHz  and  modulation  BW  of  250  Hz.  The  signal  has  its  characteristics 
defined  on  Table  14: 


Name 

Carrier 

Sampling 

Modulation 

Modulation  Period 

Frequency 

Frequency 

BW 

ifc) 

ifs) 

(AF) 

F_l_7_250_20_s 

1000  Hz 

7000  Hz 

250  Hz 

20  ms 

Table  13.  F_l_7_250_20_s  signal  characteristics. 


The  PSD  is  shown  in  Figure  58  and  shows  the  frequency  content  of  the  first  FMCW 
signal.  The  following  sequence  of  Figures  gives  an  overview  of  the  frequency  content 
(Figure  58),  time  X  frequency  plot  (Figure  59),  and  the  Estimated  SCD  (Figure  60, 
Figure  61,  Figure  62  and  Figure  63,  all  with  X=  2048,  frequency  resolution  of  16  Flz  and 
M  =  4)  of  the  signal  analyzed. 
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Figure  58  PSD  for  an  FMCW  signal  (lOOOHz  carrier,  250Hz  modulation  BW 
and  20ms  modulation  period,  only  signal). 


X  10'^ 


Figure  59  Time  X  Frequency  plot  for  an  FMCW  signal  (20ms  triangular 
modulation  up-ramp  period). 
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Figure  60  Estimated  FAM  SCD  eontour  plot  for  an  FMCW  signal  with  lOOOHz 
earner  and  estimated  modulation  BW  of  230Hz. 


2110  2115  2120  2125  2130  2135  2140 

Cycle  Frequency  (Hz) 


Figure  61  Zoomed-in  estimated  FAM  SCD  contour  plot  for  an  FMCW  signal 
with  an  estimated  “delta"'  value  of  25Hz. 
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Frequency  (Hz) 


Figure  62  Estimated  DFSM  SCD  contour  plot  for  an  FMCW  signal  with 
lOOOFlz  carrier  and  estimated  modulation  BW  of  235Flz. 


Figure  63  Zoomed-in  estimated  DFSM  SCD  contour  plot  for  an  FMCW  signal 
with  an  estimated  “delta"'  value  of  21Hz. 
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frequency  (Hz) 


Figures  64  to  69  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  64  PSD  for  an  FMCW  signal  (lOOOFlz  carrier,  250  Flz  modulation  BW 
and  20ms  modulation  period,  OdB  SNR). 


1800  1900  2000  2100  2200 

cycle  frequency  (Hz) 


Figure  65  Estimated  FAM  SCD  contour  plot  for  an  FMCW  signal  with  lOOOHz 
carrier,  0  dB  SNR  and  estimated  modulation  BW  of  200  Flz. 
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2110  2115  2120  2125  2130  2135 

cycle  frequency  (Hz) 

Figure  66  Estimated  FAM  SCD  contour  plot  for  an  FMCW  signal  with  an 


estimated  “delta  ”  value  of  26Hz. 


Figure  67  PSD  for  an  FMCW  signal  (lOOOHz  carrier,  250  Hz  modulation  BW 
and  20ms  modulation  period,  -6dB  SNR). 
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Frequency  (Hz) 


Cycle  Frequency  (Hz) 

Figure  68  Estimated  FAM  SCD  contour  plot  for  an  FMCW  signal  with  lOOOHz 
carrier,  -6  dB  SNR  and  estimated  modulation  BW  of  200  Hz. 


2105  2110  2115  2120  2125  2130  2135 


Cycle  Frequency  (Hz) 

Figure  69  Estimated  FAM  SCD  contour  plot  for  an  FMCW  signal  with  an 
estimated  “delta  ”  value  of  27Hz. 
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Table  15  shows  that  the  main  parameters  of  the  FMCW  modulation  may  be  easily 


measured  in  the  bi-frequeney  plane,  using  both  time  and  frequeney-smoothing  methods. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier  Frequency 
(fc) 

1000  Hz 

992.83  Hz  (average) 

Modulation  BW  (AF) 

250  Hz 

218.83  Hz 

^^Delta” 

25  Hz 

24.75  Hz  (average) 

Triangular 

Modulation  Up-Ramp 
Period  (tm) 

C  =  j  7^  =  20ms 

/  2*  delta 

Ki  ~  ^*24  75  “  0.02020s  =  20.20ms 

Table  14.  Comparison  between  measured  and  original  characteristics  for  an  FMCW 


signal. 

Table  16  shows  a  summary  of  all  measurements  for  the  FMCW  modulation.  Table  17  and 
graphics  on  Figure  70  show  the  detection  effectiveness  of  the  cyclostationary  processing 


for  the  FMCW  case,  comparing  with  the  original  values. 


Carrier  (Hz) 

Modulation  BW  (Hz) 

Modulation  Period  (ms) 

SNR 

992.5 

232.5 

21.74 

Only  signal 

993 

208 

19.23 

0 

993 

216 

18.52 

-6 

987.5 

221 

29.42 

Only  signal 

987.5 

234 

27.77 

0 

985 

218.72 

0 

-6 

970 

470 

20.83 

Only  signal 

970 

460 

17.85 

0 

970 

480 

17.86 

-6 

978 

470 

33.33 

Only  signal 

978 

450 

35.71 

0 

970 

480 

0 

-6 

Table  15.  Summary  of  all  measurements  for  the  FMCW  modulation. 


FMCW  Detection  Effectiveness 

Carrier 

Modulation  Bandwidth 

Modulation  Period 

Only  signal 

98% 

92% 

106% 

OdB 

98% 

90% 

99% 

(-)  6  dB 

98% 

91% 

45% 

Table  16.  Detection  effectiveness  for  the  FMCW  modulation. 
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FMCW  effectiveness  on  detection 


Parameters 


Figure  70  Graphic  demonstration  of  detection  effectiveness  for  the  FMCW 
modulation. 


D,  PI 

1,  Description 

The  first  polyphase  code  to  be  analyzed  is  called  “PI”  and  its  phase  modulation  is 
defined  by  the  equation; 


(4.4.1) 

Np 

Where  i  is  the  number  of  the  samples  in  a  given  frequency  and  j  is  the  index  of  the 
frequency,  the  phase  of  the  ith  sample  of  the  jth  frequency  is  given  by  equation  (4.4.1), 
and  i  =  \,2,...,Np  and  j  =  \,2,...,Np.  Therefore,  there’s  a  total  of  Np^  phases  in  this 
modulation.  Figure  71  reveals  the  quadratic  nature  of  the  phase  modulation  in  PI,  given 
by  (4.4.1). 
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Figure  71  PI  code  phase  shift. 


2,  Spectral  Properties  and  Results  (Pl_l_7_16_l_s) 


This  PI  signal  has  the  following  given  characteristics: 


Name 

Carrier 

Sampling 

Number  of 

Number  of  Cycles  per 

Frequency 

Frequency 

Phases 

Phase  {cpp) 

(/c) 

(fs) 

(Np^) 

PI  1  7  16  1  s 

1000  Hz 

7000  Hz 

16 

1 

Table  17.  Pl_l_7_16_l_s  signal  characteristics. 

Based  on  these  properties,  we  expect  a  BW  of  lOOOHz,  and  yields  from  (4.4.1)  that  the 
total  number  of  phases  is  equal  to  Np  and  in  this  case,  the  16  equally  spaced  regions 


inside  the  BW  are  going  to  be  separated  by  = 


=  100^^  =  62.57/z  in  the 


cycle  frequency  axis.  The  PSD  is  shown  in  Figure  72  and  shows  the  frequency  content  of 
This  PI  signal.  The  following  sequence  of  Figures  give  an  overview  of  the  frequency 
content  (Figure  72)  and  the  Estimated  SCD  (Figure  73  and  Figure  74,  with  N  =  1024, 
frequency  resolution  of  16  Hz  and  M=  2)  of  the  signal  analyzed. 
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Frequency  (Hz) 


Figure  72  PSD  for  a  PI  signal  (lOOOHz  carrier,  16  phases  and  1  cpp,  only 
signal). 
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Figure  73  Estimated  DFSM  SCD  contour  plot  for  a  PI  signal  with  900Flz  carrier 
and  estimated  BW  of  lOOOHz. 
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1850  1900  1950 

Cycle  Frequency  (Hz) 

Figure  74  Zoomed-in  estimated  DFSM  SCD  contour  plot  for  a  PI  signal  with  an 
estimated  code  rate  ifb)  of  62Hz. 

Figures  75  to  79  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  75  PSD  for  a  PI  signal  (lOOOHz  carrier,  16  phases  and  1  cpp,  OdB  SNR). 
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Figure  76  Estimated  FAM  SCD  contour  plot  for  a  PI  signal  with  900Hz  carrier 
and  estimated  BW  of  lOOOHz,  with  OdB  SNR. 
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Figure  77  Zoomed-in  estimated  FAM  SCD  contour  plot  for  a  PI  signal  with  an 
estimated  code  rate  {ft)  of  65Hz,  with  OdB  SNR. 
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Cycle  Frequency  (Hz) 

Figure  79  Estimated  FAM  SCD  contour  plot  for  a  PI  signal  with  850Hz  carrier 
and  estimated  BW  of  lOOOHz,  with  -6dB  SNR. 
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The  FAM  measurements  eould  be  performed  until  an  SNR  of-6dB,  exeept  for  the 
ff^  value.  All  results  from  the  DFSM  proeessing  were  similar. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency  (/c) 

1000  Hz 

876.67  Hz  (average) 

Bandwidth 

(BW) 

1000  Hz 

1000  Hz 

Code  Rate  (/*) 

62.5  Hz 

63.5  Hz  (average) 

Number  of 
Phases  (Np^) 

Code  Period 
itp) 

II 

II 

o 

o 

ON 

Table  18.  Comparison  between  measured  and  original  eharacteristies  for  a  PI  signal. 
Both  time  and  frequeney-smoothing  proeessing  teehniques  are  very  well  suited  for 
analyzing  high  BW  and  short  eode  period  signals.  Table  20  shows  a  summary  of  all 
measurements  for  the  PI  modulation. 


Carrier  (Hz) 

Bandwidth  (Hz) 

Number  of  Phases 

Code  period  (ms) 

SNR 

890 

1000 

4 

16.13 

Only  signal 

870 

1000 

4 

15.38 

0 

870 

1000 

0 

0 

-6 

986 

224 

4 

71.43 

Only  signal 

986 

210 

4 

71.43 

0 

986 

200 

0 

0 

-6 

986 

1088 

8 

58.82 

Only  signal 

975 

990 

8 

66.67 

0 

990 

0 

0 

0 

-6 

1000 

200 

8 

312.5 

Only  signal 

1000 

200 

8 

250 

0 

0 

0 

0 

0 

-6 

Table  19.  Summary  of  all  measurements  for  the  PI  modulation. 


Table  21  and  graphies  on  Figure  80  show  the  deteetion  effeetiveness  of  the 
cyelostationary  proeessing  for  the  PI  ease,  comparing  with  the  original  values  for  all  the 
PI  signals  analyzed. 
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PI  Detection  Effectiveness 

Carrier 

Bandwidth 

Code  period 

Phases 

Only  signal 

97% 

105% 

95% 

100% 

OdB 

96% 

101% 

92% 

100% 

(-)  6  dB 

71% 

50% 

0% 

0% 

Table  20.  Detection  effectiveness  for  the  PI  modulation. 


P1  code  effectiveness  on  detection 


Parameters 


□  Only  signal 

□  OdB 

□  (-)6dB 


Figure  80  Graphic  demonstration  of  detection  effectiveness  for  the  PI 
modulation. 


E.  P2 

1,  Description 

The  P2  code  has  almost  the  same  phase  increments  as  the  PI  code,  except  that  the 
starting  phase  is  different.  The  P2  code  is  valid  for  Np  even,  and  each  group  of  the  code 
is  symmetric  about  phase  of  index  0.  These  phases  may  be  calculated  by  the  equation 

(I>,j=\^Wp-\)I  +  (4.5.1) 

Figure  81  reveals  the  symmetric  nature  of  the  P2  code,  for  an  Np^  of  64  phases. 
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Figure  81  P2  symmetric  phase  relationship  between  the  index  in  the  matrix  and 
its  phase  shift. 


2,  Spectral  Properties  and  Results  (P2_l_7_16_l_s) 


This  P2  signal  has  the  following  given  properties: 


Name 

Carrier 

Sampling 

Number  of 

Number  of  Cycles  per 

Frequency 

Frequency 

Phases 

Phase  (cpp) 

(/c) 

ifs) 

(Np^) 

P2  1  7  16  1  s 

1000  Hz 

7000  Hz 

16 

1 

Table  21.  P2_l_7_16_l_s  signal  characteristics. 

From  these  properties  yields  a  BW  of  lOOOHz,  and  from  (4.5.1)  yields  that  the  total 
number  of  phases  is  equal  to  Np  ,  therefore  16  equally  spaced  regions  inside  the  BW  are 


going  to  be  separated  by  .2 

/{Np) 


=  62.5Hz  in  the  cycle  frequency  axis, 


which  is  related  to  the  code  period.  The  PSD  is  shown  in  Figure  82  and  shows  the 
frequency  content  of  This  P2  signal.  The  following  sequence  of  Figures  give  an  overview 
of  the  frequency  content  (Figure  82)  and  the  Estimated  SCD  (Figure  83  and  Figure  84, 
with  A^=  1024,  frequency  resolution  of  16  Hz  and  M=  2)  of  the  signal  analyzed. 
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Frequency  (Hz) 


Cycle  Frequency  (Hz) 

Figure  83  Estimated  FAM  SCD  contour  plot  for  a  P2  signal  with  lOOOHz  carrier 
and  estimated  BW  of  950  Hz. 
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Cycle  Frequency  (Hz) 

Figure  84  Zoomed-in  estimated  FAM  SCD  eontour  plot  for  a  P2  signal  with  an 
estimated  code  rate  (ft)  of  65  Hz. 

Figures  85  to  87  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  85  PSD  for  a  P2  signal  (lOOOHz  carrier,  16  phases  and  1  cpp,  OdB  SNR). 
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Figure  86  Estimated  DSFM  SCD  contour  plot  for  a  P2  signal  with  lOOOHz 
carrier  and  estimated  BW  of  850Hz,  with  OdB  SNR. 
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Cycle  Frequency  (Hz) 

Figure  87  Zoomed-in  estimated  DFSM  SCD  contour  plot  for  a  P2  signal  with  an 
estimated  code  rate  (ft)  of  65Hz,  with  OdB  SNR. 
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The  results  are  summarized  in  Table  23. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency  (fc) 

1000  Hz 

1015  Hz  (average) 

Bandwidth 

(BW) 

1000  Hz 

900  Hz  (average) 

Code  Rate  (/*) 

62.5  Hz 

65  Hz  (average) 

Number  of 
Phases  {Np^) 

=  13.84 

Code  Period  {tp) 

II 

II 

o 

o 

ON 

Table  22.  Comparison  between  measured  and  original  oharaeteristies  for  a  P2  signal. 


Table  24  shows  a  summary  of  all  measurements  for  the  P2  modulation. 


Carrier  (Hz) 

Bandwidth  (Hz) 

Number  of  Phases 

Code  period  (ms) 

SNR 

1000 

975 

3.873 

15.38 

Only  signal 

1030 

870 

3.66 

15.38 

0 

0 

0 

0 

0 

-6 

1000 

167.5 

3.52 

74.07 

Only  signal 

1007 

180 

3.59 

71.43 

0 

1000 

150 

0 

0 

-6 

1000 

962 

8.28 

71.43 

Only  signal 

1100 

1100 

4.5 

18.2 

0 

1050 

1200 

0 

0 

-6 

975 

195 

0 

0 

Only  signal 

1000 

175 

0 

0 

0 

970 

240 

0 

0 

-6 

Table  23.  Summary  of  all  measurements  for  the  P2  modulation. 


Table  25  and  graphies  on  Figure  88  show  the  deteetion  effeetiveness  of  the 


cyelostationary  proeessing  for  the  P2  ease,  eomparing  with  the  original  values. 


P2  Detection  Effectiveness 

Carrier 

Bandwidth 

Code  period 

Phases 

Only  signal 

99% 

94% 

75% 

72% 

OdB 

103% 

94% 

53% 

59% 

(-)  6  dB 

76% 

79% 

0% 

0% 

Table  24.  Deteetion  effeetiveness  for  the  P2  modulation. 
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P2  code  effectiveness  on  detection 
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Figure  88  Graphic  demonstration  of  detection  effectiveness  for  the  P2 
modulation. 


F.  P3 

1,  Description 

The  phase  changes  in  P3  (Figure  89)  occur  according  to  the  following  equation; 

(1-1)4 

I  [(/„+fe)-/,]*  =  ;rt(i-lfl)  (4.6.1) 

0 


where  i  =  1,  2,  . . Np.  Substituting  k=B/T  and  tc=l/B,  the  equation  can  be  written  as 
,(P3)  _  7r{i  —  1)  _  7r(i  — 1) 


BT 


Np 


(4.6.2) 


where  T  stands  for  pulse  length,  f  stands  for  frequency  ( /  =  f„+kt  ),k  is  a  constant,  the 
BW  is  approximately  B=kT,  tc  stands  for  compressed  pulse  length  {t^  -  H B)  and  the 
provided  pulse  compression  ratio  of  pc  =  T  !  BT  is  provided  by  the  waveform. 
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Figure  89  P3  code  phase  shift. 


2,  Spectral  Properties  and  Results  (P3_l_7_16_l_s) 


This  P3  signal  has  the  following  properties: 


Name 

Carrier 

Sampling 

Number  of 

Number  of  Cycles  per 

Frequency 

Frequency 

Phases 

Phase  {cpp) 

(/c) 

(fs) 

{Np) 

P3  1  7  16  1  s 

1000  Hz 

7000  Hz 

16 

1 

Table  25.  P3_l_7_16_l_s  Signal  characteristics. 

The  expected  BW  is  1000  Hz,  and  the  total  number  of  phases  is  equal  to  16,  therefore  16 
equally  spaced  regions  inside  the  BW  are  going  to  be  separated  by 


“  62.5Hz  in  the  cycle  frequency  axis.  The  PSD  is  shown  in  Figure 

90  and  shows  the  frequency  content  of  This  P3  signal.  The  following  sequence  of  Figures 
give  an  overview  of  the  frequency  content  (Figure  90)  and  the  Estimated  SCD  (Figure  91, 
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Figure  92,  Figure  93  and  Figure  94,  with  A^=  1024,  frequency  resolution  of  16  Hz  and  M 
=  2)  of  the  signal  analyzed. 


Figure  90  PSD  for  a  P3  signal  (lOOOHz  carrier,  16  phases  and  1  cpp,  only 


Cycle  Frequency  (Hz) 

Figure  91  Estimated  FAM  SCD  contour  plot  for  a  P3  signal  with  1  lOOHz  carrier 
and  estimated  BW  of  lOOOHz. 
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Figure  93  Estimated  DFSM  SCD  eontour  plot  for  a  P3  signal  with  1 150Hz 
carrier  and  estimated  BW  of  lOOOHz. 
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Figure  94  Zoomed-in  estimated  FAM  SCD  contour  plot  for  a  P3  signal  with  an 
estimated  fb  of  63  Hz. 


Figures  95  to  99  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  95  PSD  for  a  P3  signal  (lOOOHz  carrier,  16  phases,  1  cpp,  and  OdB 
SNR). 
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500  1000  1500  2000  2500  3000  3500 

Cycle  Frequency  (Hz) 

Figure  96  Estimated  FAM  SCD  contour  plot  for  a  P3  signal  with  1050Hz  carrier 
and  estimated  BW  of  lOOOHz,  with  OdB  SNR. 


2120  2140  2160  2180  2200 

Cycle  Frequency  (Hz) 

Figure  97  Zoomed-in  estimated  FAM  SCD  contour  plot  for  a  P3  signal  with  an 
estimated yi,  of  56  Hz  and  OdB  SNR. 
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igure  98 


Cycle  Frequency  (Hz) 

Estimated  DFSM  SCD  contour  plot  for  a  P3  signal  with  1050Hz 
carrier  and  estimated  BW  of  lOOOHz,  with  OdB  SNR. 


1980  1990  2000  2010  2020  2030  2040  2050  2060 


Cycle  Frequency  (Hz) 

Zoomed-in  estimated  DFSM  SCD  contour  plot  for  a  P3  signal  with  an 
estimated  ft  of  68  Hz  and  OdB  SNR. 
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Figure  99 


The  given  and  measured  features  are  summarized  in  Table  27. 


Feature  Extract 

ion 

Characteristic 

Original 

Measured 

Carrier 
Frequency  (fc) 

1000  Hz 

1075  Hz 

Bandwidth 

(BW) 

1000  Hz 

1000  Hz 

Code  Rate  (/*) 

62.5  Hz 

62  Hz  (average) 

Numher  of 
Phases  (Np) 

Np  =  BW/  ^100^^  =  16  ,3 

Code  Period 

II 

II 

o 

o 

ON 

t  =  ^-  =  0.0161295 
p  /62 

Table  26.  Comparison  between  measured  and  original  charaeteristies  for  a  P3  signal. 


Table  28  shows  a  summary  of  all  measurements  for  the  P3  modulation. 


Carrier  (Hz) 

Bandwidth  (Hz) 

Numher  of  Phases 

Code  period  (ms) 

SNR 

1100 

1000 

16.01 

16 

Only  signal 

1050 

1000 

16.03 

16 

0 

0 

0 

0 

0 

-6 

1000 

200 

14.71 

71.43 

Only  signal 

1000 

200 

14.71 

71.43 

0 

1000 

220 

0 

0 

-6 

1100 

1000 

71.43 

71.43 

Only  signal 

1000 

1000 

71.43 

71.43 

0 

0 

0 

0 

0 

-6 

1000 

200 

7.143 

35.71 

Only  signal 

1000 

200 

7.143 

35.71 

0 

0 

0 

0 

0 

-6 

Table  27.  Summary  of  all  measurements  for  the  P3  modulation. 


Table  29  and  graphics  on  Figure  100  show  the  detection  effectiveness  of  the 


cyclostationary  processing  for  the  P3  case,  comparing  with  the  original  values. 


P3  Detection  Effectiveness 

Carrier 

Bandwidth 

Code  period 

Phases 

Only  signal 

105% 

100% 

78% 

79% 

OdB 

101% 

100% 

78% 

79% 

(-)  6  dB 

25% 

28% 

0% 

0% 

Table  28.  Detection  effectiveness  for  the  P3  modulation. 
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P3  code  effectiveness  on  detection 


Parameters 


□  Only  signal 
nOdB 

□  (-)  6dB 

, _ 


Figure  100  Graphic  demonstration  of  detection  effectiveness  for  the  P3 
modulation. 


G.  P4 

1,  Description 

The  P4  code  is  derived  from  the  same  waveform  as  the  P3  code.  However,  in  this 
case,  the  LO  frequency  is  set  equal  to  f^+kT ! 2  in  the  1-Q  detectors.  With  this 
frequency,  the  phases  of  successive  samples  are  defined  by  the  equations;  [5] 


(1-1)4 

(l)f^^=2n  j  [{f^+kt)-f,+kTI2'\dt  =  2n  j  k{t-TI2)dt 


0 


0 


=  nk{i  -  Xf  t]  -  nkT (i  -  \)t^ 


Np 


(4.7.1) 

(4.7.2) 


Figure  101  shows  the  relationship  between  the  index  in  the  matrix  and  its  subsequent 
phase  for  a  P4-coded  signal  with  Np=  64  (phases). 
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Figure  101  Phase  shift  for  a  P4-coded  signal  with  A^=64  phases 


2,  Spectral  Properties  and  Results  (P4_l_7_16_l_s) 

The  P4  signal  has  the  following  characteristics; 


Name 

Carrier 

Sampling 

Number  of 

Number  of  Cycles  per 

Frequency 

Frequency 

Phases 

Phase  (cpp) 

(/c) 

ifs) 

(Np) 

P4  1  7  16  1  s 

1000  Hz 

7000  Hz 

16 

1 

Table  29.  P4_l_7  I6_l_s  signal  characteristics. 


The  expected  BW  is  1000  Hz,  and  the  total  number  of  phases  is  equal  to  16.  There  will  be 
16  equally  spaced  regions  inside  the  BW  separated  by  ~  ^*^*^^6  ~ 

the  cycle  frequency  axis.  The  PSD  is  shown  in  Figure  102  and  shows  the  frequency 
content  of  This  P4  signal.  The  following  sequence  of  Figures  give  an  overview  of  the 
frequency  content  (Figure  102)  and  the  Estimated  SCD  (Figure  103,  Figure  104,  Figure 
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105,  and  Figure  106  with  N  =  1024,  frequency  resolution  of  16  Hz  and  M  =  2)  of  the 
signal  analyzed. 


Figure  102  PSD  for  a  P4  signal  (lOOOHz  carrier,  16  phases  and  1  cpp,  only 
signal). 


Cycle  Frequency  (Hz) 

Figure  103  Estimated  FAM  SCD  contour  plot  for  a  P4  signal  with  1  lOOHz  carrier 
and  estimated  BW  of  1000  Hz. 
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Figure  104  Zoomed- in  estimated  FAM  SCD  contour  plot  for  a  P4  signal  with  an 
estimated  of  66  Hz. 


cycle  frequency  (Hz) 

Figure  105  Estimated  DFSM  SCD  contour  plot  for  a  P4  signal  with  1  lOOHz 
carrier  and  estimated  BW  of  1000  Hz. 
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1980  2000  2020  2040  2(B0  2080 

Cycle  Frequency  (Hz) 

Figure  106  Zoomed- in  estimated  FAM  SCD  eontour  plot  for  a  P4  signal  with  an 
estimated  of  62  Hz. 

Figures  107  to  109  inelude  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  107  PSD  for  a  P4  signal  (lOOOHz  carrier,  16  phases,  1  cpp,  and  OdB 
SNR). 
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Figure  108  Estimated  FAM  SCD  contour  plot  for  a  P4  signal  with  1  lOOHz  carrier 


and  estimated  BW  of  1000  Hz,  with  OdB  SNR. 


2240  2260  2280  2300  2320 


Cycle  Frequency  (Hz) 

Figure  109  Zoomed- in  estimated  FAM  SCD  contour  plot  for  a  P4  signal  with  an 
estimated  of  67  Hz,  with  OdB  SNR. 
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The  given  and  measured  features  are  summarized  in  the  table  below. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency  (fc) 

1000  Hz 

1030  Hz  (average) 

Bandwidth 

(BW) 

1000  Hz 

1000  Hz 

Code  Rate  (/*) 

62.5  Hz 

65.5  Hz  (average) 

Number  of 
Phases  (Np) 

Code  Period 
(4) 

II 

II 

o 

o 

ON 

C  =  .  =  0.015275 

p  /65.5 

Table  30.  Comparison  between  measured  and  original  charaeteristies  for  a  P4  signal. 


Table  32  shows  a  summary  of  all  measurements  for  the  P4  modulation. 


Carrier  (Hz) 

Bandwidth  (Hz) 

Number  of  Phases 

Code  period  (ms) 

SNR 

1100 

1000 

15.625 

16 

Only  signal 

1100 

1000 

14.93 

16 

0 

1000 

700 

0 

0 

-6 

1000 

200 

14.29 

71.43 

Only  signal 

1000 

200 

14.29 

71.43 

0 

960 

120 

0 

0 

-6 

1000 

1000 

71.43 

71.43 

Only  signal 

1000 

1000 

71.43 

71.43 

0 

1000 

1100 

0 

0 

-6 

1000 

160 

22.86 

142.86 

Only  signal 

1000 

170 

12.59 

74.1 

0 

0 

0 

0 

0 

-6 

Table  3 1 .  Summary  of  all  measurements  for  the  P4  modulation. 


Table  33  and  graphics  on  Figure  110  show  the  detection  effectiveness  of  the 


cyclostationary  processing  for  the  P4  case,  comparing  with  the  original  values. 


P4  Detection  Effectiveness 

Carrier 

Bandwidth 

Code  period 

Phases 

Only  signal 

103% 

95% 

86% 

84% 

OdB 

103% 

96% 

81% 

78% 

(-)  6  dB 

74% 

60% 

0% 

0% 

Table  32.  Detection  effectiveness  for  the  P4  modulation. 
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P4  code  effectiveness  on  detection 


Parameters 


□  Only  signal 

□  0  dB 

□  (-)6dB 


Figure  110  Graphic  demonstration  of  detection  effectiveness  for  the  P4 
modulation. 

H,  FRANK 

1,  Description 

Frank  codes  belong  to  the  family  of  polyphase  codes,  Chirp  codes  and  Barker 
codes,  and  have  been  successfully  used  in  implementing  Low  probability  of  Intercept 
(LPI)  radar  signals.  [5] 

■y 

A  Frank-coded  waveform  has  a  length  of  Ap  and  basically  consists  of  a  constant 
amplitude  signal  whose  carrier  frequency  is  modulated  by  the  phases  of  the  Frank  code 
given  by  the  following  equation; 

(4.8.1) 

^  Ap 

where  i  =  number  of  samples  (i  =  1,2...  Np)  and  j  =  number  of  frequencies  (j  =  1,2,... 
Ap).  Figure  111  clearly  displays  how  the  phases  change  in  the  Frank  modulation. 
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2,  Spectral  Properties  and  Results  (FR_l_7_16_l_s) 


This  Frank  signal  has  the  following  characteristics; 


Name 

Carrier 

Sampling 

Number  of 

Number  of  Cycles  per 

Frequency 

Frequency 

Phases 

Phase  (cpp) 

(/c) 

ifs) 

(Np^) 

FR  1  7  16  1  s 

1000  Hz 

7000  Hz 

16 

1 

Table  33.  FR_l_7_16_l_s  signal  characteristics. 

The  expected  BW  is  1000  Hz,  and  the  total  number  of  phases  is  equal  to  16.  There  will  be 
sixteen  equally  spaced  regions  inside  the  BW  range  separated  by  the  distance 


/.= 


=  62.5Hz  in  the  cycle  frequency  axis.  The  PSD  is  shown  in 


Figure  112  and  shows  the  frequency  content  of  This  Frank  signal.  The  following 
sequence  of  Figures  give  an  overview  of  the  frequency  content  (Figure  112)  and  the 
Estimated  SCD  (Figure  113  and  Figure  114,  with  N  =  1024,  frequency  resolution  of  16 
Hz  and  M  =  2)  of  the  signal  analyzed. 
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Frequency  (Hz) 


Figure  112  PSD  for  a  Frank  signal  (lOOOHz  carrier,  16  phases  and  1  cpp,  only 
signal). 


Cycle  Frequency  (Hz) 

Figure  113  Estimated  DFSM  SCD  contour  plot  for  a  Frank  signal  with  1 150Hz 
carrier  and  estimated  BW  of  1000  Hz. 
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1980  1990  2000  2010  2020  2030  2040  2050  2060  2070 
Cycle  Frequency  (Hz) 

Figure  114  Zoomed- in  estimated  DFSM  SCD  contour  plot  for  a  Frank  signal  with 
an  estimated yi,  of  61  Hz. 

Figures  115  to  120  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 


Figure  115  PSD  for  a  Frank  signal  (lOOOHz  carrier,  16  phases,  1  cpp,  and  OdB 
SNR). 
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Cycle  Frequency  (Hz) 

Figure  116  Estimated  FAM  SCD  contour  plot  for  a  Frank  signal  with  lOOOHz 
carrier  and  estimated  BW  of  1000  Hz,  with  OdB  SNR. 


2000  2020  2040  2060  2060 


Cycle  Frequency  (Hz) 

Figure  117  Zoomed- in  estimated  FAM  SCD  contour  plot  for  a  Frank  signal  with 
an  estimated  ft  of  63  Hz  and  OdB  SNR. 


100 


Frequency  (Hz) 


15 


10 


-5 


Cycle  Freauencv  (Hzl 

Figure  119  Estimated  FAM  SCD  contour  plot  for  a  Frank  signal  with  1  lOOHz 
carrier  and  estimated  BW  of  1000  Hz,  with  -6dB  SNR. 
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PSD  for  a  Frank  signal  (lOOOHz  carrier,  16  phases,  1  cpp,  and  -6dB 
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Figure  120  Zoomed- in  estimated  FAM  SCD  contour  plot  for  a  Frank  signal  with 
an  estimated  fb  of  65  Hz,  and  -6dB  SNR. 


The  given  and  measured  features  are  summarized  in  Table  35. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency  (/c) 

1000  Hz 

1083.3  Hz  (average) 

Bandwidth  {BW) 

1000  Hz 

1000  Hz 

Code  Rate  (/*) 

62.5  Hz 

63  Hz  (average) 

Number  of 
Phases  (Np^) 

Np^=BW^^  ^100^^^15,87 

Code  Period  {tp) 

r  =  ^,  =  0.015875 
p  /63 

Table  34.  Comparison  between  measured  and  origina 

characteristics  for  a  Frank 

signal. 
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Table  36  shows  a  summary  of  all  measurements  for  the  Frank  modulation. 


Carrier  (Hz) 

Bandwidth  (Hz) 

Number  of  Phases 

Code  period  (ms) 

SNR 

1150 

1000 

4.05 

16.4 

Only  signal 

1000 

1000 

3.98 

15.87 

0 

1100 

1000 

3.92 

15.38 

-6 

1000 

200 

3.65 

66.67 

Only  signal 

1000 

200 

3.65 

66.67 

0 

1000 

210 

0 

0 

-6 

1100 

1000 

7.45 

55.55 

Only  signal 

1100 

1000 

8.45 

71.43 

0 

1100 

1000 

0 

0 

-6 

1000 

200 

4.71 

111.11 

Only  signal 

1000 

200 

3.92 

76.92 

0 

0 

0 

0 

0 

-6 

Table  35.  Summary  of  all  measurements  for  the  Frank  modulation. 


The  following  table  and  graphics  on  Figure  121  show  the  detection  effectiveness  of  the 


cyclostationary  processing  for  the  Frank  case,  comparing  with  the  original  values. 


Frank  Detection  Effec 

tiveness 

Carrier 

Bandwidth 

Code  period 

Phases 

Only  signal 

106% 

100% 

77% 

86% 

OdB 

103% 

100% 

80% 

86% 

(-)  6  dB 

80% 

76% 

24% 

25% 

Table  36.  Detection  effectiveness  for  the  Frank  modu 

ation. 

Frank  code  effectiveness  on  detection 


Parameters 


□  Only  signal 
nO  dB 

□  (-)6dB 


Figure  121  Graphic  demonstration  of  detection  effectiveness  for  the  Frank 
modulation. 
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L  COSTAS  CODES 

1.  Description 

We  briefly  diseussed  the  Costas  modulation  on  Chapter  II.  If  we  don’t  phase 
modulate  the  signals  diseussed,  we  get  a  pure  Costas  FH  signal.  The  FH  signal  eonsists  of 
one  or  more  frequeneies  being  ehosen  from  a  set  ...fm)  of  available  frequeneies, 

for  transmission  at  eaeh  of  a  set  {t\,  t2,  ...  tn}  of  eonseeutives  time  intervals.  For 
modeling  purposes,  it  is  reasonable  to  eonsider  the  situation  in  whieh  m=n,  and  a 
different  one  of  n  equally  spaeed  frequeneies  ...fn)  is  transmitted  during  eaeh  of 

the  n  equal  duration  time  intervals  {t\,  tl,  ...  tn}.  Sueh  a  signal  is  represented  by  a  n  x  n 
permutation  matrix  A,  where  the  n  rows  eorrespond  to  the  n  frequeneies,  the  n  columns 
correspond  to  the  n  intervals,  and  the  entry  Uy  =  1  means  transmission  and  0  otherwise. 
[5]  Any  frequency  sequence  chosen  from  a  Costas  matrix  will  generate  a  Costas 
Frequency  Hopping  sequence  with  the  spectral  characteristics  discussed  in  [5].  In  order  to 
facilitate  the  analysis  of  these  signals,  the  complex  waveform  was  analyzed  and  this  will 
result  in  a  less  cluttered  plot,  for  this  type  of  signal. 


2.  Spectral  Properties  and  Results  (C_l_15_10_s) 

The  signal  characteristics  are  defined  on  Table  38. 


Name 

Frequency 

Sequence  (kHz) 

Sampling 

Frequency 

ifs) 

Number  of  Cycles  per 

Frequency  (cpj) 

C_l_15_10_s 

4716523 

15000  Hz 

10 

Table  37.  C_l_15_10_s  signal  characteristics. 


The  expected  BW  is  6000Hz  (from  lOOOHz  to  7000Hz),  and  the  total  number  of 
frequencies  is  equal  to  7.  The  following  Figures  give  an  overview  of  the  frequency 
content  of  the  signal  (Figure  122)  and  the  Estimated  SCD  (Figure  123,  Figure  124,  Figure 
125,  Figure  126  and  Figure  127,  with  N=  1024,  frequency  resolution  of  16  Hz  and  M  = 
2)  of  the  signal  analyzed.  Note  on  Figure  123  that  the  frequency  components  are  all  along 
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the  /  =  0  axis  (vertical  or  frequency  axis)  while  the  intermodulation  products  are  located 
symmetrically  over  the  bi-frequency  plane,  outside  the  vertical  axis. 


Figure  122  PSD  for  a  Costas  signal  (1,  2,  3,  4,  5,  6  and  7  kHz  carriers,  10  cpf, 
only  signal). 


Cycle  Frequency  (Hz) 

Figure  123  Estimated  FAM  SCD  contour  plot  for  a  complex  Costas  signal  (1,  2, 
3,  4,  5,  6  and  7000Hz  carriers  over  y  =  0  axis,  cpj=\0,  only  signal), 
with  intermodulation  products. 
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Figures  124  to  127  include  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 


measured  characteristics  of  each  signal. 


Figure  124  PSD  for  a  Costas  signal  (1,  2,  3,  4,  5,  6  and  7kHz  carriers,  10  cpf  and 
OdB  SNR). 
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Figure  125  Estimated  FAM  SCD  contour  plot  for  a  complex  Costas  signal  (1,  2, 
3,  4,  5,  6  and  7kHz  carriers  over  7  =  0  axis,  10  cpf,  OdB  SNR),  with 
intermodulation  products. 
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Frequency  (Hz) 


16 


Figure  126  PSD  for  a  Costas  signal  (1,  2,  3,  4,  5,  6  and  7kHz  carriers,  10  cpf  and 


SNR  of-6dB). 


-6000  -4000  -2000  0  2000  4000  BCOD 


Cycle  Frequency  (Hz) 

Figure  127  Estimated  FAM  SCD  contour  plot  for  a  complex  Costas  signal  (1,  2, 
3,  4,  5,  6  and  7kHz  carriers  over  7  =  0  axis,  10  cpf,  -6dB  SNR),  with 
intermodulation  products. 
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The  order  of  the  frequeneies  ean’t  be  measured  in  this  kind  of  plot  sinee  we  don’t 


have  a  '"time  x  frequency"  type  of  output,  but  the  frequency  content  is  easily  determined. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Frequency 

Components 

1,  2,  3,  4,  5,  6  and  7kHz 

1,  2,  3,  4,  5,  6  and  7kHz 

BW  (BW) 

6000Hz 

6000Hz 

Table  38.  Comparison  between  measured  and  original  characteristics  for  a  Costas 


signal. 


The  following  table  shows  a  summary  of  all  measurements  for  the  Costas  modulation. 


Sequence 

Time  in  each  freq 

Code  Period 

SNR 

4716523 

0 

0 

Only  signal 

4716523 

0 

0 

0 

4716523 

0 

0 

-6 

4716523 

0 

0 

Only  signal 

4716523 

0 

0 

0 

4716523 

0 

0 

-6 

4167532 

0 

0 

Only  signal 

4167532 

0 

0 

0 

4167532 

0 

0 

-6 

4167532 

0 

0 

Only  signal 

4167532 

0 

0 

0 

4167532 

0 

0 

-6 

Table  39.  Summary  of  all 

.  measurements  for  the  Costas  modulation. 

The  following  table  and  graphics  on  Figure  128  show  the  detection  effectiveness  of  the 


cyclostationary  processing  for  the  Costas  case,  comparing  with  the  original  values. 


Costas  Detection  Effectiveness 

Sequence 

Time  in  each  Frequency 

Code  period  (ms) 

Only  signal 

100% 

0% 

0% 

OdB 

100% 

0% 

0% 

(-)  6  dB 

100% 

0% 

0% 

Table  40. 


Detection  ef 


fectiveness  for  the  Costas  modulation. 
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Costas  code  effectiveness  on  detection 


Parameters 


□  Only  signal 
nOdB 

□  (-)6dB 


Figure  128  Graphic  demonstration  of  detection  effectiveness  for  the  Costas 
modulation. 

J.  FSK/  PSK  COSTAS 

1,  Description 

A  detailed  description  of  this  type  of  signal  is  given  in  Chapter  11. 


2,  Spectral  Properties  and  Results  (FSK_PSK_C_l_15_5_l_s) 

This  FSK/PSK  Costas  signal  has  the  following  characteristics: 


Name 

Frequency 

Sequence 

(kHz) 

Sampling 

Frequency 

ifs) 

Number  of 

Barker 

Bits  (Np) 

Number  of 

Cycles  per 

Phase 

(NPBB) 

FSK  PSK  C  1  15  5  1  s 

4716523 

15000  Hz 

5 

1 

Table  41.  FSK  PSK  C  1 

15  5  1s  signa 

characteristics. 

The  expected  BW  is  7000  Hz  (from  500  Hz  to  7.5kHz),  each  frequency  will  be 
modulated  with  as  BPSK  Barker-5  and  since  the  smallest  frequency  is  lOOOHz,  the  BW 
for  each  frequency  hop  will  be  lOOOHz.  Inside  each  hop,  we  expect  to  see  a 
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/.= 


BW/ 


1000/ 


'Np  ~  /5~  following  Figures  give  an  overview  of  the  frequency 

content  of  the  signal  (Figure  129),  the  Barker  phase  sequence  (Figure  130)  and  the 
Estimated  SCD  (Figure  131,  Figure  132  and  Figure  133,  with  N  =  1024,  frequency 
resolution  of  16  Hz  andM=  2)  of  the  signal  analyzed. 


Figure  129  PSD  for  a  FSK/PSK  Costas  signal  (1,  2,  3,  4,  5,  6  and  7kHz  carriers, 
Barker-5  and  1  NPBB,  only  signal). 


Figure  130  Barker-5  phase  sequence  used  inside  each  hop. 
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131  Estimated  FAM  SCD  contour  plot  for  a  complex  FSK/PSK  Costas 
signal  (1,  2,  4,  5,  6  and  7kHz  measured  carriers). 
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Figure  132  Zoomed- in  estimated  FAM  SCD  contour  plot  for  a  complex 
FSK/PSK  Costas  signal  (4,  5  and  6kHz  measured  carriers  and 
estimated  BW  of  1000  Hz  for  each  frequency  hop). 


Ill 


Figure  133  Estimated  fb  value  of  200  Flz  for  the  embedded  Barker-5  BPSK 
modulation. 

Figures  134  to  137  inelude  the  analysis  of  the  same  signal  with  added  White  Gaussian 
Noise.  At  the  end  of  this  analysis,  a  table  is  included  comparing  the  original  and  the 
measured  characteristics  of  each  signal. 
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Figure  134  PSD  for  a  FSK/PSK  Costas  signal  (1,  2,  3,  4,  5,  6  and  7kHz  carriers, 
Barker-5  and  1  NPBB,  OdB  SNR). 
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Frequency  (Hz)  31  Frequency  (Hz) 
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135  Estimated  FAM  SCD  contour  plot  for  a  complex  FSK/PSK  Costas 
signal  (1,  2,  4,  5,  6  and  7kHz  measured  carriers,  OdB  SNR). 
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Figure  136  Zoomed-in  estimated  FAM  SCD  contour  plot  for  a  complex 
FSK/PSK  Costas  signal  (5,  6  and  7kHz  measured  carriers  and 
estimated  BW  of  1000  Hz  for  each  frequency  hop,  OdB  SNR). 


113 


6040 

6020 

6000 

5980 

I  5960 

& 

I  5940 

O’ 

o 

^  5920 
5900 


5880 

5860 

-50  0  50  100  150  200  250 

Cycle  Frequency  (Hz) 

Figure  137  Estimated  ft  value  of  200  Hz  for  the  embedded  Barker-5  BPSK 
modulation. 


The  given  and  measured  features  are  summarized  in  the  table  below. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Carrier 
Frequency  (fc) 

1,  2,  3,  4,  5,  6  and  7kHz 

1,  4,  5,  6  and  7kHz 

BW  for  Each 
Frequency  Hop 
(BW) 

1000  Hz 

1000  Hz 

Code  Rate  (/*) 

200  Hz 

200  Hz 

Number  of  Bits  in 
Barker  Code  (Np) 

m^-BW/  -1000/  -5 

/fb~  /200  ^ 

Code  Period 
Inside  Each  Hop 

=0,005. 

o 

o 

o 

II 

o 

o 

II 

o. 

4.. 

Table  42.  Comparison  between  measured  and  original  characteristics  for  a  FSK/PSK 
Costas  signal. 
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The  large  BW  within  eaeh  frequeney  hop  may  increase  the  error  in  the  estimation  of  the 
carrier  frequencies.  For  the  example  above,  it  was  possible  to  measure  practically  all 
characteristics  of  the  signals,  with  SNR  until  0  dB.  These  signals  present  a  modulation 
complexity  level  much  higher  than  any  modem  receiver  would  be  able  to  analyze.  The 
cyclostationary  processing  algorithms  implemented  reveal  their  real  processing  power 
when  facing  these  types  of  signals:  frequency  hopping  with  large  intra-hop  BWs  and 
small  code  periods.  Both  cyclostationary  algorithms  are  very  powerful  tools  to  perform 
analysis  of  LPI  signals  like  the  FSK/PSK  Costas.  Still  the  biggest  constraints  reside  in 
processing  power  of  the  analyzer  (computer  speed  and  memory).  Table  44  shows  a 


summary  of  all  measurements  for  the  Costas  modulation. 


Sequence 

Bandwidth 

Np 

Code  Period 

SNR 

4  7  1  65  2 

7000 

5 

5 

Only  signal 

4  7  1  65  2 

7000 

5 

5 

0 

47  1 

6200 

5 

25 

Only  signal 

47  1 

6200 

5 

25 

0 

47  1 

7000 

11 

11 

Only  signal 

47  1 

7000 

11 

11 

0 

4 

200 

11 

55 

Only  signal 

4 

200 

11 

55 

0 

Table  43.  Summary  of  all  measurements  for  the  FSK/PSK  Costas  modulation. 

Table  45  and  graphics  on  Figure  138  show  the  detection  effectiveness  of  the 
cyclostationary  processing  for  the  FSK/PSK  Costas  case,  comparing  with  the  original 
values. 


FSK/PSK  Costas  Detection  Effectiveness 

Sequence 

Bandwidth  (Hz) 

Code  period  (ms) 

Only  signal 

47% 

76% 

100% 

OdB 

47% 

76% 

100% 

Table  44.  Detection  effectiveness  for  the  FSK/PSK  Costas  modulation. 
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Figure  138  Graphic  demonstration  of  detection  effectiveness  for  the  FSK/PSK 
Costas  modulation. 

K.  FSK/  PSK  TARGET 

1,  Description 

Vide  Chapter  11  for  a  detailed  description  of  these  types  of  signals. 


2,  Spectral  Properties  and  Results  (FSK_PSK_T_15_128_5_s) 

This  FSK/PSK  Target  signal  has  the  following  characteristics: 


Name 

Sampling 

Frequency 

ifs) 

Number  of 

Frequency  Hops 

m 

Number  of  Cycles 

per  Phase  (cpp) 

FSK_PSK_T_15_128_5_s 

15000  Hz 

128 

5 

Table  45.  FSK_PSK_T_15_128_5_s  Signal  characteristics. 

The  expected  BW  is  around  4200Flz  (from  2400Hz  to  6600Hz),  each  frequency  will  be 
modulated  with  a  random  phase  sequence  of  length  equal  to  the  number  of  cycles  per 
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phase.  The  following  Figures  reveal  the  frequency  content  of  the  signal  (Figure  139),  the 
phase  plots  (Figure  140),  the  histogram  of  the  frequency  hops  (Figure  141)  and  the 
Estimated  SCD  (Figure  142,  with  A^=  1024,  frequency  resolution  of  16  Hz  and  M=  2)  of 
the  signal  analyzed. 


Figure  139  PSD  for  a  FSK/PSK  Target  signal  (4200Hz  BW,  random  phase  with 
length  5  and  5  cpp,  only  signal). 


Figure  140  Random  phase  sequence  of  length  5  used  inside  each  hop. 
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Figure  141  Frequency  hops  histogram  after  random  tiring  order 
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Figure  142  Estimated  FAM  SCD  contour  plot  for  a  complex  FSK/PSK  Target 
signal  with  an  estimated  BW  of  4800Flz. 
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The  given  and  measured  features  are  summarized  in  the  table  below. 


Feature  Extraction 

Characteristic 

Original 

Measured 

Total  BW  for  Entire  Frequency  Hopping  Sequence 

4200Hz 

4800Hz 

Table  46.  Comparison  between  measured  and  original  oharaeteristies  for  a  FSK/PSK 
Target  signal. 


The  highest  peaks  of  the  speetrum  of  this  signal  (as  seen  in  the  PSD  plot)  may  also  be 
deteeted  in  the  bi-frequeney  plane.  The  firing  order  of  the  frequeney  hopping  sequenee 
eannot  be  measured  beeause  we  don’t  have  the  time  information  in  the  output  of  the 
eyelostationary  proeessing. 

This  signal  presented  the  highest  degree  of  eomplexity  of  the  whole  set  of  signals 
analyzed.  The  eause  of  this  high  eomplexity  may  be  linked  to  the  presenee  of  two 
pseudo-random  faetors: 

•  The  frequeney  firing  order;  and 

•  The  phase  modulation  for  each  frequency  hop. 

Since  the  emitter  knows  the  pseudo-random  sequence  used,  it  can  coherently  demodulate 
the  signal  with  little  effort  but  for  a  non-cooperative  receiver  this  task  becomes  less 
feasible.  Table  48  shows  a  summary  of  all  measurements  for  the  FSK/PSK  Target 
modulation. 


Sequence 

Bandwidth 

Code  Period 

SNR 

0 

4800 

0 

Only  signal 

0 

3000 

0 

Only  signal 

0 

4200 

0 

Only  signal 

0 

3000 

0 

Only  signal 

Table  47.  Summary  of  all 

.  measurements  for  the  FSK/PSK  Target  modulation. 

Table  49  and  graphics  on  Figure  143  show  the  detection  effectiveness  of  the 
eyelostationary  processing  for  the  FSK/PSK  Target  case,  comparing  with  the  original 
values.  The  analysis  was  performed  for  the  only  signal  case,  since  the  results  from  the 
added  White  Gaussian  Noise  cases  were  all  inconclusive. 
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FSK/PSK  Target  E 

►election  Effectiveness 

Sequence 

Bandwidth 

(Hz) 

Code  period 
(ms) 

Only  signal 

0% 

95% 

0% 

Table  48.  Detection  effectiveness  for  the  FSK/PSK  Target  modulation. 


FSK/PSK  Target  effectiveness  on  detection 


□  Only  signal 

□  OdB 


Parameters 


Figure  143  Graphic  demonstration  of  detection  effectiveness  for  the  FSK/PSK 
Target  modulation. 

L.  COMPARISON  BETWEEN  POLYPHASE  CODES 

The  following  figures  show  a  comparison  between  Frank  (Figure  144),  PI  (Figure 
145),  P2  (Figure  146),  P3  (Figure  147)  and  P4  (Figure  148),  with  the  same  number  of 
phases  {l^p  or  l^p  =  16),  same  carrier  frequency  (IkHz)  and  same  number  of  carrier 
cycles  per  phase  {cpp=\').  It  is  possible  to  see  that  there’s  some  difference  between  the 
SCD  estimated  results  for  the  various  modulations.  The  general  shape  of  Frank,  P3  and 
P4  are  very  similar  but  PI  and  P2  stand  out  very  easily  from  the  group.  The  absence  of 
time  in  the  plots  prevents  us  to  see  the  phase  changes  along  the  code.  The  number  of 
phases  may  be  calculated,  as  seen  on  the  previous  sections,  by  measuring  the  carrier 
frequency,  the  bandwidth  and  the  code  period.  Both  cyclostationary  algorithms  were  very 
effective  on  determining  all  parameters  for  the  LPI  Radar  signals  with  larger  bandwidth 
(cp/>=l)  and  shorter  code  periods. 
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Figure  144  Estimated  DFSM  SCD  contour  plot  for  a  Frank  signal  with  1 150Hz 
carrier,  BW  of  1000  Hz,  Np^=\6  and  cpp=\. 
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Figure  145  Estimated  DESM  SCD  contour  plot  for  a  PI  signal  with  900Hz 
carrier,  BW  of  1000  Hz,  Np^=\6  and  cpp=\. 
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Estimated  DFSM  SCD  contour  plot  for  a  P2  signal  with  1050Hz 
carrier,  BW  of  950  Hz,  Np^=\6  and  cpp=\. 
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Figure  147  Estimated  DFSM  SCD  contour  plot  for  a  P3  signal  with  1 150Hz 
carrier,  BW  of  1000  Hz,  Np=\6  and  cpp=\. 
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cycle  frequency  (Hz) 

Figure  148  Estimated  DFSM  SCD  contour  plot  for  a  P4  signal  with  1 150Hz 
carrier,  BW  of  1000  Flz,  Np=\6  and  cpp=\. 


The  DFSM  plots  were  chosen  just  to  establish  a  standard  observation  but  the  FAM 
generated  plots,  as  seen  on  the  previous  sections,  show  the  same  results  and  take  less 
computing  time  to  process. 


M.  CHAPTER  SUMMARY 

In  this  last  Chapter  we  analyzed  and  extracted  parameters  from  a  set  of  FPl  Radar 
Signals.  The  extracted  parameters  were  compared  to  the  original  ones  in  tables  at  the  end 
of  every  processing. 

The  analyzed  modulations  were 

•  BPSK  (Binary  Phase  Shift  Keying); 

•  FMCW  (Frequency  Modulated  Continuous  Wave); 

•  Polyphase  Codes  (P4,  P3,  P2  PI  and  Frank  Codes); 
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•  Costas  Codes  (Frequency  Hopping); 

•  FSK/PSK  (Combined  Frequency  Shift  Keying  and  Phase  Shift  Keying)  with  a 
Costas  frequency  distribution;  and 

•  FSK/PSK  (Combined  Frequency  Shift  Keying  and  Phase  Shift  Keying)  with  a 
target  matched  frequency  distribution 

It  was  possible  to  see  that  both  time  and  frequency-smoothing  methods  of 
cyclostationary  processing  may  be  applied  to  analyze  these  types  of  signals  yielding 
fairly  appropriate  results  for  the  majority  of  the  signals.  A  comparison  between  the 
polyphase  codes  is  also  presented. 
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V.  CONCLUSIONS  AND  RECOMMENDATIONS 


This  thesis  effort  was  developed  in  order  to  implement  two  cyclostationary 
processing  techniques  (Time  and  Frequency-Smoothing  algorithms).  The  Time 
Smoothing  FFT  Accumulation  Method  and  the  Direct  Frequency  Smoothing  Method 
were  implemented  in  MATLAB®  Version  6.1  [4]  and  tested  in  simulations  using  a 
standardized  data  set.  During  this  work  it  was  possible  to  verify  that  the  cyclostationary 
processing  approach  is  very  suitable  for  analyzing  radar  signals  with  LPI  characteristics. 
Depending  on  the  type  of  modulation,  the  processing  and  feature  extraction  can  be  more 
or  less  difficult.  The  analysis  of  LPI  signals  can  be  improved  for  any  type  of  signals  with 
the  combination  of  the  algorithms  implemented  here  and  other  techniques  such  as: 

•  Parallel  Filter  Arrays  and  Higher  Order  Statistics  [5]; 

•  Quadrature  Mirror  Filtering  [6];  and 

•  Wigner-Ville  Distribution  [7]. 

The  similarities  between  the  DFSM  results  and  the  FAM  results  go  until  a  certain 
level;  in  the  zoomed  plots  we  see  that  the  channel  pair  regions  are  a  little  different  in 
shape  and  size,  although  they  occur  in  the  same  values  for  frequency  and  cycle 
frequency.  This  may  be  the  result  of  the  different  windows  applied  in  each  method 
(Hamming  window  for  FAM  and  square  window  for  DFSM).  The  computing  time  is  also 
noticeable  larger,  two  or  three  times  more,  for  the  DFSM  routine  in  comparison  to  the 
FAM  routine.  The  FAM  implementation  is  recommended  for  long  signals  with  a  large 
number  of  samples. 

Both  methods  were  very  effective  for  phase  modulated  signals  (BPSK,  Frank,  PI, 
P2,  P3,  P4  and  FSK/PSK  Costas)  and  for  FH  signals  (Costas)  until  signal-to-noise  ratios 
of  OdB  and  -6dB.  The  FMCW  results  showed  more  variation  and  errors  in  the  results, 
comparing  to  the  original  characteristics  of  the  signals,  even  on  the  signal-only  cases.  The 
FSK/PSK  Target  signal  was  also  very  difficult  to  analyze  but  mainly  due  to  its  pseudo¬ 
random  characteristic  and  complexity.  In  the  technical  report  [10],  signals  with  different 
bandwidths  and  different  code  rates  were  included.  It  is  possible  to  see  that  the 
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cyclostationary  methods  are  very  effeetive  for  signals  with  larger  bandwidth  and  larger 
eode  rates.  In  order  to  improve  the  eomputation  time,  a  parallel  proeessing 
implementation  is  reeommended  for  future  work,  using  clusters  or  PC-based  parallel 
proeessors  for  academic  purposes. 

Starting  from  the  extraeted  features,  we  reeommend  a  future  researeh  on  the 
automatie  elassifieation  of  these  types  of  signals  as  a  form  of  maintaining  the  eontinuity 
of  this  work.  Neural  networks  have  being  exhaustively  pointed  as  a  viable  solution  for 
aoeomplishing  this  task.  The  use  of  MATLAB®  Version  6.1  and  GUIDE  Version  2.0  [4] 
was  very  important  to  the  sueeess  of  this  research  and  helped  the  development  of  a  very 
useful  tool  for  future  academic  applications  in  Electronic  Warfare  and  Signal  Processing 
oriented  courses  of  the  Naval  Post  Graduate  Sehool. 
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APPENDIX  A.  CYCLOSTATIONARY  IMPLEMENTATION 
CODES  (CYCLO.M,  FAM.M  AND  DFSM.M) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  CYCLO.M 

% 

%  Use:  GUI  implementation  for  inputing  parameters  in  Cyelostationary  Signal  Proeessing  Methods: 
DFSM.m  and  FAM.m 

% 

%  Inputs:  Path  of  the  file  to  be  analyzed;  Name  of  the  fde;  Sampling  Frequeney  used  in  the  generation  of 
the  fde  to  be  analyzed, 

%  desired  frequeney  resolution  and  desired  method  to  be  used  (DFSM  or  FAM) 

% 

%  Output:  Various  different  plots  exploring  the  bi-frequeney  plane  representation  typieal  of  eyelostationary 
proeessing 

% 

%  Created  by  Antonio  F.  Lima,  Jr.-lo  Ten.  Av.  Brazilian  Air  Foree  -  Mareh/2002 
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  These  are  default  eommands  and  eomments  from  the  MATLAB®  GUI  guide  tool 
%  DO  NOT  EDIT  ! ! 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

funetion  varargout  =  eyelo(varargin) 

%  CYCLO  Applieation  M-fde  for  eyelo.fig 
%  FIG  =  CYCLO  launeh  eyelo  GUI. 

%  CYCLO('eallbaek_name',  ...)  invoke  the  named  eallbaek. 

%  Last  Modified  by  GUIDE  v2.0  Ol-May-2002  12:56:05 

if  nargin  ==  0  %  LAUNCH  GUI 

fig  =  openfig(mfilename,'reuse'); 

%  Use  system  eolor  seheme  for  figure: 
set(fig,'Color',get(0,'defaultUieontrolBaekgroundColor')); 

%  Generate  a  strueture  of  handles  to  pass  to  eallbaeks,  and  store  it. 
handles  =  guihandles(fig); 
guidata(fig,  handles); 

if  nargout  >  0 

varargout{l}  =  fig; 
end 

elseif  isehar(varargin{l})  %  INVOKE  NAMED  SUBFUNCTION  OR  CALLBACK 
try 

if  (nargout) 

[varargout {Lnargout}]  =  feval(varargin{:});  %  FEVAL  switehyard 

else 

feval(varargin{:});  %  FEVAL  switehyard 
end 

eateh 
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end 


disp(lasterr); 


end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  These  are  the  GUI  functions  definitions  and  global  variables 
% - 


function  varargout  =  edit3_Callback(h,  eventdata,  handles,  varargin) 
d  =  get(h, 'String'); 

handles. edits  =  d;  %  Save  file  path  in  global  variable  handles. editS  (path  =  handles.editS) 
guidata(h, handles) ; 


% - 

function  varargout  =  editl_Callback(h,  eventdata,  handles,  varargin) 
name  =  get(h,'String'); 

handles. editl  =  name;  %  Save  file  name  in  global  variable  handles. editl  (filename  =  handles. editl) 

if  isstr(name) 
else 

errordlg('Filename:  Please  enter  a  valid  file  name  (B_l_7_7_s  is  default)',  'Bad  input',  'modal') 
end 

guidata(h,handles) ; 


% - 

function  varargout  =  edit2_Callback(h,  eventdata,  handles,  varargin) 
fs  =  str2double(get(h,'String')); 

handles. edit2  =  fs;%  Save  sampling  frequency  in  global  variable  handles. edit2  (fs  =  handles. edit2) 
if  isnan(fs) 

errordlg('Sampling  Frequency:  Please  enter  a  numeric  value  (Hz)',  'Bad  input',  'modal') 
end 

guidata(h,handles) ; 


% - 

function  varargout  =  popupmenul_Callback(h,  eventdata,  handles,  varargin) 
val  =  get(h, 'Value'); 

switch  val  %  Switches  between  the  values  chosen  by  user  on  popupmenu 


case  1 

handles.popupmenul  =  128;%  Save  frequency  resolution  in  global  variable  handles.popupmenuS  (df  = 
handles  .popupmenuS ) 

case  2 

handles.popupmenul  =  64;%  Save  frequency  resolution  in  global  variable  handles.popupmenuS  (df  = 
handles  .popupmenuS ) 

case  S 

handles.popupmenul  =  S2;%  Save  frequency  resolution  in  global  variable  handles.popupmenuS  (df  = 
handles  .popupmenuS ) 
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case  4 

handles. popupmenul  =  16;%  Save  frequency  resolution  in  global  variable  handles.popupmenu3  (df  = 
handles  .popupmenu3 ) 

case  5 

handles.popupmenul  =  8;%  Save  frequency  resolution  in  global  variable  handles. popupmenu3  (df  = 
handles  .popupmenu3 ) 

% 

case  6 

handles.popupmenul  =  4;%  Save  frequency  resolution  in  global  variable  handles. popupmenu3  (df  = 
handles  .popupmenu3 ) 


case  7 

handles.popupmenul  =  2;%  Save  frequency  resolution  in  global  variable  handles. popupmenu3  (df  = 
handles  .popupmenu3 ) 


end 


%  More  options  may  be  added  following  the  same  pattern. 
guidata(h, handles); 


% - 

function  varargout  =  popupmenu2_Callback(h,  eventdata,  handles,  varargin) 

Mval  =  get(h, 'Value'); 

switch  Mval  %  Switches  between  the  values  chosen  by  user  on  popupmenu 
case  1 

handles.popupmenu2  =  2;%  Save  reliability  condition  M  in  global  variable  handles.popupmenu2  (M  = 
handles  .popupmenu2) 

case  2 

handles.popupmenu2  =  4;%  Save  reliability  condition  M  in  global  variable  handles.popupmenu2  (M  = 
handles  .popupmenu2) 

case  3 

handles.popupmenu2  =  8;%  Save  reliability  condition  M  in  global  variable  handles.popupmenu2  (M  = 
handles. popupmenu2) 

end 

guidata(h, handles); 


%% - 

function  varargout  =  togglebuttonl_Callback(h,  eventdata,  handles,  varargin) 

dfsm;%  Runs  frequency  smoothing  method  (DFSM.m),  using  the  values  entered  -  Direct  Frequency 
Smoothing  Method 


% - 

function  varargout  =  togglebutton2_Callback(h,  eventdata,  handles,  varargin) 

fam;%  Runs  time  smoothing  method  (FFTAM.m),  using  the  values  entered  -  FFT  Accumulation  Method 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
%  FAM.M 

%  Cyclostationary  Estimator  using  the  time  smoothing  FFT  Accumulation  Algorithm. 

%  Part  of  the  LPI  oriented  thesis  work  with  Prof  P.E.  PACE  and  Prof  Ft.  LOOMIS 
%  Student:  Antonio  F.  LIMA  Jr.  -  1st.  Lt.  Brazilian  Air  Force 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 


%  Getting  input  data  from  the  global  variables  generated  by  GUI  cyclo  s.m  and  calculating  constants 
%  present  directory  =  cd;  %  Save  present  directory  to  a  variable  to  make  the  GUI  go  back  to  it  in  the  end 
path  =  handles. edits ;  %  Saving  the  path  information  in  a  variable 
cd  (path);%  Loading  the  fde  from  the  specified  path 

filename  =  handles. editl;  %  File  to  be  analyzed  from  GUI  generated  global  variable  handles. editl 
fs  =  handle s.edit2;  %  Sampling  frequency  from  GUI  generated  global  variable  handles. edit2 
df  =  handles.popupmenul;  %  Frequency  Resolution  from  GUI  generated  global  variable 
handles  .popupmenuS 

M  =  2;  %  df/dalpha  =  M  »  1  (reliability  condition) 
dalpha  =  df/M; 

load  (filename);  %loading  file  to  be  analyzed  in  workspace 

X  =  1+ j.*Q;%  using  complex  signal  (for  Real  part  only,  use  only  I  -  input  data  file  needs  to  be  formatted  as 
%  described  in  (4) 

%  Defining  parameters 

Np  =  pow2(nextpow2(fs/df));  %  Initially,  it  is  the  number  of  rows  in  the  channelization  matrix  (X)  or 
%  the  number  of  input  channels  defined  by  fs/df 

L  =  Np/4;  %  Overlap  factor  in  order  to  reduce  the  number  of  short  time  fft's 

%  L  is  the  offset  between  points  in  the  same  column  at  consecutive  rows 
%  should  be  less  than  or  equal  to  Np/4  (Prof.  Loomis  paper) 

P  =  pow2(nextpow2(fs/dalpha/L));  %  Number  of  columns  formed  in  the  channelization  matrix  (X) 

N  =  P*L;  %  Total  number  of  points  in  the  input  data  to  be  processed 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Input  channelization  -  this  part  limits  the  input  total  number  of  points  to  be  analyzed.  It  also  generates  a 
Np  X  P  matrix 

%  (X)  with  shifted  versions  of  the  input  vector  in  each  column. 

if  length(x)<N 
x(N)  =  0; 

elseif  length  (x)>N 
x  =  x(l:N); 
end 


NN  =  (P-1)*L+Np; 

XX  =  x; 

xx(NN)  =  0; 

XX  =  xx(:); 

X  =  zeros(Np,P); 
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for  k  =  0:P-1 

X(:,k+1)  =  xx(k*L+l:k*L+Np);  %  X  =  (Np  X  P)  Channelization  input  matrix 
end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Windowing  (Hamming  window)  -  a  vector  of  length  Np  is  created  with  the  MATLAB®  'hamming' 
function,  then  this  vector  is 

%  inserted  in  a  Np  X  Np  matrix  diagonal  and  this  result  is  multiplied  by  the  channelized  input  matrix  (X). 

a  =  hamming  (Np); 

XW  =  diag(a)*X; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  First  FFT  -  this  is  the  first  short  time  fft  (Np  points  long)  applied  to  each  column  of  matrix  X. 

XFl  =  fft(XW); 

XFl  =  fftshift(XFl); 

XFl  =  [XF1(:, P/2+1  :P)  XFl(:,l:P/2)]; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%Downconversion  -  the  short  sliding  fffs  results  are  shifted  to  base  band  to  obtain  decimated  complex 
demodulate  sequences  [2] 

E  =  zeros(Np,P); 
for  k  =  -Np/2:Np/2-l 
for  m  =  0:P-1 

E(k+Np/2+l,  m+1)  =  exp(-i*2*pi*k*m*L/Np); 
end 
end 

XD  =  XF1.*E; 

XD  =  conj(XD');  %  Transposing  the  matrix  and  taking  the  complex  conjugate  of  the  signal 
clear  ('XFl',  'E',  'XW',  'X',  'x','!',  'Q');  %  cleaning  up  memory  space 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Multiplication  -  the  product  sequences  between  each  one  of  the  complex  demodulates  and  the  complex 
conjugate  of  the  others 

%  are  formed.  This  forms  the  area  in  the  bi-frequency  plane  as  explained  in  [2]  figures  4 , 7  and  10 

XM  =  zeros(P,Np^2); 
for  k  =  1  :Np 
for  c  =  1  :Np 

XM(:,(k-l)*Np+c)  =  (XD(:,k).*conj(XD(:,c))); 
end 
end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Second  FFT  -  a  P  point  FFT  is  applied  to  XM  (in  each  of  its  columns) 

XF2  =  fft(XM); 

XF2  =  fftshift(XF2); 

XF2  =  [XF2(:,Np^2/2+l:Np^2)  XF2(:,l:Np^2/2)]; 

XF2  =  XF2  (P/4:3*P/4,:); 

MM  =  abs(XF2); 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 
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%Output  matrices  range 

alphaO  =  -fs:fs/N:fs; 
fD  =  -fs/2:fs/Np:fs/2; 

Sx  =  zeros(Np+l,  2*N+1); 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Display  only  the  data  inside  the  range  of  interest  -  centralizes  the  bi-frequency  plane  according  to  alphaO 
aird  ID  vectors 

forkl  =  l:?/2+l 
fork2=  l:Np^2 
if  rem(k2,Np)  ==  0 
c  =  Np/2-  1; 
else 

c  =  rem(k2,Np)  -  Np/2  -  1; 
end 

k  =  ceil(k2/Np)-Np/2-l; 
p  =  kl -P/4-1; 
alpha  =  (k-c)/Np-l-(p-l)/N; 
f  =  (k-l-c)/2/Np; 
ifalpha<-l  |  alpha>l 
k2=k2-H; 
elseif  f<-.5  I  f>.5 
k2=k2-H; 
else 

kk=  RNp*(f-K.5); 

11  =  l-l-N*(alpha  -I-  1); 

Sx(rouird(kk),  rouird(ll))  =  M]Vl(kl,k2); 
eird 
end 
end 

clear  ('alpha',  'XM',  'XF2',  'MM',  'f);  %  cleaning  up  memory  space 

Sx  =  Sx./max(max(Sx));  %  Normalizes  the  magnitudes  of  the  values  in  output  matrix  (maximum  =  1) 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

% 

%  Plotting  Routines 
%  saves  fde  containing  results 
%  savefile  =  'FFTAM  result  s'; 

%  save(savefde,'fdename','N','alphaO',  'fO',  'Sx',  'M',  'df) 

%  clear  all; 

%  clc; 

% 

%  load  ('FFTAM  result  s.mat'); 

%  figure  (df  -I-  M); 

%  mesh 

%  view  (-55,70); 

%  xlabel  ('cycle  frequency');  ylabel('frequency');  zlabel  ('SCD  magnitude'); 

%  title  (['Time  Smoothing  (FFTAM)  generated  SCD  for  ',  filename,  ',  df  =  ',  int2str(df),  '  and  N  =  ', 
int2str(N)]); 

% 

%  %  figure 

%  %  mesh  (alpha0,f0,  Sx,  'EdgeColor',  'interp'); 
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%  %  view  (0,90); 

%  %  xlabel  ('cycle  frequency');  ylabel('ffequency');  zlabel  ('SCD  magnitude'); 

%  %  title  (['SCD  for  ',  filename, '  N=',  int2str(N)]); 

%  %  % 

%  figure 

%  subplot(2,2,l:2); 

%  plot  (fO,  Sx);  grid; 

%  title  (['Time  Smoothing  Processing  for  ',  filename, ',  df  =  ',  int2str(df), '  and  N  =  ',  int2str(N)]); 

%  xlabel  ('frequency  (Hz)');  ylabel('SCD  magnitude'); 

%  subplot(2,2,3:4); 

%  plot  (alphaO,  Sx);  grid; 

%  xlabel  ('cycle  frequency  (Hz)');  ylabel  ('SCD  magnitude'); 

%  %  figure 
%%subplot(2,2,l:2); 

%  %  contour  (alphaO,  10,  Sx);  grid; 

%  %  xlabel('cycle  frequency  (Hz)');  ylabel('frequency  (Hz)'); 

%  %  title  (['Time  Smoothing  (FFTAM)  generated  SCD  for  ',  filename,  ',  df  =  ',  int2str(df),  '  and  N 
int2str(N)]); 

%  %  zoom(2);colorbar; 

%  %  subplot(2,2,3:4); 

%  %  contour  (alphaO,  10,  Sx);  grid; 

%  %  xlabel('cycle  frequency  (Hz)');  ylabel('frequency  (Hz)'); 

%  %  zoom(2);colorbar; 

% 

figure 

contour  (alphaO,  10,  Sx);  grid; 

xlabel('cycle  frequency  (Hz)');  ylabel('frequency  (Hz)'); 

title  (['Time  Smoothing  (FFTAM)  generated  SCD  for  ',  filename,  ',  df  =  ',  int2str(dl),  '  and  N 

int2str(N)]); 

colorbar; 

delete  ('FFTAM  result  s.maf); 

%  newname  =  ['FFTAM  '  filename]; 

%  cd('C:\Documents  and  Settings\aflima\Tese\Result  Figures'); 

%  saveas  (gcf,  newname,  'jpg'); 

cd('C:\Documents  and  Settings\aflima\Tese\cyclo'); 
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%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  ***  Cyclostationary  Estimator  using  the  Direct  Frequency  Smoothing  Method  algorithm.*** 

%  Part  of  the  LPI  oriented  thesis  work  with  Prof  P.E.  PACE  and  Prof  Ft.  LOOMIS 
%  Student:  Antonio  F.  LIMA  Jr.  -  1st.  Lt.  Brazilian  Air  Force 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Getting  input  data  from  the  global  variables  generated  by  GUI  cyclo  s.m 
%  function  f  =  dfsm  s  (x); 

%  present  directory  =  cd;  %  Save  present  directory  to  a  variable  to  make  the  GUI  go  back  to  it  in  the  end 

path  =  handles. edits ;  %  Saving  the  path  information  in  a  variable 

cd  (path);%  Loading  the  fde  from  the  specified  path 

filename  =  handle s. edit  1;  %  File  to  be  analyzed 

Fs  =  handle s.edit2;  %  Sampling  frequency 

df  =  handles.popupmenul;  %  Frequency  Resolution 

M  =  2;  %  Smoothing  window  range 

load  (filename); 

X  =  I;%  Using  real  part  of  signal 
N  =  (M*Fs)/df; 

N  =  pow2  (nextpow2(N));  %  windowing  record  for  FFT 


X  =  fft(x,N);  %  fft  of  the  truncated  (or  zero  padded)  time  series 
X  =  fftshift(X);%  shift  components  of  fft 

Xc  =  conj(X);  %  precompute  the  complex  conjugate  vector 


S  =  zeros  (N,N); 
f  =  zeros  (N,N); 
alfa  =  zeros  (N,N); 
F  =  Fs/(2*N); 

G  =  Fs/N; 
m  = -M/2+1  :M/2; 


%  size  of  the  Spectral  Correlation  Density  matrix 
%  size  of  the  frequency  matrix; 

%  size  of  the  cycle  frequency  matrix 
%  precompute  constants  -  F  =  Fs/(2*N); 

%  precompute  constants  -  G  =  Fs/N; 

%  set  frequency  smoothing  window  index 


for  k  =  1  :N  %  fix  k 

%  computes  vectors  of  f  and  alfa, 

%  store  frequency  and  cycle  frequency  data  for  given  k. 


kl  =  1:N; 

f(k,kl)  =  F*(k+kl-l)  -  Fs/2;  %  Computes  f  values  and  shift  them  to  center  in  zero  (f  =  (K+L)/2N) 

[1] 

alfa(k,kl)  =  G*(k-kl  +  N-1)  -  Fs;  %  Computes  alfa  values  and  shift  them  to  center  in  zero  (alfa  =  (K- 
L)/N)  [1] 


for  kl  =  1  :N  %fix  kl  =  J 

%calculate  X(K+m)  &  conj  (X(J+m))  for  arguments  of  X(1:N)  only 
B  =  max(l-k,  1-kl);  %  Largest  min  of  1  <=  (K+m)|  (J+m)  <=  N 

A  =  min  (N-k,  N-kl);  %  Smallest  max  of  1  <=  (K+m)|  (J+m)  <=  N 
n  =  m((m<=A)  &  (m>=B));  %fix  the  index  out  of  range  problem  by 

%  truncating  the  window 

if  isempty(n) 

S(k,kl)  =  0; 
else 

p  =  k+n;  q  =  kl+n; 

Y  =  X(p).*Xc(q); 
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S(k,kl)  =  sum(Y); 
end 
end 
end 

S  =  abs(S./max(max(S)));%  normalize  output  matrix 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  Saving  variables,  cleaning  up  memory  and  running  plots 

savefde  =  'DFSMresults'; 

%save(savefde,'r,  'Q',  'alfa',  'f,  'S',  'N',  'Fs',  'fdename',  'M',  'df) 
save(savefde,'alfa',  'f,  'S',  'N','fdename',  'M',  'df) 
clear  all; 

load  ('DFSM  result  s.maf); 

%  figure  (2*(df+M)); 

%  mesh  (alfa,f,  S,  'EdgeColor',  'interp'); 

%  view  (-55,70); 

%  xlabel  ('cycle  frequency');  ylabel('frequency');  zlabel  ('SCD  magnitude'); 

%  title  (['Frequency  Smoothing  (DFSM)  generated  SCD  for  ',  filename,  '  df  =  ',  int2str(df),  '  N  =  ', 
int2str(N)]); 

%  % 

%  figure  (2*(df+M)  -I-  1); 

%  mesh  (alfa,f,  S,  'EdgeColor',  'interp'); 

%  view  (0,90); 

%  xlabel  ('cycle  frequency');  ylabel('frequency');  zlabel  ('SCD  magnitude'); 

%  title  (['SCD  for',  filename, '  df  =',  int2str(df), '  N=',  int2str(N)]); 

% 

%  figure(2*(df+M)  -I-  2) 

%%  subplot  (2,2, 1:2) 

%  psd(I,  N,  Fs); 

%  title(['I  -  real  part  for  nfft  =  ',  int2str(N), '  and  for  default  nfft  =  256']); 

%  subplot  (2,2, 3:4) 

%psd(I,  [],Fs); 

%  subplot  (2,2,3) 

%  psd(I-t-j.*Q,  N/2,  Fs); 

%  title('I-tj.Q  -  complex'); 

%  subplot  (2,2,4) 

%  psd(I-j.*Q,  N/2,  Fs); 

%  title('I-j.Q  -  complex'); 

%  % 

%  figure 

%  subplot(2,2,l:2); 

%  plot  (f,  S);  grid; 

%  title  (['Frequency  Smoothing  Processing  for  ',  filename, ',  df  =  ',  int2str(df), '  and  N  =  ',  int2str(N)]); 

%  xlabel  ('frequency  (Flz)');  ylabel('SCD  magnitude'); 

%  subplot(2,2,3:4); 

%  plot  (alfa,  S);  grid; 

%  xlabel  ('cycle  frequency  (Flz)');  ylabel  ('SCD  magnitude'); 

%  %  % 

%  figure  (2*(df+M)  -I-  5) 

%  subplot(2,2,l:2); 

%  contour  (alfa,  f,  S);  grid; 

%  xlabel('cycle  frequency  (Flz)');  ylabel('frequency  (Flz)'); 
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%  title  (['Frequency  Smoothing  (DFSM)  generated  SCD  for  filename,  df  =  int2str(df),  '  and  N 
int2str(N)]); 

%  zoom(2);colorbar; 

%  subplot(2,2,3:4); 

%  contour  (alfa,  f,  S);  grid; 

%  xlabel('cycle  frequency  (Hz)');  ylabel('frequency  (Hz)'); 

%  zoom(2);colorbar; 
figure 

contour  (alfa,  f,  S);  grid; 

xlabel('cycle  frequency  (Hz)');  ylabel('frequency  (Hz)'); 

title  (['Frequency  Smoothing  (DFSM)  generated  SCD  for  ',  filename,  ',  df  =  ',  int2str(df),  '  and  N 

int2str(N)]); 

colorbar; 

delete  ('DFSM  result  s.mat'); 


newname  =  ['DFSM  '  filename]; 

cd('C:\Documents  and  Settings\afbma\Tese\Result  Figures'); 
saveas  (gcf,  newname,  'jpg'); 

cd('C:\Documents  and  Settings\aflima\Tese\cyclo'); 


136 


APPENDIX  B.  FSK/PSK  GENERATION  CODES 


%  FREQUENCY  AND  PHASE  SHIFT  COMBINED  CODE  using  Costas  Frequency  modulation 
%  Developed  by  Antonio  Lima  and  Harsha  Tummala 
%  July,  2002 

%  Phase  code  for  FSK/PSK  from  IEEE  International  Radar  Conference  Paper  -  1990 
%  “The  ambiguity  properties  of  FSK/PSK  signals” 

%  by  J.  Patrick  Donohoe,  and  Franklin  M.  Ingels 


clear  all; 
clc; 

disp('FREQUENCY  AND  PHASE  SHIFT  CODE  (FSK/PSK)  USING  COSTAS  FREQUENCIES******') 


dispC 


') 


%  DEFAULT  VARIABLES 


A=l; 

fs=15e3; 

SNRdB  =  0; 

cpf=l; 

scale=5; 

barker  =  5; 

j=sqrt(-l); 


%  Amplitude  of  CW 
%  Sample  Rate 

%  Signal  to  Noise  Ratio 
%Cycles  per  frequency  (>  10) 

%  Scaling  for  plotting  time  domain  graphs 
%  Number  of  Barker  Bots  for  phase  modulation 

%j 


%  NEW  INPUT 
newvar  =  1 ; 
while  newvar  ==  1 ; 
dispC  ') 

disp('WHICH  PARAMETER  DO  YOU  WANT  TO  SET  ?  ') 
disp(' ') 

fprintf('l.  Amplitude  of  frequencies  -  A=  %g.\n',  A) 
tprintf('2.  Sampling  frequency  -  fs  (Hz)=  %g.\n',  fs) 
t^rintf('3.  Signal  to  noise  ratio  -  SNR  dB  (dB)  =  %g.\n',  SNR  dB) 
tprintf('4.  Cycles  per  frequency  =  %gAn',  cpf) 

tprintf('5.  Number  of  bits  per  Barker  code  for  phase  modulation  -  barker  (13/1 1/7/5)=  %g.\n',  barker) 
tprintf('6.  No  changes\n') 
dispC  ') 

option=  input('Select  a  option:  '); 

switch  option 
case  1 

A=input('New  amplitude  of  the  carrier  signaU  '); 
case  2 

fs=input('New  sampling  frequency  (Hz)=  '); 
case  3 

SNR_dB=input('New  signal  to  noise  ratio  (dB)=  '); 
case  4 

cpf=input('New  number  of  cycles  per  frequency  =  '); 
case  5 

barker=input('New  number  of  bits  for  Barker  Code  =  '); 
case  6 

newvar  =  0; 
end 
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clc; 

end 

%  FREQUENCY  CHOICES 
newvar  =  1 ; 
while  newvar  ==  1 ; 
dispC  ') 

disp('WHICH  COSTAS  FREQUENCY  SEQUENCE  WOULD  YOU  LIKE  TO  USE  ?  ') 
disp(' ') 

disp(T.  4,  7,  1,  6,  5,  2,  3 '); 

disp('2.  2,  6,  3,  8,  7,  5,  1  '); 

disp('3.  4,  1,  6,  7,  5,  3,  2  '); 

disp('4.  5,  6,  8,  2,  4,  3,  1  '); 

disp('5.  4,  9,  2,  5,  1,9,8'); 
disp('6.  2,4,1,1,7,11,9'); 
disp('7.  2,  2,  5,  1,3,  7,4'); 
dispC ') 

option2=  input('Select  an  option:  '); 

freq=[4  7  1  6  5  2  3; 

2  6  3  8  7  5  1; 

4  1  6  7  5  3  2; 

5  6  8  24  3  1; 

492  5  1  9  8; 

24117119; 

22  5  1  3  7  4]*1000 

switch  option2 
case  1 

seq=freq(l,:); 
case  2 

seq=freq(2,:); 
case  3 

seq=freq(3,:); 
case  4 

seq=freq(4,:); 
case  5 

seq=freq(5,:); 
case  6 

seq=freq(6,:); 
case  7 

seq=freq(7,:); 


end 

newvar=0; 

clc; 

end 

%  Create  spanning  Barker  frequency 
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minimum=min(seq) ; 

SAR=ceil(fs/minimum);  %  Sampling  ratio 
tb=l/(fs);  %  Sampling  period 


%  This  section  generates  I  &  Q  without  COSTAS  phase  shift  and  I  &  Q  with  Phase  shift.  The  signals  are 
generated 

%  five  times  hy  the  outer  loop.  The  variable  'index'  is  used  to  generate  a  time  vector  for  time  domain  plots. 
%  The  signal  is  generated  at  seven  samples  per  phase  change. 

pw=  1; 

if  barker==13 

phase  =  [ones(l,pw*5),-(ones(l,pw*2)),ones(l,pw*2),-ones(l,pw),ones(l,pw),- 

ones(l,pw),ones(l,pw)];%  13  bits 
elseif  barker==l  1 

phase  =  [ones(l,pw*3),-(ones(l,pw*3)),ones(l,pw),-ones(l,pw),-ones(l,pw),ones(l,pw),-ones(l,pw)];% 
11  bits 
elseif  barker==7 

phase  =  [ones(l,pw*3),-(ones(l,pw*2)),ones(l,pw),-ones(l,pw)];%  7  bits 
elseif  barker==5 

phase  =  [ones(l,pw*3),-(ones(l,pw)),ones(l,pw)];%  5  bits 
end 

index=0; 
for  pp  =  1:5 

for  xx=l:7 

for  nn  =  1  :barker 

for  n=l:SAR*cpf 

I W  0(index+ 1  )=A*cos(2  *pi  *  seq(xx)  *  (n- 1 )  *  tb) ; 

QWO(index+l)=A*sin(2*pi*seq(xx)*(n-l)*tb); 

I(index+ 1  )=A*cos(2  *pi  *  seq(xx)  *  (n- 1 )  *  tb)  *phase(nn) ; 

Q(index+ 1  )=A*  sin(2  *pi*  seq(xx)  *  (n- 1 )  *tb)  *phase(nn) ; 
time(index+l)=index*tb;  %time  vector  cumulation 
index  =  index  +1; 
end 
end 
end 
end 

%Power  Spectral  Density  for  I  with  phase  shift  &  with  WGN  with  Signal  to  noise  ratios  (SNR)  =  [0,- 
5,5,10,-10,-20] 
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%for  loop  makes  calculations  and  plots  for  each  value  of  SNR  for  WGN 
[a,b]=size(I); 

SNR=10^(SNR_dB/10); 

power=10*logl0(A^2/(2*SNR));%calculate  SNR  indB  for  WGN  function 
noise=wgn(a,b,power);%calculate  noise  at  specified  SNR 
IN=I+noise;  %add  noise  to  I  with  FSK/PSK  phase  shift 

IPWON=I;  %I  with  phase  shift  without  noise 

QN=Q+noise;  %add  noise  to  Q  with  FSK/PSK  phase  shift 

QPWON=Q;  %Q  with  phase  shift  without  noise 

ffs  =  fs/1000; 

%PLOTS 

dispC ') 

pit  =  input('Do  you  want  to  generate  plots  of  the  signal  (Y/y  or  N/n)  ?Vs'); 
dispC  ') 

if(pit=='Y')|(pit==y) 
dispC ') 

%Plot  Power  Spectral  Density  for  I  without  phase  shift 
%  figure  ;  %  open  new  figure  for  plot 

%  psd(IWO,[],fs);  %Power  Spectral  Density  of  I  without  Phase  shift 
%  title(['PSD  of  I  without  Phase  Shift  or  Noise']); 


%Plot  PSD  of  1+  FSK/PSK  Phase  +  WGN  and  Time  Domain  of  I  +  FSK/PSK  Phase 
figure  ;%  open  new  figure  for  plot 
signal  =  I+j*Q; 

psd(signal,[],fs);%plot  PSD  for  specified  noise  SNR 

title(['PSD  of  FSK-PSK-C-'  num2str(option2)  num2str(ffs)  num2str(barker)  num2str(cpf)  '-s 

i+j*Q']); 

figure  ;%  open  new  figure  for  plot 
signal  =  IN+j*QN; 

psd(signal,[],fs);%plot  PSD  for  specified  noise  SNR 

title(['PSD  of  FSK-PSK-C-'  num2str(option2)  '-'  num2str(ffs)  '-'  num2str(barker)  '-'  num2str(cpf)  '-' 
num2str(SNR_dB) '  I+j*Q']); 

%  %plot  time  domain  signal  I  with  FSK/PSK  phase  shift  and  WGN  at  specified  SNR 
%  figure  ;%open  new  figure  for  plot 

%  plot(time(l  :floor(size(time,2)/scale)),I(l  :floor(size(time,2)/scale))); 

%  title(['FSK-PSK-C-'  num2str(option2) '-'  num2str(ffs) '-'  num2str(barker)  '-s  Time  Domain']); 

%  xlabel('{\itTime  -  Seconds}  '); 

%  ylabel('Amplitude'); 

%  grid  on; 

% 

%  Now  check  to  see  if  signal  is  correct  by  plotting  phase  shift  alone  and  then  determining  phase  shift 
from  I+jQ. 

%  To  determine  phase  shift,  look  at  the  phase  angle  of  I+iO  at  every  7th  time  interval.  Expect  to  see  the 
FSK/PSK  phase 

%  function  plot  repeated  5  times  after  unwrapping  and  detrending  the  phase  angle. 

figure;%open  new  figure  for  plot 
plot(phase); 
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title(['FSK-PSK-C-'  num2str(option2)  num2str(ffs)  num2str(barker)  num2str(cpf)  '-s  Phase  Shift 
for  each  Costas  frequency  hop']); 
xlabel('i  -  index  for  phase  change'); 
ylabel('FSK/PSK  Phase  Shift  -  Theta'); 
grid  on; 
else 

disp('Signal  not  plotted') 

^rintf('\n\n') 

end 

%  This  section  generates  the  fdes  for  analysis 

INP=IN';%transpose  I  with  noise  and  FSK/PSK  phase  shift  for  text  fde 
QNP=QN';%transpose  Q  with  noise  and  FSK/PSK  phase  shift  for  text  fde 
IPWONT=IPWON';%transpose  I  with  phase  without  noise  for  text  file 
QPWONT=QPWON';%transpose  Q  with  phase  without  noise  for  text  fde 

%  %  save  results  in  data  fdes 

I=INP(:,1); 

Q=QNP(:,1); 

11=  IPWONT(:,l); 

QQ=QPWONT(:,l); 

dispC ') 

saveresult  =  input('Do  you  want  to  save  the  new  signal  (Y/y  or  N/n)  ?','s'); 

if  (saveresult  ==  'Y')  |  (saveresult  =='y') 
ffs=floor(fs/le3); 

save(['FSK_PSK_C_'  num2str(option2)  '_'  num2str(ffs)  '_'  num2str(barker)  '_'  num2str(cpf)  '_' 
num2str(SNR_dB)],'r,'Q'); 

I=II; 

Q=QQ; 

save(['FSK_PSK_C_'  num2str(option2) '_'  num2str(ffs) '_'  num2str(barker) '_'  num2str(cpf)  '_s'],T,'Q'); 
dispC  '); 

disp(['Signal  and  noise  save  as  :  FSK  PSK  C  '  num2str(option2) '_'  num2str(ffs) '_'  num2str(barker) '_' 
num2str(cpf) '_'  num2str(SNR_dB)]); 

disp(['Signal  only  save  as  :  FSK  PSK  C  '  num2str(option2)  '_'  num2str(ffs)  '_'  num2str(barker)  '_' 

num2str(cpf)  '_s']); 

disp(['Directory:  '  num2str(cd)]); 

else 
dispC ') 

disp('Signal  not  saved') 

^rintf('\n\n') 

end 
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%  FREQUENCY  AND  PHASE  SHIFT  CODE  (FSK/PSK)  for  a  target  frequency  distribution 
%  Developed  by  Antonio  Lima  and  Harsha  Tummala 
%  July,  2002 

%  Phase  code  for  FSK/PSK  from  IEEE  International  Radar  Conference  Paper 
%  MATCHED  FSK/PSK  RADAR 

%  by  B.  Jeffrey  Skinner,  J.  Patrick  Donohoe,  and  Franklin  M.  Ingels 


clear  all; 
clc; 

disp(  *********************************************************************** ^^ 
jjgp(,H=************pREQUENCY  AND  PHASE  SHIFT  CODE  (PSK/PSK)*************'); 
disp(  *********************************************************************** ^^ 


%DEFAULT  VARIABLES 


A=l; 

fs  =15000; 
SNRdB  =  0; 
scale=20; 
j=sqrt(-l); 
global  N; 
N=128; 
cpp  =  5; 


%  Amplitude  of  CW 
%  Sampling  Frequency 
%  Signal  to  Noise  Ratio 
%  Scaling  for  plotting  time  domain  graphs 

%j 

%  Number  of  fsk/psk  sections  (frequency  hops  and  phase  changes) 
%  Number  of  cycles  per  phase 


%  NEW  INPUT 
newvar  =  1 ; 
while  newvar  ==  1 ; 
dispC  ') 

disp('WHICH  PARAMETER  DO  YOU  WANT  TO  SET  ?  ') 
dispC ') 

fprintf(T.  Amplitude  of  the  carrier  signal  -  A=  %g.\n'.  A) 
fprintf('2.  Sampling  frequency  -  fs  (Hz)=  %g.\n',  fs) 
fprintf('3.  Signal  to  noise  ratio  -  SNR  dB  (dB)  =  %g.\n',  SNR  dB) 
fprintf('4.  Number  of  phase/frequency  hops  -  N  =  %g.\n',  N) 
fprintf('5.  Number  of  cycles  per  phase  -  cpp  =  %g.\n',  cpp) 
fprintf('6.  No  changes\n') 
dispC  ') 

option=  input('Select  a  option:  '); 


switch  option 
case  1 

A=input('New  amplitude  of  the  carrier  signaU  '); 
case  2 

fs=input('New  sampling  frequency  (Hz)=  '); 
case  3 

SNR_dB=input('New  signal  to  noise  ratio  (dB)=  '); 
case  4 

N=input('New  number  of  phase/frequency  hops  ='); 
case  5 

cpp=input('New  number  of  cycles  per  phase='); 
case  6 

newvar  =  0; 
end 
clc; 
end 


airplane;  %  Compute  frequency  hops  values  according  to  an  airplane's  frequency  response 
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syn  test;  %  Generates  frequeney  distribution  aceording  to  previous  probability  distribution 


global  target;  %  Compute  Frequency,  first  call  up  airplane  variable 
global  detection;  % 

tb=l/(fs);  %  Sampling  period 

SAR=ceil(fs/min(detection));  %  Sampling  ratio  (for  smaller  frequency) 

phase  =  rand(l,cpp);  %Compute  the  phase  encoding  random  sequence  for  each  frequency  burst 

greater  =  fmd(phase>=0.5); 

phase(greater)  =  pi; 

lesser  =  fmd(phase<0.5); 

phase(lesser)  =  0; 


%  This  section  generates  I  &  Q  without  phase  shift  and  I  &  Q  with  Phase  shift.  The  signals  are  generated 
%  five  times  the  number  of  frequency  hops  by  the  outer  loop.  The  variable  'index'  is  used  to  generate  a 
time  vector  for  time  domain  plots. 

%  The  signal  is  generated  at  seven  samples  per  phase  change. 
index=0;  %  Time  vector  for  time  domain  plots. 

for  p=l  :2%  Generate  the  signal  two  times  and  store  sequentially  in  corresponding  vectors 
for  ii  =  (1  :N)  %Loop  for  each  frequency  hop 
f  =  detection(ii,:); 
for  kk  =  (l:size(target,l)) 
if  f  ==  target(kk,l) 

init_phase  =  target(kk,3);  %defmes  initial  phase  of  FSK  sequence 
end 
end 

for  n=l:SAR*cpp  %Loop  to  increment  time  for  single  frequency  value. 

I(index+l)=A*cos(2*pi*P(n-l)*tb+phase(ceil(n/SAR))+init_phase);  %Calculate  in  phase 
component  of  signal  with  phase  shift 

IWO(index+l)=A*cos(2*pi*P(n-l)*tb);  %  Calculate  in  phase  component  of  signal  without  phase 
shift 

Q(index+l)=A*sin(2*pi*P(n-l)*tb+phase(ceil(n/SAR))+init_phase);  %  Calculate  quadrature 
component  of  signal  with  phase  shift 

QWO(index+l)=A*sin(2*pi*P(n-l)*tb);  %Calculate  quadrature  component  of  signal  without 
phase  shift 

time(index+l)=index*tb;  %time  vector  cumulation 
index  =  index  +1; 
end 
end 
end 


%Power  Spectral  Density  for  I  with  phase  shift  &  with  WGN  with  Signal  to  noise  ratios  (SNR)  =  [0,- 
5,5,10,-10,-20] 

%for  loop  makes  calculations  and  plots  for  each  value  of  SNR  for  WGN 
[a,b]=size(I); 

SNR=10^(SNR_dB/10); 

power=10*logl0(A^2/(2*SNR));%calculate  SNR  indB  for  WGN  function 
noise=wgn(a,b,power);%calculate  noise  at  specified  SNR 
IN=H-noise;  %add  noise  to  I  with  FSK/PSK  phase  shift 

IPWON=I;  %I  with  phase  shift  without  noise 

QN=Q-t-noise;  %add  noise  to  Q  with  FSK/PSK  phase  shift 
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QPWON=Q;  %Q  with  phase  shift  without  noise 

ffs  =  fs/1000; 

%PLOTS 

dispC ') 

pit  =  input('Do  you  want  to  generate  plots  of  the  signal  (Y/y  or  N/n)  ?Vs'); 
dispC  ') 

if(pit=='Y')|(pit==y) 
dispC ') 
global  range; 
global  airplane; 

%Plot  Target  Frequeney  Distribution  and  Range/Magnitude  plot 
figure;%  open  new  figure  for  plot 
plot(targetC,l),  target(:,2));  title('Ship  FFT  ABS');  grid  on 
title(['Original  Target  Frequeney  Probability  Distribution']); 
xlabel('Frequeney'); 

ylabel('Normalized  Magnitude  =  Probability'); 

figure;%  open  new  figure  for  plot 

plot(range,  real(airplane(l:64)));  title('Ship');  grid  on 

title(['Original  Target  Range  /  Magnitude  Plot']); 

xlabel('Range  (ft)'); 

ylabel('Magnitude'); 

%  Plot  original  frequeney  distribution  histogram  and  frequeney  random  firing  distribution 

figure;%  open  new  figure  for  plot 

orient  tall; 

subplot(2,l,l), 

hist(deteetion( : ,  1  ),N) ; 

xlabel('Deteetion  Index'); 

ylabel('Number  of  Oeeurenees'); 

fidl=['Target  SYNTHETIC']; 

title(fidl); 

subplot(2,l,2), 
bar(target(:,l),  target(:,2)); 
xlabel('Deteetion  Index') 
ylabel('Probabibty') 
fid2=['Target  ORIGINAL']; 
title(fid2); 

%Plot  Power  Speetral  Density  for  I  without  phase  shift 
%  figure  ;  %  open  new  figure  for  plot 

%  psd(IWO,[],fs);  %Power  Speetral  Density  of  I  without  Phase  shift 
%  title(['PSD  of  I  without  Phase  Shift  or  Noise']); 

%time  domain  plot  of  in  phase  signal  I  with  phase  shift 
%  figure  ;  %open  new  figure  for  plot 

%  %  plot  small  portion  of  time  domain  signal  I  so  that  data  will  fit  meaningfully  in  figure. 

%  %floor(size(time,2)/seale)  seleets  a  small  sample  of  the  veetors  to  plot 
%  plot  (time(l  :floor(size(time,2)/seale)),I(l  :floor(size(time,2)/seale))); 

%  title(['Time  Domain  of  I  with  Phase  Shift  &  no  Noise']); 
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%  xlabel('{\itTime  -  Seconds}  '); 

%  ylabel('Amplitude'); 

%  grid  on; 

% 

%Plot  PSD  of  1+  FSK/PSK  Phase  +  WGN  and  Time  Domain  of  I  +  FSK/PSK  Phase 
figure  ;%  open  new  figure  for  plot 
signal  =  I+j*Q; 

psd(signal,[],fs);%plot  PSD  for  specified  noise  SNR 

title(['PSD  of  FSK-PSK-T-'  num2str(ffs) num2str(N) num2str(cpp)  '-s  I+j*Q']); 

%plot  time  domain  signal  I  with  FSK/PSK  phase  shift  and  WGN  at  specified  SNR 
figure  ;%open  new  figure  for  plot 

plot(time(l:floor(size(time,2)/scale)),I(l:floor(size(time,2)/scale))); 
title(['FSK-PSK-T-'  num2str(ffs) num2str(N) num2str(cpp)  '-s  Time  Domain']); 
xlabelC  {\itTime  -  Seconds}  '); 
ylabel('Amplitude'); 
grid  on; 

%  Now  check  to  see  if  signal  is  correct  by  plotting  phase  shift  alone  and  then  determining  phase  shift 
from  I+jQ. 

%  To  determine  phase  shift,  look  at  the  phase  angle  of  I+jQ  at  every  7th  time  interval.  Expect  to  see  the 
FSK/PSK  phase 

%  function  plot  repeated  5  times  after  unwrapping  and  detrending  the  phase  angle. 

figure;%open  new  figure  for  plot 
plot(phase); 

title(['FSK-PSK-T-'  num2str(ffs) num2sfr(N) num2sfr(cpp)  '-s  Phase  Shift ']); 
xlabel('i  -  index  for  phase  change'); 
ylabel('FSK/PSK  Phase  Shift  -  Theta'); 
grid  on; 
else 

disp('Signal  not  plotted') 

^rintf('\n\n') 

end 

%  This  section  generates  the  files  for  analysis 

INP=IN';%transpose  I  with  noise  and  FSK/PSK  phase  shift  for  text  file 
QNP=QN';%transpose  Q  with  noise  and  FSK/PSK  phase  shift  for  text  file 
IPWONT=IPWON';%franspose  I  with  phase  without  noise  for  text  file 
QPWONT=QPWON';%transpose  Q  with  phase  without  noise  for  text  file 

%  %  save  results  in  data  files 

I  =  INP(:,1); 

Q  =  QNP(:,1); 

II  =  IPWONT(:,l); 

QQ  =  QPWONT(:,l); 

dispC ') 

saveresult  =  input('Do  you  want  to  save  the  new  signal  (Y/y  or  N/n)  ?','s'); 

if  (saveresult  ==  'Y')  |  (saveresult  =='y') 
ffs=floor(fs/le3); 
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save(['FSK_PSK_T_'  num2str(ffs)  num2str(N)  num2str(cpp)  num2str(SNR_dB)],'r,'Q', 
'detection'); 

I=II; 

Q=QQ; 

save(['FSK_PSK_T_'  num2str(ffs) '_'  num2str(N) '_'  num2str(cpp)  '  s'], T/Q', 'detection'); 
dispC '); 

disp(['Signal  and  noise  save  as  :  FSK  PSK  T  '  num2str(ffs)  '_'  num2str(N)  '_'  num2str(cpp)  '_' 
num2str(SNR_dB)]); 

disp(['Signal  only  save  as  :  FSK  PSK  T  '  num2str(ffs) '_'  num2str(N) '_'  num2str(cpp)  '_s']); 

disp(['Directory:  '  num2str(cd)]); 

else 
dispC ') 

disp('Signal  not  saved') 

^rintf('\n\n') 

end 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  This  part  may  be  implemented  on  a  different  fde 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

function  airplane  =  airplane; 

global  target; 
global  airplane; 
global  range; 

%  Target  range  Profile 

airplane  =  [complex(0, -0.00000032:0.00000001  :-0.00000016),  complex(-0. 00 1,-0. 000000 15),  complex(- 
0.002,-0.00000014),  complex(-0.003,-0.00000013),  complex(-0.0027,-0.00000012),  complex(- 
0.0025,-0.00000011), ... 

complex(-0. 0022, -0.00000010),  complex(-0.0020,  -0.00000009),  complex(-0.0017,  -0.00000008), 
complex(-0.00 15,-0.00000007),  complex(-0.0023,  -0.00000006),  complex(-0.003,  -0.00000005), 
complex(-0.0015,  -0.00000004), ... 

complex(0, -0.00000003),  complex(0,-0. 00000002),  complex(-0.0015,  -0.00000001),  complex(- 
0. 003,0), complex(0.028,  0.00000001),  complex(-0.0075,  0.00000002),  complex(0.014, 

0.00000003),  complex(0.0 13,0.00000004),  complex(-0.038,0.00000005), ... 

complex(-0.003,  0.00000006),  complex(-0.004,  0.00000007),  complex(-0.005,  0.00000008), 
complex(-0.015,  0.00000009),  complex(0.003,  0.0000001),  complex(-0.0015, 

0.000000 1 1 :0.0000000 1 :0.000000 1 6),  complex(0,  0.000000 1 7:0.0000000 1 :0.0000003 1), 

complex(0,-0.00000032:0.00000001:0.00000031)]; 

%  Target  Time  pro  fde 
airplane  t  =  flipdim(airplane,  2); 

airplane  t  =  real(airplane_t).*(-l)  -I- j.*(imag(airplane_t)); 

%  FFT  of  the  time  record  fo  the  target  response 
airplanel  =  fft(airplane_t); 
airplane_2  =  abs(fftshift(airplane_l)); 
phairplane  =  phase(airplane_l); 

freq_step  =  64;  %  Number  of  frequency  steps  to  be  used 

fss  =  12e6; 

fs  =  fss/2000;  %Downsized  sampling  frequency  (from  12MF1Z  to  6000Ftz) 
frequency  =  fs/2:fs/2/freq_step:fs  -  1; 


146 


airplanes  =  airplane_2(64:127); 
total  =  sum(airplane_3); 
airplane_4  =  airplane_3/total; 

target  =  [(frequency)',  (airplane_4)',  (pli_airplane(l:freq_step))']; 

range_bin  =  (fss*2e-9)^(-l); 

range  =  l:range_bin:frec|_step*range_bin-l; 


%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  This  part  may  be  implemented  on  a  different  fde 

%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%%  FREQUENCY  SYNTHESIS  ALGORITHM  FOR  FSK/PSK  TARGET 
%%  Frequencies  are  in  the  vector  "detection" 
function  syn  test  =  syn  test; 

fid=['Ship']; 
global  N; 
global  target; 
new_rec=  target(:,l:2); 

%NO  EDITS  BELOW 

%  PULL  OFF  THE  DECIMAL  DIGITS 

[nn,jo]=size(new_rec);  %  Determine  number 
%  of  detections  or  rows 

%  Pull  off  the  density 
%  decimal  digits  as  strings 
“delta”=sprintf('%6.4f,new_rec(:,2)); 
ss=reshape(“delta”,  6,  nn); 
ss=ss'; 

%  Turn  back  into  a  number 
%  with  each  digit  a  column 
“delta”s(:,[l])=str2num(ss(:,3)); 

“delta”s(:,[2])=str2num(ss(:,4)); 

“delta”s(:,[3])=str2num(ss(:,5)); 

“delta”s( : ,  [4])=str2num(ss(: ,  6)) ; 

%%  kk  =4;  %  NUMBER  OF  DECIMAL  DIGITS 

%  Define  probabilities 
%  pi  through  p4 

p0=0; 

pl=10^(-l)  *  sum(“delta”s(:,l)); 
p2=10^(-2)  *  sum(“delta”s0,2)); 
p3=10^(-3)  *  sum(“delta”s0,3)); 
p4=10^(-4)  *  sum(“delta”s0,4)); 

%  Determine  pi  values 
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piO=  0; 
pil=  10*pl; 

pi2=  (10*pl)+(10^2  *  p2); 

pi3=  (10*pl)+(10^2  *  p2)+(10^3  *  p3); 

pi4=  (10*pl)+(10^2  *  p2)+(10^3  *  p3)  +  (lOM  *  p4); 

%  Determine  test  values  for 
%  the  uniform  random  variable 

ptest0=0; 

ptestl=pl; 

ptest2=ptestl+p2; 

ptest3=ptest2+p3; 

ptest4=ptest3+p4; 

%  FILL  MEMORY  LOCATION 
%  SET  1 

ok=0; 

for  j=l  :nn  %  RANDOM  VARIABLE  INDEX 
if“delta”sG,l)~=0 

fori=l:”delta”sG,l) 

memG+ok,l)=j; 

end 

ok=i+ok; 

end 

end 

if“delta”sG,l)~=0 

[mem  size  xx]=size(mem); 
else 

memsize  =  0; 

XX  =  0; 

end 


%  FILL  MEMORY  LOCATIONS 
%  SET  2 

ok=0; 

for  j=l  :nn  %  RANDOM  VARIABLE  INDEX 
if“delta”sG,2)~=0 

for  i=l:”delta”sG,2) 

mem(i+ok+mem_size,  1  )=j ; 
end 

ok=i+ok; 

end 

end 

[mem  size  xx]=size(mem); 

%  FILL  MEMORY  LOCATIONS 
%  SET  3 

ok=0; 

for  j=l  :nn  %  RANDOM  VARIABLE  INDEX 
if“delta”sG,3)~=0 

for  i=l:”delta”sG,3) 

mem(i+ok+mem_size,  1  )=j ; 
end 

ok=i+ok; 
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end 


end 


[mem  size  xx]=size(mem); 

%  FILL  MEMORY  LOCATIONS 
%  SET  4 

ok=0; 

for  j=I  :nn  %  RANDOM  VARIABLE  INDEX 
if  “delta”s(j,4)~=0 

for  i=I:”delta”s(j,4) 

mem(i+ok+mem_size,  I  )=j ; 
end 

ok=i+ok; 

end 

end 

% _ 

%  Now  that  mem  is  fdled 
%  generate  the  detections 

%  nn  number  of  detections 

for  gi=I:N 

%  uni  is  the  uniform  RV 
%  pull  the  decimal  digits  off 

uni=rand; 

uni_str=sprintf('%6.4f,uni)  ; 
sss=reshape(uni_str,  6,  1); 
sss=sss'; 

dl=str2num(sss(:,3)); 
d2=str2num(sss( :  ,4)) ; 
d3 =str2num(sss( : ,  5)) ; 
d4=str2num(sss( : ,  6)) ; 

%  Test  the  RV  to  find  correct 
%  index.  Then  generate  detection 

%  DETECTION 

global  detection; 

if  uni  >=  0  &  uni  <  ptestl 

rv_index=mem(round(d  1 + 1 ),  1 ) ; 
detection(gi,  [  1  ])=new_rec(rv_index,  (1:1)); 
elseif  uni  >=  ptestl  &  uni  <  ptest2 
const=pi  l-(100*(pl)); 

rv_index=mem(round(dl*10  +  d2  +  const  +1),1); 
detection(gi,  [  1  ])=new_rec(rv_index,  (1:1)); 
elseif  uni  >=  ptest2  &  uni  <  ptest3 
const=pi2-(  1 000 *  (p  1  +p2)) ; 

rv_index=mem(round(dl*100  +  d2*10  +  d3  +  const  +1),1); 
detection(gi,  [  1  ])=new_rec(rv_index,  (1:1)); 
elseif  uni  >=  ptest3  &  uni  <  ptest4 

const=pi3 -( 1 0000 *  (p  1  +p2+p3)) ; 

rv_index=mem(round(dl*1000  +  d2*100  +  d3*10  +d4  +  const  +1),1); 
detection(gi,  [  1  ])=new_rec(rv_index,  (1:1)); 
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end 


%  % _ PLOTS 

%  % 

%  figure 

% 

%  orient  tall; 

%  subplot(2,l,l), 

%  liist(detection(:,l),N); 

%  xlabel('Detection  Index'); 

%  ylabel('Number  of  Occurences'); 

%  fidl=[fid,'  SYNTHETIC']; 

%  title(fidl); 

% 

%  subplot(2,I,2), 

%  bar(new_rec(:,I),  new_rec(:,2)); 

%  xIabel('Detection  Index') 

%  ylabel('Probability') 

%  fid2=[fid,'  ORIGINAL']; 

%  title(fid2); 
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%PAF  for  FREQUENCY  AND  PHASE  SHIFT  CODE  (FSK/PSK) 

%  Developed  by  Antonio  Lima 
%  July,  2002 

%  Phase  code  for  FSK/PSK  from  IEEE  International  Radar  Conference  Paper 
%  “MATCHED  FSK/PSK  RADAR” 

%  by  B.  Jeffrey  Skinner,  J.  Patrick  Donohoe,  and  Franklin  M.  Ingels 

clear  all; 
clc;  close  all; 

load  FSK_PSK_T_15_256_10_s;  %  Loads  the  I  and  Q  for  the  P4  signal,  1  cycle  N^2  =  128 

cpp  =  10; 

NN  =  256; 


%  load  FSK_PSK_C_l_15_5_l_s;  %  Loads  the  I  and  Q  for  the  P4  signal,  1  cycle  N^2  =  128 

%  seq  =  1 ; 

%  barker  =  5; 

%cpf=  1; 

SignalI=I(l:size(I,l)/2)';  %  Signal  repeates  after  896  intervals,  128  phase  codes*l  samples/phase*3 

SAR*5 

SignalQ=Q(l:size(Q,l)/2)'; 

SignalQJ=SignalQ.*j;  %  Make  Q  the  complex  part 
ul=SignalI+SignalQ  j;  %  Create  signal  and  reference  signal 
u2=SignalI-  SignalQ  J ; 


tan  =  1 ; 

T=l; 

N=l; 

tx  =  length(u2); 
tb  =  T/tx; 
vstep  =  500; 
t=[0:tx-l]*tau/tx; 
dt  =  t(2)-t(l); 


%  Adjust  number  of  pulses  in  the  train 
%  Number  of  sambles  in  the  time  domain 
%  Code(chip)  period 

%  Steps  in  the  freq  shift  dimension  -  (600  FOR  COSTAS) 
%  Steps  in  the  time  dimension,  positive 


z=[];  %  Empty  result  matrix 

v=linspace(-50/tau,50/tau,vstep);  %  Width  in  the  freq  shift  dimension  -  (TARGET,  50  FOR  COSTAS) 
for  m=  1  :vstep 
mult=[]; 

mult  =  exp(j*2*pi*v(m)*t); 
u2ex=  u2.*mult; 

u3  =  conj(u2ex);  %  Transfer  in  order  to  use  built  in  function 

c  =  xcorr(u3,ul).*dt;  %  Using  the  built  in  correlation  function 

universaU[];  %  Multiplying  with  the  universal  function 

dummy  =  (pi*v(m)*T); 

dummy  =  dummy  +  (dummy==0)*eps; 

universal  =  abs(sin(dummy*N)/(N*sin(dummy))); 

e=  (c.*  universal); 

z=[z;abs(e)]; 

v(m)  =  v(m)*tb;  %  Normalizing  the  freq  shift  axis 

end 

t=[fliplr(-t),t(2:tx)];  %  Creating  a  negative  time  axis 

t=(t/(tb*tx));  %  Normalizing  the  time  axis 

z  =  z/max(max(z)); 
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%  FOR  FSKPSKTARGET,  USE  THE  FOLLOWING 

PLOTS%%%%%%%%%%%%%%%%%%%%%% 

figure 

contour(t,v,z),grid  on 
colorbar 

title(['FSK/PSK  Target  Code,  Cpp  =  num2str(cpp) N  hops  =  '  num2str(NN)]) 
xlabel('Normalized  \tau') 
ylabel('Normalized  \nu') 

%  axis  tight 

figure 

plot(t,z((vstep+l)/2,:)),  grid  on 

title('PSK/FSK  Target  code,  Cut  along  the  0  doppler  axis') 
xlabel('Normahzed  \tau') 
ylabel('Magnitude') 
axis  tight 

figure 

plot(v,z(:,tx)),  grid  on 

title('PSK/FSK  Target  code.  Cut  along  the  0  delay  axis') 
xlabel('Normahzed  \nu') 
ylabel('Magnitude') 
axis  tight 


%  FOR  FSKPSKCOSTAS,  USE  THE  FOLLOWING 

PLOTS%%%%%%%%%%%%%%%%%%%%%%% 

%  figure 

%  contour(t,v,z),grid  on 
%  colorbar 

%  title(['PAF  for  a  FSK/PSK  Costas  Code,  Costas  sequence  =  ',  num2str(seq)  ',  Barker  Bits  =  ' 
num2str(barker) ',  Cycles  per  Phase  =  '  num2str(cpf)]) 

%  xlabel('Normahzed  \tau') 

%  ylabel('Normahzed  \nu') 

%  axis  tight 
% 

%  figure 

%  plot(t,z((vstep+l)/2,:)),  grid  on 

%  title('PSK/FSK  Costas  code.  Cut  along  the  0  doppler  axis') 

%  xlabel('Normahzed  \tau') 

%  ylabel('Magnitude') 

%  axis  tight 

% 

%  figure 

%  plot(v,z(:,tx)),  grid  on 

%  title('PSK/FSK  Costas  code.  Cut  along  the  0  delay  axis') 

%  xlabel('Normahzed  \nu') 

%  ylabel('Magnitude') 

%  axis  tight 

% 

% 

%  %Doing  it  again 
%  index  =  0; 

%  for  ii  =  l:(cpf*375):((size(I,l))/2) 

%  for  pp  =  0:74 
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%  Signall(index+1)  =  I(ii+pp)'; 

%  SignalQ(index+l)  =  Q(ii+pp)'; 

%  index  =  index  +  1 ; 

%  end 
%  end 
% 

%  SignalQJ=SignalQ.*j;  %  Make  Q  the  complex  part 
%  ul=SignalI+SignalQJ;  %  Create  signal  and  reference  signal 
%  u2=SignalI-SignalQJ; 

% 

%  tan  =  1 ; 

%T  =  1; 

%  N  =  1 ;  %  Adjust  number  of  pulses  in  the  train 

%  tx  =  length(u2);  %  Number  of  sambles  in  the  time  domain 
%  tb  =  T/tx;  %  Code(chip)  period 

%  vstep  =  500;  %  Steps  in  the  freq  shift  dimension  -  (600  FOR  COSTAS) 

%  t=[0:tx-l]*tau/tx;  %  Steps  in  the  time  dimension,  positive 
%dt  =  t(2)-t(l); 

%  z=[];  %  Empty  result  matrix 

%  v=linspace(-20/tau,20/tau, vstep);  %  Width  in  the  freq  shift  dimension  -  (TARGET,  50  FOR  COSTAS) 
%  for  m=  1  :vstep 


% 

mult=[]; 

% 

mult  =  exp(j*2*pi*v(m)*t); 

% 

u2ex=  u2.*mult; 

% 

u3  =  conj(u2ex); 

%  Transfer  in  order  to  use  built  in  function 

% 

c  =  xcorr(u3,ul).*dt; 

%  Using  the  built  in  correlation  function 

% 

universaU[]; 

%  Multiplying  with  the  universal  function 

%  dummy  =  (pi*v(m)*T); 

%  dummy  =  dummy  +  (dummy==0)*eps; 

%  universal  =  abs(sin(dummy*N)/(N*sin(dummy))); 

%  e=  (c.*  universal); 

%  z=[z;abs(e)]; 

%  v(m)  =  v(m)*tb;  %  Normalizing  the  freq  shift  axis 

%  end 

%  t=[fliplr(-t),t(2:tx)];  %  Creating  a  negative  time  axis 

%  t=(t/(tb*tx));  %  Normalizing  the  time  axis 

%  z  =  z/max(max(z)); 

% 

%  %%%%%%%%%%%%%%  FOR  FSK_PSK_COSTAS,  USE  THE  FOLLOWING 

PLOTS%%%%%%%%%%%%%%%%%%%%%%% 

%  figure 

%  contour(t,v,z),grid  on 
%  colorbar 

%  title(['PAF  for  a  FSK/PSK  Costas  Code,  Costas  sequence  =  ',  num2str(seq)  ',  Barker  Bits  =  ' 
num2str(barker) ',  Cycles  per  Phase  =  '  num2str(cpf)]) 

%  xlabel('Normahzed  \tau') 

%  ylabel('Normahzed  \nu') 

%  axis  tight 
% 

%  figure 

%  plot(t,z((vstep+l)/2,:)),  grid  on 

%  title('PSK/FSK  Costas  code.  Cut  along  the  0  doppler  axis') 

%  xlabel('Normahzed  \tau') 

%  ylabel('Magnitude') 

%  axis  tight 

% 
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%  figure 

%  plot(v,z(:,tx)),  grid  on 

%  title('PSK/FSK  Costas  code,  Cut  along  the  0  delay  axis') 

%  xlabel('Normalized  \nu') 

%  ylabel('Magnitude') 

%  axis  tight 

% 

%  %%%%%%%%%%%  FOR  OTHER  PLOTS,  USE  THE 

FIGURES%%%%%%%%%%%%%%%%%%%%%%%%%%% 

%  %  figure 

%  %  colormap(jet);  caxis=[3  1]; 

%  %  mesh(t,v,z),grid  on 
%  %  colorbar; 

%  %  title('PAF  for  a  FSK/PSK  Code') 

%  %  xlabel('Normalized  \tau') 

%  %  ylabel('Normalized  \nu') 

%  %  zlabel('Magnitude') 

%  %  axis  tight 
%  % 

% 

%  %  figure 

%  %  mesh(t,v,z),grid  on 
%  %  view(0,0) 

%  %  colorbar 

%  %  title('Side-view  towards  the  delay  axis') 

%  %  xlabel('Normalized  \tau') 

%  %  ylabel('Normalized  \nu') 

%  %  zlabel('Magnitude') 

%  %  axis  tight 
%  % 

%  %  figure 

%  %  mesh(t,v,z),grid  on 
%  %  view(90,0) 

%  %  colorbar 

%  %  title('Side-view  towards  the  doppler  axis') 

%  %  xlabel('Normalized  \tau') 

%  %  ylabel('Normalized  \nu') 

%  %  zlabel('Magnitude') 

%  %  axis  tight 
%  % 


FOLLOWING 


154 


APPENDIX  C.  LIST  OF  LPI  RADAR  SIGNALS  ANALYZED 


SIGNAL  NAME 

CARRIER 
FREQUENCY 
OR  HOPPING 
SEQUENCE 
(kHz) 

BW 

(Hz) 

NPBB/ 

CPP 

NUMBER  OF 
BITS  /  PHASES 
AND 

MODULATION 

PERIOD 

MODULATION  TYPE 

B  1  7  11  1  -6.mat 

1 

1000 

1 

11 

BPSK 

B  1  7  11  1  O.mat 

1 

1000 

1 

11 

BPSK 

B  1  7  11  1  s.mat 

1 

1000 

1 

11 

BPSK 

B  1  7  11  5  -6.mat 

1 

200 

5 

11 

BPSK 

B  1  7  11  5  O.mat 

1 

200 

5 

11 

BPSK 

B  1  7  11  5  s.mat 

1 

200 

5 

11 

BPSK 

B  1  7  7  1  -6.mat 

1 

1000 

1 

7 

BPSK 

B  1  7  7  1  O.mat 

1 

1000 

1 

7 

BPSK 

B  1  7  7  1  s.mat 

1 

1000 

1 

7 

BPSK 

B  1  7  7  5  -6. mat 

1 

200 

5 

7 

BPSK 

B  1  7  7  5  O.mat 

1 

200 

5 

7 

BPSK 

B  1  7  7  5  s.mat 

1 

200 

5 

7 

BPSK 

FR  1  7  16  1  -6.mat 

1 

1000 

1 

16 

FRANK 

FR  1  7  16  1  O.mat 

1 

1000 

1 

16 

FRANK 

FR  1  7  16  1  s.mat 

1 

1000 

1 

16 

FRANK 

FR  1  7  16  5  -6.mat 

1 

200 

5 

16 

FRANK 

FR  1  7  16  5  O.mat 

1 

200 

5 

16 

FRANK 

FR  1  7  16  5  s.mat 

1 

200 

5 

64 

FRANK 

FR  1  7  64  1  -6.mat 

1 

1000 

1 

64 

FRANK 

FR  1  7  64  1  O.mat 

1 

1000 

1 

64 

FRANK 

FR  1  7  64  1  s.mat 

1 

1000 

1 

64 

FRANK 

FR  1  7  64  5  -6.mat 

1 

200 

5 

64 

FRANK 

FR  1  7  64  5  O.mat 

1 

200 

5 

64 

FRANK 

FR  1  7  64  5  s.mat 

1 

200 

5 

64 

FRANK 

F  1  7  250  20  -6.mat 

1 

250 

X 

20 

FMCW 

F  1  7  250  20  O.mat 

1 

250 

X 

20 

FMCW 

F  1  7  250  20  s.mat 

1 

250 

X 

20 

FMCW 

F  1  7  250  30  -6.mat 

1 

250 

X 

30 

FMCW 

F  1  7  250  30  O.mat 

1 

250 

X 

30 

FMCW 

F  1  7  250  30  s.mat 

1 

250 

X 

30 

FMCW 

F  1  7  500  20  -6.mat 

1 

500 

X 

20 

FMCW 

F  1  7  500  20  O.mat 

1 

500 

X 

20 

FMCW 

F  1  7  500  20  s.mat 

1 

500 

X 

20 

FMCW 

F  1  7  500  30  -6.mat 

1 

500 

X 

30 

FMCW 

F  1  7  500  30  O.mat 

1 

500 

X 

30 

FMCW 

F  1  7  500  30  s.mat 

1 

500 

X 

30 

FMCW 

Table  49.  Test  matrix  of  LPI  radar  signals  analyzed. 
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PI  1  7  16  1  -6.mat 

1 

1000 

1 

16 

PI 

PI  1  7  16  1  O.mat 

1 

1000 

1 

16 

PI 

PI  1  7  16  1  s.mat 

1 

1000 

1 

16 

PI 

PI  1  7  16  5  -6.mat 

1 

200 

5 

16 

PI 

PI  1  7  16  5  O.mat 

1 

200 

5 

16 

PI 

PI  1  7  16  5  s.mat 

1 

200 

5 

16 

PI 

PI  1  7  64  1  -6.mat 

1 

1000 

1 

64 

PI 

PI  1  7  64  1  O.mat 

1 

1000 

1 

64 

PI 

PI  1  7  64  1  s.mat 

1 

1000 

1 

64 

PI 

PI  1  7  64  5  -6.mat 

1 

200 

5 

64 

PI 

PI  1  7  64  5  O.mat 

1 

200 

5 

64 

PI 

PI  1  7  64  5  s.mat 

1 

200 

5 

64 

PI 

P2  1  7  16  1  -6.mat 

1 

1000 

1 

16 

P2 

P2  1  7  16  1  O.mat 

1 

1000 

1 

16 

P2 

P2  1  7  16  1  s.mat 

1 

1000 

1 

16 

P2 

P2  1  7  16  5  -6.mat 

1 

200 

5 

16 

P2 

P2  1  7  16  5  O.mat 

1 

200 

5 

16 

P2 

P2  1  7  16  5  s.mat 

1 

200 

5 

16 

P2 

P2  1  7  64  1  -6.mat 

1 

1000 

1 

64 

P2 

P2  1  7  64  1  O.mat 

1 

1000 

1 

64 

P2 

P2  1  7  64  1  s.mat 

1 

1000 

1 

64 

P2 

P2  1  7  64  5  -6.mat 

1 

200 

5 

64 

P2 

P2  1  7  64  5  O.mat 

1 

200 

5 

64 

P2 

P2  1  7  64  5  s.mat 

1 

200 

5 

64 

P2 

P3  1  7  16  1  -6.mat 

1 

1000 

1 

16 

P3 

P3  1  7  16  1  O.mat 

1 

1000 

1 

16 

P3 

P3  1  7  16  1  s.mat 

1 

1000 

1 

16 

P3 

P3  1  7  16  5  -6.mat 

1 

200 

5 

16 

P3 

P3  1  7  16  5  O.mat 

1 

200 

5 

16 

P3 

P3  1  7  16  5  s.mat 

1 

200 

5 

16 

P3 

P3  1  7  64  1  -6.mat 

1 

1000 

1 

64 

P3 

P3  1  7  64  1  O.mat 

1 

1000 

1 

64 

P3 

P3  1  7  64  1  s.mat 

1 

1000 

1 

64 

P3 

P3  1  7  64  5  -6.mat 

1 

200 

5 

64 

P3 

P3  1  7  64  5  O.mat 

1 

200 

5 

64 

P3 

P3  1  7  64  5  s.mat 

1 

200 

5 

64 

P3 

P4  1  7  16  1  -6.mat 

1 

1000 

1 

16 

P4 

P4  1  7  16  1  O.mat 

1 

1000 

1 

16 

P4 

P4  1  7  16  1  s.mat 

1 

1000 

1 

16 

P4 

P4  1  7  16  5  -6.mat 

1 

200 

5 

16 

P4 

P4  1  7  16  5  O.mat 

1 

200 

5 

16 

P4 

P4  1  7  16  5  s.mat 

1 

200 

5 

16 

P4 

P4  1  7  64  1  -6.mat 

1 

1000 

1 

64 

P4 

P4  1  7  64  1  O.mat 

1 

1000 

1 

64 

P4 

Table  50.  Test  matrix  of  LPI  radar  signals  analyzed. 
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P4  I  7  64  I  s.mat 

1 

1000 

P4  I  7  64  5  -6.mat 

1 

200 

P4  I  7  64  5  O.mat 

1 

200 

P4  I  7  64  5  s.mat 

1 

200 

C  I  15  10  -6.mat 

4,  7,  1,  6,  5,  2,  3 

- 

C  I  15  10  O.mat 

4,  7,  1,  6,  5,  2,  3 

- 

C_l_15_I0_s.mat 

4,  7,  1,6,  5,  2,3 

- 

C  I  15  20  -6.mat 

4,  7,1,6,  5,2,3 

- 

C_I_I5_20_0.mat 

4,  7,  1,6,  5,2,3 

- 

C  I  15  20  s.mat 

4,  7,  1,  6,  5,  2,  3 

- 

C  2  17  10  -6.mat 

2,  6,  3,  8,  7,  5,  1 

- 

C  2  17  10  O.mat 

2,  6,  3,  8,  7,  5,  1 

- 

C  2  17  10  s.mat 

2,  6,  3,  8,  7,  5,  1 

- 

C  2  17  20  -6.mat 

2,  6,  3,  8,  7,  5,  1 

- 

C  2  17  20  O.mat 

2,  6,  3,  8,  7,  5,  1 

- 

C  2  17  20  s.mat 

2,  6,  3,  8,  7,  5,  1 

- 

FSK  PSK  C  I  15  II  I  O.mat 

4,  7,  1,  6,  5,  2,  3 

1000 

FSK  PSK  C  I  15  II  I  s.mat 

4,  7,  1,  6,  5,  2,  3 

1000 

FSK  PSK  C  I  15  II  5  O.mat 

4,  7,  1,6,  5,  2,  3 

200 

FSK  PSK  C  I  15  II  5  s.mat 

4,  7,  1,  6,  5,  2,  3 

200 

FSK  PSK  C  I  15  5  I  O.mat 

4,  7,  1,6,  5,2,3 

1000 

FSK  PSK  C  I  15  5  I  s.mat 

4,  7,  1,6,  5,2,3 

1000 

FSK  PSK  C  I  15  5  5  O.mat 

4,  7,  1,6,  5,2,3 

200 

FSK  PSK  Cl  15  5  5  s.mat 

4,  7,  1,  6,  5,  2,  3 

200 

FSK_PSK_T_I5_I28_I0_0.mat 

- 

3600 

FSK  PSK  T  15  128  10  s.mat 

- 

3600 

FSK  PSK  T  15  128  5  O.mat 

- 

4200 

FSK  PSK  T  15  128  5  s.mat 

- 

4200 

FSK  PSK  T  15  256  10  O.mat 

- 

3600 

FSK  PSK  T  15  256  10  s.mat 

- 

3600 

FSK_PSK_T_I5_256_5_0.mat 

- 

4200 

FSK  PSK  T  15  256  5  s.mat 

- 

4200 

T  I  7  I  s.mat 

1000 

T  12  7  2  s.mat 

1000  and  2000 

Table  51.  Test  matrix  of  LPI  radar  signals  ai 
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64 

P4 

64 

P4 

64 

P4 

64 

P4 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

1 

COSTAS 

11 

FSK/PSK  COSTAS 

11 

FSK/PSK  COSTAS 

11 

FSK/PSK  COSTAS 

11 

FSK/PSK  COSTAS 

5 

FSK/PSK  COSTAS 

5 

FSK/PSK  COSTAS 

5 

FSK/PSK  COSTAS 

5 

FSK/PSK  COSTAS 

128  Hops 

FSK/PSK  TARGET 

128  Hops 

FSK/PSK  TARGET 

128  Hops 

FSK/PSK  TARGET 

128  Hops 

FSK/PSK  TARGET 

256  Hops 

FSK/PSK  TARGET 

256  Hops 

FSK/PSK  TARGET 

256  Hops 

FSK/PSK  TARGET 

256  Hops 

FSK/PSK  TARGET 

TEST  SIGNAL  1 

- 

CARRIER 

TEST  SIGNAL  2 
CARRIERS 
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